Bosonic Spectral Function and The Electron-Phonon Interaction in HTSC Cuprates 
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In Part I we discuss accumulating experimental evidence related to the structure and origin of 
the bosonic spectral function a^F{uj) in high-temperature superconducting (HTSC) cuprates at and 
near optimal doping. Some global properties of a^F{uj), such as number and positions of peaks, are 
extracted by combining optics, neutron scattering, ARPES and tunnelling measurements. These 
methods give convincing evidence for strong electron-phonon interaction (EPI) with 1 < X^p < 3 in 
cuprates near optimal doping. Here we clarify how these results are in favor of the Eliashberg-like 
theory for HTSC cuprates near optimal doping. 

In Part II we discuss some theoretical ingredients - such as strong EPI, strong correlations - 
which are necessary to explain the experimental results related to the mechanism of d-wave pairing 
in optimally doped cuprates. These comprise the Migdal-Eliashberg theory for EPI in strongly 
correlated systems which give rise to the forward scattering peak. The latter is further supported 
by the weakly screened Madelung interaction in the ionic- metallic structure of layered cuprates. In 
this approach EPI is responsible for the strength of pairing while the residual Coulomb interaction 
(by including spin fluctuations) triggers the d-wave pairing. 



Part I 

Experimental Evidence 
for Strong EPI 

I. INTRODUCTION 

In spite of an unprecedented intensive experimen- 
tal and theoretical study after the discovery of high- 
temperature superconductivity (HTSC) in cuprates there 
is, even twenty-three years after, no consensus on the 
pairing mechanism in these materials. At present there 
are two important experimental facts which are not un- 
der dispute: (1) the critical temperature Tc in cuprates 
is high, with the maximum T™^'^ ~ 160 K in the Hg- 
1223 compounds; (2) the pairing in cuprates is d-wave 
like, i.e. A(k, cj) « As(fc,w) -|- Ad{u}){coskx — cosky) 
with As < O.lArf. On the contrary there is a dispute 
concerning the scattering mechanism which governs nor- 
mal state properties and pairing in cuprates. To this 
end, we stress that in the HTSC cuprates, a number 
of properties can be satisfactorily explained by assum- 
ing that the quasi-particle dynamics is governed by some 
electron-boson scattering and in the superconducting 
state bosonic quasi-particles are responsible for Cooper 
pairing. Which bosonic quasi-particles are dominating 
in the cuprates is the subject which will be discussed in 
this work. It is known that the electron-boson (phonon) 
scattering is well described by the Migdal-Eliashberg the- 
ory if the adiabatic parameter A = a ■ A(aJs/H4) fulfills 
the condition A <^ 1, where A is the electron-boson cou- 
pling constant, is the characteristic bosonic energy 
and Wb is the electronic band width and a depends on 
numerical approximations Jj]. The important character- 



istic of the electron-boson scattering is the Eliashberg 
spectral function a^FCk^k' ,uj) (or its average a'^F{ijj)) 
which characterizes scattering of quasi-particle from k 
to k' by exchanging bosonic energy w. Therefore, in sys- 
tems with electron-boson scattering the knowledge of this 
function is of crucial importance. There are at least two 
approaches differing in assumed pairing bosons in the 
HTSC cuprates. The first one is based on the electron- 
phonon interaction (EPI), with the main proponents in 
H, 0i Hi @; @; where mediating bosons are phonons 
and where the average spectral function a'^F{uj) is similar 
to the phonon density of states Fph{uj). Note, a'^F{uj) is 
not the product of two functions although sometimes one 
defines the function a'^{u}) = F {u)) / F {oj) which should 
approximate the energy dependence of the strength of 
the EPI coupling. There are numerous experimental evi- 
dence in cuprates for the importance of the EPI scatter- 
ing mechanism with a rather large coupling constant in 
the normal scattering channel 1 < Aep ^ 3, which will be 
discussed in detail below. In the EPI approach a'^Fph{'^) 
is extracted from tunnelling measurements in conjunction 
with IR optical measurements. We stress again that the 
Migdal-Eliashberg theory is well justified framework for 
EPI since in most superconductors the condition A < 1 
is fulfilled. The HTSC cuprates are on the borderline and 
it is a natural question - under which condition can high 
Tc be realized in the non- adiabatic limit A ^ 11 The 
second approach [tJ assumes that EPI is too weak to be 
responsible for high Tc in cuprates and it is based on a 
phenomenological model for spin-fluctuation interaction 
(SFI) as the dominating scattering mechanism, i.e. it 
is a non-phononic mechanism. In this (phenomenologi- 
cal) approach the spectral function is proportional to the 
imaginary part of the spin susceptibility /mx(k — k',a;), 
i.e. a2F(k, k',w) ^ /mx(k — k',aj). NMR spectroscopy 
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and magnetic neutron scattering give evidence that in 
HTSC cuprates x(q, w) is peaked at the antiferromag- 
netic wave vector Q = {Tr/a,Tr/a) and this property is 
very favorable for d-wave pairing. The SFI theory roots 
basically on the strong electronic repulsion on Cu atoms, 
which is usually studied by the Hubbard model or its 
(more popular) derivative the t-J model. Regarding the 
possibility to explain high Tc solely by strong correla- 
tions, as it is reviewed in we stress two facts. First, at 
present there is no viable theory as well as experimental 
facts which can justify these (non-phononic) mechanisms 
of pairing with some exotic pairing mechanism such as 
RVB pairing 0, fractional statistics and anyon super- 
conductivity, etc. Therefore we shall not discuss these, 
in theoretical sense, very interesting approaches. Second, 
the central question in these non-phononic approaches 
is - do models based solely on the Hubbard Hamilto- 
nian show up superconductivity at sufficiently high crit- 
ical temperatures {T^ ~ 100 K) 1 Although the answer 
on this important question is not definitely settled there 
are a number of numerical studies of these models which 
offer rather convincing negative answers. For instance, 
the sign-free variational Monte Carlo algorithm in the 
2D repulsive {U > 0) Hubbard model gives no evidence 
for superconductivity with high Tc, neither the BCS- nor 
Berezinskii-Kosterlitz-Thouless (BKT)-like |?|. At the 
same time, similar calculations show that there is a strong 
tendency to superconductivity in the attractive ([/ < 0) 
Hubbard model for the same strength of J7, i.e. at finite 
temperature in the 2D model with [/ < the BKT super- 
conducting transition is favored. Concerning the possi- 
bility of HTSC in the t — J model, various numerical cal- 
culations such as Monte Carlo calculations of the Drude 
spectral weight [To] and high temperature expansion for 
the pairing susceptibility have shown that there is 
no superconductivity at temperatures characteristic for 
cuprates and if it exists Tc must be rather low - few 
Kelvins. These numerical results tell us that the lack of 
high Tc (even in 2D BKT phase) in the repulsive [U > 0) 
single-band Hubbard model and in the t — J model is not 
only due to thermodynamical 2L)-fluctuations (which at 
finite T suppress and destroy superconducting phase co- 
herence in large systems) but it is mostly due to an inher- 
ent ineffectiveness of strong correlations to produce solely 
high Tc in cuprates. These numerical results signal that 
the simple single-band Hubbard and its derivative the t-J 
model are insufficient to explain solely the pairing mech- 
anism in cuprates and some additional ingredients must 
be included. 

Since EPI is rather strong in cuprates, then it must 
be accounted for. As it will be argued in the follow- 
ing, the experimental support for the importance of EPI 
in cuprates comes from optics, tunnelling, and recent 
ARPES measurements [l^. It is worth mentioning that 
recent ARPES activity was a strong impetus for re- 
newed experimental and theoretical studies of EPI in 
cuprates. However, in spite of accumulating experimen- 
tal evidence for importance of EPI with Xcp > 1, there 



are occasionally reports which doubt its importance in 
cuprates. This is the case with recent interpretation of 
some optical measurements in terms of SFI only (T^l , [13] , 
its'] and with LDA-DFT (local density approximation- 
density functional theory) band structure calculations 
[11, [ll, where both claim that EPI is negligibly small, 
i.e. Aep < 0.3. The inappropriateness of these calcula- 
tions will be discussed in the following Sections. 

The paper is organized as follows. In Part I we will 
mainly discuss experimental results in cuprates at and 
near optimal doping by giving also minimal theoreti- 
cal explanations which are related to the bosonic spec- 
tral function a'^F(uj) as well to the transport spectral 
function af^F{u!) and their relations to EPI. The rea- 
son that we study only cuprates at and near optimal 
doping is, that in these systems there are rather well 
defined quasi-particles - although strongly interacting, 
while in highly underdoped systems the superconduc- 
tivity is perplexed and possibly masked by other phe- 
nomena, such as pseudogap effects, formation of small 
polarons, interaction with spin and (possibly charge) or- 
der parameters, pronounced inhomogeneities of the scat- 
tering centers, etc. As ARPES experiments confirm 
there are no polaronic effects in systems at and near 
optimal doping, while there are pronounced polaronic 
effects due to EPI in undoped and very underdoped 
HTSC [1], [I]- In this work we consider mainly those 
direct one-particle and two-particles probes of low en- 
ergy quasi-particle excitations and scattering rates which 
give information on the structure of the spectral functions 
a^F(k, k',a;) and afj.F{u}) in systems near optimal dop- 
ing. These are angle-resolved photoemission (ARPES), 
various arts of tunnelling spectroscopy such as super- 
conductor/insulator/ normal metal (SIN) junctions and 
break junctions, scanning-tunnelling microscope spec- 
troscopy (STM), infrared (IR) and Raman optics, in- 
elastic neutron and x-ray scattering, etc. We shall argue 
that these direct probes give evidence for a rather strong 
EPI in cuprates. Some other experiments on EPI are 
also discussed in order to complete the arguments for the 
importance of EPI in cuprates. The detailed contents of 
Part I is the following. In Section II we discuss some prej- 
udices related to the strength of EPI as well as on the 
Fermi-liquid behavior of HTSC cuprates. We argue that 
any non-phononic mechanism of pairing should have very 
large bare critical temperature Tco S> Tc in the presence 
of the large EPI coupling constant, Aep > 1, if the EPI 
spectral function is weakly momentum dependent, i.e. 
if a^F{k., k', cj) « a'^F{uj) like in low temperature super- 
conductors. The fact that EPI is large in the normal state 
of cuprates and the condition that it must conform with 
d-wave pairing implies inevitably that EPI in cuprates 
must be strongly momentum dependent. In Section III 
we discuss direct and indirect experimental evidence for 
the importance of EPI in cuprates and for the weakness 
of SFI in cuprates. These are: 

(A) Magnetic neutron scattering measurements - 
These measurements provide dynamic spin susceptibil- 
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ity x(q, w) which is in the SFI phenomenological ap- 
proach related to the Eliashberg spectral function, 
i.e. a^Fs/(k, k', oj) ~ g'^flmxil = k — k',a;). We stress 
that such an approach can be theoretically justified only 
in the weak coupling limit, gsf Wb, where Wb is the 
band width and gsf is the phenomenological SFI coupling 
constant. Here we discuss experimental results on YBCO 
which give evidence for strong rearrangement (with re- 
spect to uj) of /mx(q, w) (with q at and near Q = (n^n)) 
by doping toward the optimal doped HTSC [ll], [l9| . 
It turns out that in the optimally doped cuprates with 
Tc — 92.5 K /mx(Q,i^) is drastically suppressed com- 
pared to that in slightly underdoped ones with Tc — 91 
K. This fact implies that the SFI coupling constant gsf 
must be small. 

(B) Optical conductivity measurements - From these 
measurements one can extract the transport relaxation 
rate 7f^(a;) and indirectly an approximative shape of 
the transport spectral function afj,F{uj). In the case 
of systems near optimal doping we discuss the follow- 
ing questions: (i) the physical and quantitative differ- 
ence between the optical relaxation rate 7j^(w) and the 
quasi-particle relaxation rate 7(0;). It was shown in 
the past that by equating these two (unequal) quanti- 
ties is dangerous and brings incorrect results concerning 
the quasi-particle dynamics in most metals by includ- 
ing HTSC cuprates too 0, 0, [l^, [Hi, [13, [111; (ii) 
methods of extraction of the transport spectral function 
af^F{uj). Although these methods give at finite temper- 
ature T a blurred aj^F(aj) which is (due to the ill-defined 
methods) temperature dependent, it turns out that the 
width and the shape of the extracted a^^F^uj) are in fa- 
vor of EPI; (iii) the restricted sum-rule for the opti- 
cal weight as a function of T which can be explained by 
strong EPI [13], [l^]; (iv) good agreement with experi- 
ments of the T-dependence of the resistivity p{T) in op- 
timally doped YBCO, where p(T) is calculated by using 
the spectral function from tunnelling experiments. Re- 
cent femtosecond time-resolved optical spectroscopy in 
La2-xSrxCu04^ which gives additional evidence for im- 
portance of EPI [26J will be shortly discussed. 

(C) ARPES measurements and EPI - From these 
measurements the self-energy E(k, w) is extracted as well 
as some properties of a^i^(k, k',a;). Here we discuss the 
following items: (i) existence of the nodal and anti-nodal 
kinks in optimally and slightly underdoped cuprates, as 
well as the structure of the ARPES self-energy (S(k, w)) 
and its isotope dependence, which are all due to EPI; 
(ii) appearance of different slopes of I](k, oj) at low {uj <C 
LOph) and high energies (w 3> ujph) which can be explained 
by strong EPI; (iii) formation of small polarons in the un- 
doped HTSC was interpreted to be due to strong EPI - 
this gives rise to phonon side bands which are clearly seen 
in ARPES of undoped HTSC @. 

(D) Tunnelling spectroscopy - It is well known that 
this method is of an immense importance in obtain- 
ing the spectral function a^i^(a;) from tunnelling con- 
ductance. In this part we discuss the following items: 



(i) the extracted Eliashberg spectral function a^F(cLi) 
with the coupling constant A*-*""'' = 2 — 3.5 from 
the tunnelling conductance of break-junctions in opti- 
mally doped YBCO and Bi-2212 which gives 
that the maxima of a'^F{u!) coincide with the max- 
ima in the phonon density of states Fph{ijj)] (ii) exis- 
tence of eleven peaks in —d^I/dV^ in superconducting 
Lai.s,ASrQ,iQCuOi films }34] . where these peaks match 
precisely with the peaks in the intensity of the exist- 
ing phonon Raman scattering data [35]; (iii) the pres- 
ence of the dip in dl/dV in STM which shows the pro- 
nounced oxygen isotope effect and important role of these 
phonons. 

{E) Inelastic neutron and x-ray scattering measure- 
ments - From these experiments one can extract the 
phonon density of state Fph{uj) and in some cases 
strengths of the quasi-particle coupling with various 
phonon modes. These experiments give sufficient evi- 
dence for quantitative inadequacy of LDA-DFT calcula- 
tions in HTSC cuprates. Here we argue, that the large 
softening and broadening of the half-breathing Cu — O 
bond-stretching phonon, of apical oxygen phonons and of 
oxygen B-^g buckling phonons (in LSCO, BISCO,YBCO) 
can not be explained by LDA-DFT. It is curious that the 
magnitude of the softening can be partially obtained by 
LDA-DFT but the calculated widths of some important 
modes are an order of magnitude smaller than the neu- 
tron scattering data show. This remarkable fact confirms 
additionally the inadequacy of LDA-DFT in strongly cor- 
related systems and a more sophisticated many body the- 
ory for EPI is needed. The problem of EPI will be dis- 
cussed in more details in Part II . 

In Section IV brief summary of the Part I is given. 
Since we are dealing with electron-boson scattering in 
cuprate near optimal doping, then in Appendix A (and 
in Part II) we introduce the reader briefiy into the 
Migdal-Eliashberg theory for superconductors (and nor- 
mal metals) where the quasi-particle spectral func- 
tion a^i<"(k, k', cj) and the transport spectral function 
a1j.F{ijj) are defined. 

Finally, one can pose a question - do the experimental 
results of the above enumerated spectroscopic methods 
allow a building of a satisfactory and physically reason- 
able microscopic theory for basic scattering and pairing 
mechanism in cuprates? The posed question is very mod- 
est compared with a much stringent request for the theory 
of everything - which would be able to explain all proper- 
ties of HTSC materials. Such an ambitious project is not 
realized even in those low-temperature conventional su- 
perconductors where it is definitely proved that in most 
materials the pairing is due to EPI and many properties 
are well accounted for by the Migdal-Eliashberg theory. 
Let us mention only two examples: - First, the exper- 
imental value for the coherence peak in the microwave 
response CTs(T < Tc,ui = const) at w = 17 GHz in the 
superconducting Nh is much higher than the theoretical 
value obtained by the strong coupling Eliashberg theory 
[36t . So to say, the theory explains the coherence peak 
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at 17 GHz in Nb qualitatively but not quantitatively. 
However, the measurements at higher frequency a; ~ 60 
GHz are in agreement with the Ehashberg theory [37l |. 
Then one can say that instead of the theory of every- 
thing we deal with a satisfactory theory, which allows 
us qualitative and in many aspects quantitative explana- 
tion of phenomena in superconducting state. - Second 
example is the experimental boron (B) isotope effect in 
MgB2 {Tc ~ 40 K) which is much smaller than the the- 
oretical value, i.e. « 0.3 < = 0.5, although 
the pairing is due solely by EPI for boron vibrations [38j . 
Since the theory of everything is impossible in the com- 
plex materials such as HTSC cuprates in Part I we shall 
not discuss those phenomena which need much more mi- 
croscopic details and/or more sophisticated many-body 
theory. These are selected by chance: (i) large ratio 
2A/Tc which is on optimally doped YBCO and BISCO 
~ 5 and 7, respectively, while in underdoped BISCO one 
has even (2A/Tc) ~ 20; [ii) peculiarities of the coherence 
peak in the microwave response <t(T) in HTSC cuprates, 
which is peaked at T much smaller than Tc, contrary to 
the case of LTSC where it occurs near Tc; {in) the de- 
pendence of Tc on the number of Cu02 in the unit cell; 
(w) temperature dependence of the Hall coefficient; [v) 
distribution of states in the vortex core, etc. 

The microscopic theory of the mechanism for super- 
conducting pairing in HTSC cuprates will be discussed 
in Part H. In Section V we introduce an ah initio many- 
body theory of superconductivity which is based on the 
fundamental (microscopic) Hamiltonian and the many- 
body technique. This theory can in principle calcu- 
late measurable properties of materials such as the crit- 
ical temperature Tc, critical fields, dynamic and trans- 
port properties, etc. However, although this method 
is in principle exact, which needs only some fundamen- 
tal constants e, fi-, mg, Mio„, fc^ and the chemical com- 
position of superconducting materials, it was practically 
never realized in practice due to the complexity of many- 
body interactions - electron-electron and electron-lattice, 
as well as of structural properties. Fortunately, the 
problem can be simplified by using the fact that su- 
perconductivity is a low-energy phenomenon character- 
ized by very small energy parameters {Tc/Ep, A/Ep, 
Wph/Ep) ^1. It turns out, that one can integrate 
high-energy electronic processes (which are not changed 
by the appearance of superconductivity) and then solve 
the low-energy problem by the (so called) strong-coupling 
Migdal-Eliashberg theory. It turns out that in such an 
approach the physics is separated into: (1) solving an 
ideal hand- structure Hamiltonian with the nonlocal exact 
crystal potential (sometimes called excitation potential) 
^/Bs(r, r') {IBS - ideal band structure) which includes 
the static self-energy (E^Q^(r,r',w — 0)) due to high- 
energy electronic processes, i.e. VfB5(r, r') = [T4_i(r) -|- 

VH{r)]5{T - r') + E^;;)(r,r',w = 0), with Vh the 

electron-ion and Hartree potential, respectively; (2) solv- 
ing the low-energy Eliashberg equations. However, the 
calculation of the (excited) potential V7B5'(r,r') and the 



real EPI coupling gep(r, r') = 5V/Bs(r, r')/^^-": which 
include high-energy many-body electronic processes - for 
instance large Hubbard U effects, is extremely difficult at 
present, especially in strongly correlated systems such as 
HTSC cuprates. Due to this difficulty the calculations of 
the EPI coupling in the past was usually based on the 
LDA-DFT method which will be discussed in Section 
VI in the contest of HTSC cuprates, where the nonlo- 
cal potential is replaced by the local potential Vlda(j^) - 
the ground state potential, and the real EPI coupling by 
the "local" LDA one gep{^) = ^VLDA(r)/<5R„. Since the 
exchange-correlation effects enter VLOAiic) = V'e-i(r) + 
yf/(r) -I- Vxc(r) via the local exchange-correlation po- 
tential Vxc(i") it is clear that the LDA-DFT method de- 
scribes strong correlations scarcely and it is inadequate 
in HTSC cuprates (and other strongly correlated systems 
such as heavy fermions) wher one needs an approach be- 
yond the LDA-DFT method. In Section VII we discuss 
a minimal theoretical model for HTSC cuprates which 
takes into account a minimal number of electronic or- 
bitals and strong correlations in a controllable manner 
Q . This theory treats the interplay of EPI and strong 
correlations in systems with finite doping in a system- 
atic and controllable way. The minimal model can be 
further reduced in some parameter range to the single- 
band t — J model, which allows the approximative cal- 
culation of the excited potential Vibs{^i r') and the non- 
local EPI coupling gep(r,r'). As a result one obtains the 
momentum dependent EPI coupling (7ep(k_F,q) which is 
for small hole-doping [6 < 0.3) strongly peaked at small 
transfer momenta (more precisely at q = 0) - the for- 
ward scattering peak. In the framework of this minimal 
model it is possible to explain some important proper- 
ties and resolve some puzzling experimental results, for 
instance: {a) Why is d-wave pairing realized in the pres- 
ence of strong EPI? (&) Why is the transport coupling 
constant (Xtr) rather smaller than the pairing one A, i.e. 
Xtr ^ A/3? (c) Why is the mean- field (one-body) LDA- 
DFT approach unable to give reliable values for the EPI 
coupling constant in cuprates and how many-body effects 
help; (d) Why is d-wave pairing robust in the presence 
of non-magnetic impurities and defects? (e) Why are the 
ARPES nodal and antinodal kinks differently renormal- 
ized in the superconducting states, etc? In spite of the 
encouraging successes of this minimal model, at least in 
a qualitative explanation of numerous important proper- 
ties of HTSC cuprates, we are still at present stage far 
from a fully microscopic theory of HTSC cuprates which 
is able to explain high Tc. In that respect at the and 
of Section VII we discuss possible improvements of the 
present minimal model in order to obtain at least a semi- 
quantitative theory for HTSC cuprates. 

Finally, we would like to point out that in real HTSC 
materials there are numerous experimental evidence for 
nanoscale inhomogeneities. For instance recent STM ex- 
periments show rather large gap dispersion at least on 
the surface of BISCO crystals |39i] giving rise for a pro- 
nounced inhomogeneity of the superconducting order pa- 
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rameter, i.e. A(k, R) where k is the relative momen- 
tum of the Cooper pair and R is the center of mass of 
Cooper pairs. One possible reason for the inhomogene- 
ity of A(k, R) and disorder on the atomic scale can be 
due to extremely high doping level of ~ (10 — 20) % in 
HTSC cuprates which is many orders of magnitude larger 
than in standard semiconductors (10^^ vs 10^^ carrier 
concentration). There are some claims that high Tc is 
exclusively due to these inhomogeneities (of an extrinsic 
or intrinsic origin) which may effectively increase pair- 
ing potential :40J , while some others try to explain high 
Tc solely within the inhomogeneous Hubbard or i — J 
model. Here we shall not discuss this interesting prob- 
lem but mention only that the concept of Tc increase 
by inhomogeneity is ill-defined, since the increase of Tc 
is defined with respect to the average value Tc. How- 
ever, Tc is experimentally not well defined quantity and 
the hypothesis of an increase of Tc by material inhomo- 
geneities cannot be tested at all. In studying and analyz- 
ing HTSC cuprates near optimal doping we assume that 
basic effects are realized in nearly homogeneous systems 
and inhomogeneities are of secondary role, which deserve 
to be studied and discussed separately. 

II. EPI VS NON-PHONONIC MECHANISMS 

Concerning the high Tc in cuprates, two dilemmas have 
been dominating after its discovery: [i) which interaction 
is responsible for strong quasi-particle scattering in the 
normal state - this question is related also to the dilemma 
Fermi vs no n- Fermi liquid; {ii) what is the mediating 
(gluing) boson responsible for the superconducting pair- 
ing, i.e. there is a dilemma phonons or non-phononsl In 
the last twenty-three years, the scientific community was 
overwhelmed by numerous proposed pairing mechanisms, 
most of which are hardly verifiable in HTSC cuprates. 
1. Fermi vs non-Fermi liquid in cuprates 
After discovery of HTSC in cuprates there was a large 
amount of evidence on strong scattering of quasi-particles 
which contradicts the canonical (popular but narrow) 
definition of the Fermi liquid, thus giving rise to numer- 
ous proposals of the so called non-Fermi liquids, such 
as Luttinger liquid, RVB theory, marginal Fermi liquid, 
etc. In our opinion there is no need for these radical ap- 
proaches in explaining basic physics in cuprates at least 
in optimally, slightly underdoped and overdoped metallic 
and superconducting HTSC cuprates. Here we give some 
clarifications related to the dilemma of Fermi vs non- 
Fermi liquid. The definition of the canonical Fermi liquid 
(based on the Landau work) in interacting Fermi systems 
comprises the following properties: (1) there are quasi- 
particles with charge q = ±e, spin s — 1/2 and low-laying 
energy excitations ^k(= fk — m) which are much larger 
than their inverse life-times, i.e. ^j^. >• l/rk ~ Ck/^b- 
Since the level width F = 2/rk of the quasi-particle is 
negligibly small, this means that the excited states of 
the Fermi liquid are placed in one-to-one correspondence 



with the excited states of the free Fermi gas; (2) at T = 
K there is an energy level with the Fermi surface at 
which 'CkF = ^i^d the Fermi quasi-particle distribu- 
tion function np{^j^) has finite jump at kp; (3) the num- 
ber of quasi-particles under the Fermi surface is equal to 
the total number of conduction particles (we omit here 
other valence and core electrons) - the Luttinger theo- 
rem; (4) the interaction between quasi-particles are char- 
acterized with a few (Landau) parameters which describe 
low-temperature thermodynamics and transport proper- 
ties. Having this definition in mind one can say that 
if fermionic quasi-particles interact with some bosonic 
excitation, for instance with phonons, and if the cou- 
pling is sufficiently strong, then the former are not de- 
scribed by the canonical Fermi liquid since at energies 
and temperatures of the order of the characteristic (De- 
bye) temperature kBQD{= ficojj) (for the Debye spec- 
trum ~ 0£)/5), i.e. for ^ Qd one has t^^ > and 
the quasi-particle picture (in the sense of the Landau 
definition) is broken down. In that respect an electron- 
boson system can be classified as a non- canonical Fermi 
liquid for sufficiently strong electron-boson coupling. It 
is nowadays well known that for instance Al, Zn are weak 
coupling systems since for ^j^. ~ Qd one has rj^^ <C ^k 
and they are well described by the Landau theory. How- 
ever, in (the non-canonical) cases, where for higher en- 
ergies fk ~ ©-D one has t^^ > ^j^, the electron-phonon 
system is satisfactory described by the Migdal-Eliashberg 
theory and the Boltzmann theory, where thermodynamic 
and transport properties depend on the spectral function 
a^Fs/(k, k',a;) and its higher momenta. Since in HTSC 
cuprates the electron-boson (phonon) coupling is strong 
and Tc is large, i.e. of the order of characteristic bo- 
son energies (w_b), Tc ^ ujb/^, then it is natural that 
in the normal state (at T > Tc) we deal with a strong 
interacting non-canonical Fermi liquid which is for mod- 
est non-adiabaticity parameter A < 1 described by the 
Migdal-Eliashberg theory, at least qualitatively and semi- 
quantitatively. In order to justify this statement we shall 
in the following elucidate some properties in more de- 
tails by studying optical, ARPES, tunnelling and other 
experiments in HTSC oxides. 

2. Is there limitation of the strength of EPI? 

In spite of reach experimental evidence in favor of 
strong EPI in HTSC oxides there was a disproportion in 
the research activity (especially theoretical) in the past, 
since the investigation of the SFI mechanism of pairing 
prevailed in the literature. This trend was partly due 
to an incorrect statement in [4l| on the possible upper 
limit of Tc in the phonon mechanism of pairing. Since 
in the past we have discussed this problem thoroughly in 
numerous papers - for the recent one see 42?], we shall 
outline here the main issue and results only. 

It is well known that in an electron-ion crystal, besides 
the attractive EPI, there is also repulsive Coulomb inter- 
action. In case of an isotropic and homogeneous system 
with weak quasi-particle interaction, the effective poten- 
tial Vi3//(k, w) in the leading approximation looks like as 
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for two external charges (e) embedded in the medium 
with the total longitudinal dielectric function e^(lcw) 
(k is the momentum and uj is the frequency) [4J|, |44| . 
i.e. 



47re^ 



(1) 



In case of strong interaction between quasi-particles, 
the state of embedded quasi-particles changes signifi- 
cantly due to interaction with other quasi-particles, giv- 
ing rise to ye//(k,w) ^ 47re^//c^etot(k, w). In that case 
Veff depends on other (than etoti)<.,uj)) response func- 
tions. However, in the case when Eq.([T|) holds, i. e. 
when the weak-coupling limit is realized, Tc is given by 
T, « cSexp(-l/(Aep - fi*) Hi, ig, m. Here, X,p is 
the EPI coupling constant, uj is an average phonon fre- 
quency and fi* is the Coulomb pseudo-potential, /i* = 
/i/(l + iilnEp/ui) (Ep is the Fermi energy). The cou- 
plings Xep and /i are expressed by etot (k, lu = 0) 

t^-X,p^{N{0)Veff{k,LO = 0)) 



N{0) 



2kf 



kdk 



2kl fc2etot(k,L^ = 0)' 



(2) 



where iV(0) is the density of states at the Fermi surface 
and kp IS the Fermi momentum - see more in Q. In f4l| 
it was claimed that lattice stability of the system with 
respect to the charge density wave formation implies the 
condition etot (k, w = 0) > 1 for all k. If this were correct 
then from Eq.(I2]) it follows that /i > Aep, which limits the 
maximal value of T^ to the value T™^'^ w Bp exp(— 4 — 
3/Aep). In typical metals Ep < (1 — 10) eV and if one 
accepts the statement in [41'] that Agp < /i(< 0.5), one 
obtains Tc ^ (1 — 10) K. The latter result, if it would 
be correct, means that EPI is ineffective in producing not 
only high-Tc superconductivity but also low-temperature 
superconductivity (LTS with Tc < 20 K). However, this 
result is in conflict first of all with experimental results 
in LTSC, where in numerous systems one has ^ < Xep 



and Aep > 1. 



For instance, Aep ~ 2.6 is realized in PhBi 



alloy which is definitely much higher than /i(< 1), etc. 

Moreover, the basic theory tells us that etot(k ^ 0,w) 
is not the response function [i^, [i^ (contrary to the 
assumption in [4l|). Namely, if a small external po- 
tential (JVextlk, oj) is applied to the system (of electrons 
and ions in solids) it induces screening by charges of the 
medium and the total potential is given by ^Vtot(k, u) — 
(5Ve2:t(k, a;)/etot(k, w) which means that l/£tot(k, w) is 
the response function. The latter obeys the Kramers- 
Kronig dispersion relation which implies the following 
stability condition [il], [H] 



i.e. either 



l/etot(k,cj = 0)<l, k^O, 



etot(k^O,L^ = 0) > 1 



(3) 



(4) 



etot(ky6 0,cj = 0) <0. 



(5) 



This important theorem invalidates the restriction on the 
maximal value of Tc in the EPI mechanism given in [4l| . 
We stress that the condition etot(k ^ 0, w = 0) < is not 
in conflict with the lattice stability at all. For instance, in 
inhomogeneous systems such as crystal, the total longi- 
tudinal dielectric function is matrix in the space of recip- 
rocal lattice vectors (Q), i.e. etot(k + Q, k -f Q', uj), and 
etof(k,w) is defined by £~^^{]<i,uj) = e("^J(k -f 0, k -f 0, w). 
In dense metallic systems with one ion per cell (such 
as metallic hydrogen) and with the electronic dielectric 
function eei(k, 0) and the macroscopic total dielectric 
function etot (k, 0) is given by 



etot(k, 0) = 



£ei(k,0) 



l-l/ee/(k,0)Gep(k)- 



(6) 



At the same time the energy of the longitudinal phonon 
LOi (k) is given by 



^?(k) 



J 'p 



ee/(k,0)' 



,(k,0)Gep(k)], 



(7) 



Gep is the local 



where fi^ is the ionic plasma frequenc 
(electric) field correction - see Ref. [46|. The right con- 
dition for lattice stability requires that w^(k) > 0, which 
implies that for ee/(k, 0) > one has ee/(k, O)G'ep(k) < 1. 
The latter condition gives automatically etot(k, 0) < 0. 
Furthermore, the calculations 46] show that in the metal- 
lic hydrogen (H) crystal, £(ot(k, 0) < for all k 7^ 0. 
Note, that in metallic H the EPI coupling constant is 
very large, i.e Aep « 7 and Tc may reach very large value 
Tc ~ 600 K fi^. Moreover, the analyzes of crystals with 
more ions per unit cell [1^ gives that etot(k ^ 0, 0) < is 
more a rule than an exception - see Fig. [1] The physical 
reason for etot(k 7^ 0, 0) < are local field effects de- 
scribed by Gep (k) . Whenever the local electric field E;oc 
acting on electrons (and ions) is different from the aver- 
age electric field E, i.e. E/oc 7^ E, there are corrections 
to etot(k, 0) which may lead to etot(k, 0) < 0. 



The above analysis tells us that in real crystals 
etot(k, 0) can he negative in the large portion of the Bril- 
louin zone giving rise to Aep — ^ > in Eq.®. This 
means that analytic properties of the dielectric function 
etot do not limit Tc in the phonon mechanism of pairing. 
This result does not mean that there is no limit on Tc 
at all. We mention in advance that the local field effects 
play important role in HTSC oxides, due to their layered 
structure with very unusual ionic-metallic binding, thus 
opening a possibility for large EPL 

In conclusion, we point out that there are no serious 
theoretical and experimental arguments for ignoring EPI 
in HTSC cuprates. To this end it is necessary to answer 
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FIG. 1. Inverse total static dielectric function £^^(p) for nor- 
mal metals (K, Al, Pb) and metallic H in p = (1,0,0) direc- 
tion. G is the reciprocal lattice vector. 



several important questions which are related to exper- 
imental findings in HTSC cuprates (oxides): (1) if EPI 
is important for pairing in HTSC cuprates and if super- 
conductivity is of d— wave type, how are these two facts 
compatible? (2) Why is the transport EPI coupling con- 
stant \tr (entering resistivity) rather smaller than the 
pairing EPI coupling constant Aep(> 1) (entering Tc), 
i.e. why one has \tr[^ 0.6 - 1.4) < Aep(- 2 - 3.5)? (3) 
If EPI is ineffective for pairing in HTSC oxides, in spite 
of Aep > 1, why it is so? 

3. Is a non-phononic pairing realized in HTSC? 

Regarding EPI one can pose a question - whether it 
contributes significantly to d-wave pairing in cuprates? 
Surprisingly, despite numerous experiments in favor of 
EPI, there is a believe that EPI is irrelevant for pairing 
[3]. This belief is mainly based first, on the above dis- 
cussed incorrect lattice stability criterion related to the 
sign of etot(k, 0), which implies small EPI and second, on 
the well established experimental fact that d-wave pair- 
ing is realized in cuprates , which is believed to be in- 
compatible with EPI. Having in mind that EPI in HTSC 
at and near optimal dopimng is strong with 2 < Aep < 3.5 
(see below), we assume for the moment that the lead- 
ing pairing mechanism in cuprates, which gives d-wave 
pairing, is due to some non-phononic mechanism. For 
instance, let us assume an exitonic mechanism due to 
the high energy pairing boson {flnph ^ ^ph) and with 
the bare critical temperature Tcq and look for the effect 
of EPI on Tc- If EPI is approximately isotropic, like in 
most LTSC materials, then it would be very detrimental 
for d-wave pairing. In the case of dominating isotropic 
EPI in the normal state and the exitonic-like pairing. 



then near Tc the linearized Eliashberg equations have an 
approximative form for a weak non-phonon interaction 
(with the large characteristic frequency flnph) 



Z(a;„)A„(k) w ttT^ ^ ^ 14p/i(k, q, n, to 



Z{uJn) ~ 1 + Fep/Wn 



Am(q) 



(8) 



For pure d-wave pairing with the pairing potential 

Vnph = Vnphi6w,0q) ' Q{flnph — | Wn | )6 (il„p/i — \uJn'\) 

with Kp/.(k,q) = Vo ■ Yd{9i,)Yd{9^) and YM = 
7r~-^/^ cos20k one obtains A„(k) = • Q{flnph — 
|w„|)Fd(0k) and the equation for Tc - see Q] 



T 1 
In^ « ^lf(-) 



2ttT, 



■)• 



(9) 



Here 5* is the di-gamma function. At temperatures 
near Tc one has Fep ~ 27rAepTc and the solution of 
Eq. (|ni) is approximately Tc ~ Tco exp{— Agp} with 
Tco ~ ^nphexp{-X„ph}, Xnph = ^(0)Vb. This means 
that for Tc"^^'' - 160 K and Aep > 1 the bare Tco due 
to the non-phononic interaction must be very large, i.e. 
Tco > 500 K. 

Concerning other non-phononic mechanisms, such as 
the SFI one, the effect of EPI in the framework of Eliash- 
berg equations was studied numerically in 49]. The 
latter is based on Eqs. (|99m0ip in Appendix A. with 
the kernels in the normal and superconducting chan- 
nels Afp (*i^n) and A^p, respectively. Usually, the spin- 
fluctuation kernel Xsf,kp(ii^n) is taken in the FLEX ap- 
proximation [5^. The calculations [i^ confirm the very 
detrimental effect of the isotropic (k-indepcndent) EPI 
on d-wave pairing due to SFI. For the bare SFI critical 
temperature TcO ~ 100 K and for Aep > 1 the calcula- 
tions give very small (renormalized) critical temperature 
Tc <C 100 K. These results tell us that a more realis- 
tic pairing interaction must be operative in cuprates and 
that EPI must be strongly momentum dependent [5l| . 
Only in that case is strong EPI conform with d-wave 
pairing, either as its main cause or as a supporter of 
a non-phononic mechanism. In Part II we shall argue 
that the strongly momentum dependent EPI is important 
scattering mechanism in cuprates providing the strength 
of the pairing mechanism, while the residual Coulomb in- 
teraction (by including weaker SFI) triggers it to d-wave 
pairing. 



III. EXPERIMENTAL EVIDENCE FOR 
STRONG EPI 

In the following we discuss some important experi- 
ments which give evidence for strong electron-phonon in- 
teraction (EPI) in cuprates. However, before doing it, 
we shall discuss some indicative inelastic magnetic neu- 
tron scattering measurements (IMNS) in cuprates whose 
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results in fact seriously doubt in the effectiveness of the 
phenomenological SFI mechanism of pairing which is ad- 
vocated in [s^l- First, the experimental results re- 
lated to the pronounced imaginary part of the suscep- 
tibility ImxCkjkzjio) in the normal state at and near 
the AF wave vector k = Q = (vr, tt) were interpreted in 
a number of papers as a support for the SFI mechanism 
for pairing [7] , 52] . Second, the existence of the so called 
magnetic resonance peak of /TOx(k, fc^, w) (at some en- 
ergies uj < 2A) in the superconducting state was also 
interpreted in a number of papers either as the origin of 
superconductivity or as a mechanism strongly affecting 
superconducting gap at the ant- nodal point. 



A. Magnetic neutron scattering and the spin 
fluctuation spectral function 

a. Huge rearrangement of the SFI spectral function 
and little change of Tc 

Before discussing experimental results in cuprates on 
the imaginary part of the spin susceptibility /mx(k, u) 
we point out that in (phenomenological) theories based 
on spin fluctuation interaction (SFI) the quasi-particle 
self-energy Es^(k, cli„) (a;„ is the Matsubara frequency 
and f is the Nambu matrix) in the normal and super- 
conducting state and the effective (repulsive) pairing po- 
tential VsfjM., to) (where «w„ — > w + iyy) are assumed in 
the form \j\ 

i;s/(k,w„) = X! ^s/(k- k',a;„-Wm)ToG(k',a;m)fo, 

k' m 



dv /mx(q, v + iO^) 



(10) 



Although the form of Vgf can not be justified theoret- 
ically, except in the weak coupling limit {gsf ^ Wb) 
only, it is used in the analysis of the quasi-particle 
properties in the normal and superconducting state of 
cuprates where the spin susceptibility (spectral function) 
Imx(q, w) is strongly peaked at and near the AF wave 
vector Q — (7r/a,7r/a). Can the pairing mechanism 
in HTSC cuprates be explained by such a phenomenol- 
ogy and what is the prise for it? The best answer is 
to look at the experimental results related to inelastic 
magnetic neutron scattering (IMNS) experiments which 
give /mx(q, w). In that respect very indicative and im- 
pressive IMNS measurements on YBa2Cu30Q+x, which 
are done by Bourges group Ts'] , demonstrate clearly that 
the normal state susceptibility /mx^°'''^^(q, w) (the odd 
part of the spin susceptibility in the bilayer system) at 
q = Q = (tt, 7r) is strongly dependent on the hole doping 
as it is shown in Fig. 

The most pronounced result for our discussion is that 
by varying doping there is a huge rearrangement of 
, uj) in the normal state, especially in the en- 
ergy (frequency) region which might be important for 
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FIG. 2. Magnetic spectral function /mx'~Hk,a;) in the 
normal state of Y Ba2Cu30e+x a,t T — 100 K and at 
Q = (7r,7r). 100 counts in the vertical scale corresponds to 
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superconducting pairing, let say meV < w < 60 meV. 
This is clearly seen in the last two curves in Fig. [2] 
where this rearrangement is very pronounced, while at 
the same time there is only small variation of the crit- 
ical temperature Tc- It is seen in Fig. [2] that in the 
underdoped YBa2Cu30Q,<)2 crystal Imx(°'^''^(Q,a;) and 
S{Q) = N{^i)g^^j J^" dujImx^°'^''HQ,uj) are much larger 
than that in the near optimally doped YBa^Cu^O^.g'j, 
i.e. one has S'6.92(Q) ^ <5'6.97(Q), although the dif- 
ference in the corresponding critical temperatures Tc is 
very small, i.e. yj^'^^' = 91 if (in Y Ba2Cu^OQ,g2) and 
rp{fi.9i) _ g2 g YBa2Cu30Q,g7). This pronounced 

rearrangement and suppression 

of Imx^^'^'^HQ-'^) (in the 
normal state of YBCO) by doping toward the optimal 
doping but with negligible change in Tc is strong evidence 
that the SFI pairing mechanism is not the dominating 
one in HTSC cuprates. This insensitivity of T^, if inter- 
preted in term of the SFI coupling constant Xsf{^ dif): 



means that the latter is small, i.e. A 



(cxp) 



f) 

< 1. We 
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stress that the explanation of high Tc in cuprates by 
the SFI phenomenological theory |7j] assumes very large 
SFI coupling energy with g^^^ « 0.7 eV while the fre- 
quency (energy) dependence of /tox(Q,w) is extracted 
from the fit of the NMR relaxation rate Tj"^ which gives 



(NMR) 



100 K 0. To this point, the NMR mea- 
surements (of T^^) give that there is an anti- correlation 
between the decrease of the NMR spectral function Iq = 
limi^^o Imx^^^^^'^ (Q, a;)/w and the increase of Tc by in- 
creasing doping toward the optimal one - see [1] and 
References therein. The latter result additionally dis- 
favors the SFI model of pairing [7| since the strength of 
pairing interaction is little affected by SFI. Note, that if 
instead of taking /TOx(Q,a;) from NMR measurements 
one takes it from IMNS measurements, as it was done 
in [HI], than for the same value g'fj^^one obtains much 
smaller Tc- For instance, by taking the experimental val- 
ues for /tox^^^'^^'^HQi ^) ill underdoped YBa2Cu30e q 
with Tc ^ 60 K one obtains ri^^'^^^ < ri^^^^Vs 
p^, while ry^'^^'' ^ 50 K for gf/^ » 1. The situ- 
ation is even worse if one trays to fit the resistivity with 
,uj) in YBa2Cu306,6 since this fit gives 
rpiiMNS) ^ Y K. These results point to a deficiency of 



the SFI phenomenology (at least that based on Eq.lTC 
to describe pairing in HTSC cuprates. 

Having in mind the results in ^53^, the recent theo- 
retical interpretation in [52j of IMNS experiments 
and ARPES measurements Q on the underdoped 
Y Ba2Cu^0Q,Q in terms of the SFI phenomenology de- 
serve to be commented. The IMNS experiments [s^l 
give evidence for the "hourglass" spin excitation spec- 
trum (in the superconducting state) for the momenta 
q at, near and far from Q, which is richer than the 
common spectrum with magnetic resonance peaks mea- 
sured at Q. In [53| the self-energy of electrons due to 
their interaction with spin excitations is calculated by 
using Eq. pUj) with g^j = (3/2)[7^ and /mx(q, w) taken 
from ^5^]. However, in order to fit the ARPES self- 
energy and low-energy kinks (see discussion in Section 
C) the authors of 52] use very large value U = 1.59 
eV , i.e. much larger than the one used in [53j . Such a 
large value of U has been used earlier within the Monte 
Carlo simulation and the fit of the Hubbard model [5^ . 
In our opinion this value for U is unrealistically large 
in the case of strongly correlated systems where spin- 
fluctuations are governed by the effective electron ex- 
change interaction Jcu-Cu ^ 0.15 eV [13]. This implies 
that [/ < 1 eV" and Tc < 10 K. Note, that this value 
for Jcu-Cu comes out also from the theory of strongly 
correlated electrons in the three-band Emery model with 
Jcu-Cu ~ (4i^d/(Adp + Upd)^)[il/Ud) + 2/ (Up + 2A)] - 
for parameters see Part II Section VII. 

Concerning the problem related to the rearrangement 
of the SFI spectral function /r7ix(Q,a;) in YB2CU3OQ+X 
[l8i] we would like to stress, that in spite that the lat- 
ter results were obtained ten years ago they are not dis- 



puted by the new IMNS measurements [l^] on high qual- 
ity samples of the same compound (where much longer 
counting times were used in order to reduce statistical 
errors). In fact the results in Ts*] are confirmed in [l9j 
where the magnetic intensity /(q, w)(^ /mx(q, w)) (for 
q at and in the broad range of Q) for the optimally 
doped Y Ba2Cu30e.95 (with Tc = 93 K) is at least three 
times smaller than in the underoped Y Ba2Cu^OQ,Q with 
Tc — 60 K. This result is again very indicative sign 
of the weakness of SFI since such a huge reconstruction 
would decrease Tc in the optimally doped Y Ba2Cuy,OQ,Qz 
if analyzed in the framework of the phenomenological SFI 
theory based on Eq. . It also implies that due to the 
suppression of /mx(q, uj) by increasing doping toward 
the optimal one a straightforward extrapolation of the 
theoretical approach in [52.] to the explanation of Tc in 
the optimally doped Y Ba2Cu^OQ,Qz would require an in- 
crease of JJ to the value even larger than 4 eV , what is 
highly improbable. 

h. Ineffectiveness of the magnetic resonance peak 
A less direct argument for smallness of the SFI cou- 
pling constant, i.e. g^^^ < 0.2 eV and g^^^ ^ gsf 
comes from other experiments related to the magnetic 
resonance peak in the superconducting state, and this 
will be discussed next. In the superconducting state of 
optimally doped YBCO and BISCO, /r7ix(Q,a;) is sig- 
nificantly suppressed at low frequencies except near the 
resonance energy ujres ~ 41 meV where a pronounced 
narrow peak appears - the magnetic resonance peak. We 
stress that there is no magnetic resonance peak in some 
families of HTSC cuprates, for instance in LSCO, and 
consequently one can question the importance of the res- 
onance peak in the scattering processes. The experiments 
tell us that the relative intensity of this peak (compared 
to the total one) is small, i.e. /q ~ (1 — 5)% - see Fig [31 
In underdoped cuprates this peak is present also in the 
normal state as it is seen in FigO 

After the discovery of the resonance peak there were 
attempts to relate it first, to the origin of the super- 
conducting condensation energy and second, to the kink 
in the energy dispersion or the peak-dimp structure in 
the ARPES spectral function. In order that the con- 
densation energy is due to the magnetic resonance it is 
necessary that the peak intensity Iq is small [5^. Jq 
is obtained approximately by equating the condensa- 
tion energy Econ ~ iV(0)A^/2 with the change of the 
magnetic energy E^ag in the superconducting state, i.e. 



Efnaq — J 



dud ^ cosfc^ — COS ky)S{'k,uj). (11) 



(27r)3 



By taking A « 2Tc and the realistic value A^(0) ^ 
1/(10J) ~ 1 states/eV ■ spin, one obtains Iq ~ 
10"^(rc/J)^ ^ 10^^. However, such a small intensity 
can not be responsible for the anomalies in ARPES and 
optical spectra since it gives rise to small coupling con- 
stant Xsf^res for the interaction of holes with the reso- 
nance peak, i.e. Xsf,res ~ {2IoN {0)g1j j iOres) ^ 1- Such 
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FIG. 3. Magnetic spectral function Imx^~'' in the su- 

perconducting state of 5^-802(711306+2: at T = 5 A' and at 
Q = (tt, tt). 100 counts in the vertical scale corresponds to 
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a small coupling does not affect superconductivity at all. 
Moreover, by studying the width of the resonance peak 
one can extract order of magnitude of the SFI coupling 
constant gsf ■ Since the magnetic resonance disappears in 
the normal state of the optimally doped YBCO it can be 
qualitatively understood by assuming that its broadening 
scales with the resonance energy ujres, i-e. Jres < ^res, 
where the line-width is given by 7^^^ = 4Tr{N{0)gsf)'^uJres 
[58| . This condition limits the SFI coupling to gsf < 0.2 
eV. We stress that the obtained gsf is much smaller 
(at least by factor three) than that assumed in the phe- 
nomenological spin-fluctuation theory 17] and [52] where 
gsf ^ 0.6 — 0.7 eV and U ~ 1.6 eV, but much larger 
than estimated in [58| (where gsf < 0.02 eV). The small- 
ness of gsf comes out also from the analysis of the anti- 
ferromagnetic state in underdoped metals of LSCO and 
YBCO t59|, where the small (ordered) magnetic moment 
/i(< 0.1 Hg) points to an itinerant antiferromagnetism 
with small coupling constant gsf < 0.2 eV. The con- 



clusion from this analysis is that in the optimally doped 
YBCO the sharp magnetic resonance is a consequence of 
the onset of superconductivity and not its cause. There 
is also one principal reason against the pairing due to 
the magnetic resonance peak at least in optimally doped 
cuprates. Since the intensity of the magnetic resonance 
near Tc is vanishingly small, though not affecting pair- 
ing at the second order phase transition at Tc, then if it 
would be solely the origin for superconductivity the phase 
transition at Tc would be first order, contrary to exper- 
iments. Recent ARPES experiments give evidence that 
the magnetic resonance cannot be related to the kinks in 
ARPES spectra [g^l, (Sli - see the discussion below. 

Finally, we would like to point out that the recent 
magnetic neutron scattering measurements on optimally- 
doped large-volume crystals Bi2Sr2CaCu20s+s [62| . 
where the absolute value of /mx(q, 1^) is measured, are 
questioning the interpretation of the electronic mag- 
netism in cuprates in terms of itinerant magnetism. This 
experiment shows a lack of temperature dependence of 
the local spin susceptibility /mx(a;) = J2q -f'7^x(q, w) 
across the superconducting transition Tc — 91 K, i.e 
there is only a minimal change in /mx(w) between 10 K 
and 100 K. Note, if the magnetic excitations were due 
to itinerant quasi-particles we should have seen dramatic 
changes of /tox(w) as a function of T over the whole 
energy range. This T-independence of Im~x{ijj) strongly 
opposes many theoretical results in [l^, [ij], [l^ which 
assume that the bosonic spectral function is proportional 
to /mx(a;) that can be extracted from optic measure- 
ments. This procedure gives that Imx{ijj) is strongly 
T-dependent contrary to the experimental results in |62| 
- see more in Subsection B on optical conductivity. 



B. Optical conductivity and EPI 

Optical spectroscopy gives information on optical con- 
ductivity cr(w) and on two-particle excitations, from 
which one can indirectly extract the transport spectral 
function a1^.F[uj). Since this method probes bulk sam- 
ple (on the skin depth), contrary to ARPES and tun- 
nelling methods which probe tiny regions (10 - 15 A) 
near the sample surface, this method is indispensable. 
However, one should be careful not over-interpreting the 
experimental results since cr(w) is not a directly measured 
quantity but it is derived from the reflectivity i?(w) — 

(yi~H-i)/(VM^ + i) 

electric tensor £m(w) 
Here, 



with the transversal di- 



£m,oo + £ii,iatt + ^TTicrii{u:)/uj. 
£m,oo is the high frequency dielectric function, 
£ii,iatt describes the contribution of the lattice vibrations 
and aii{uj) describes the optical (dynamical) conductivity 
of conduction carriers. Since R{uj) is usually measured 
in the limited frequency interval Wmin < a; < Wmax some 
physical modelling for R{uj) is needed in order to guess it 
outside this range - see more in reviews 0,0- This was 
the reason for numerous misinterpretations of optic mea- 
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surements in cuprates, that will be uncover below. An 
illustrative example for this claim is large dispersion in 
the reported value oicupi - from 0.06 to 25 eV, i.e. almost 
three orders of magnitude. However, it turns out that IR 
measurements of R{uj) in conjunction with elipsometric 
measurements of eii(w) at high frequencies allows more 
reliable determination of (j{uj) f63|. 

1. Transport and quasi-particle relaxation rates 
The widespread misconception in studying the quasi- 
particle scattering in cuprates was an ad hoc assumption 
that the transport relaxation rate "/^^.(lu) is equal to the 
quasi-particle relaxation rate 7(w), in spite of the well 
known fact that the inequality itAoj) ^ 7(w) holds in 
a broad frequency (energy) region [20|. This (incorrect) 
assumption was one of the main arguments against EPI 
as relevant scattering mechanism in cuprates. Although 
we have discussed this problem several times before, we 
do it again due to the importance of this subject. 

The dynamical conductivity cr(cLi) consists of two parts, 
i.e. cr(w) = cr™*'^'~(w) + CT™*™(cj) where cr'"*'=''(cj) de- 
scribes inter-hand transitions which contribute at higher 
than intra-band energies, while a™^^°-(Lij) is due to intra- 
band transitions which are relevant at low energies uj < 
(1 — 2) eV. (Note, that in the IR measurements the fre- 
quency is usually given in cm~^, where the following con- 
version holds Icm-i = 29.98 GHz = 0.123985 meV = 
1.44 K.) The experimental data for a{Lu) — ai -\- ia2 
in cuprates are usually processed by the generalized (ex- 
tended) Drude formula [l^l, (HI, [22], [gl], 



V UJ -\- V , , U} — V , , ^ 

-2 nB{'jJ + v)^ nsiuJ-v)]- (14) 

(jj U) id 

Here nsiuj) is the Bose distribution function. For com- 
pleteness we give also the explicit form of the transport 
mass mtriuj) see 0, g, [H, [H, [H, . 



oo 



U> 



with the Kernel K{x,y) = [i/y) -\- {{{y — x)/x)['ili{l — 
ix -\- iy) — -\- iy)]} — {y ^ —y} where ip is the di- 
gamma function. In the presence of impurity scattering 
one should add Jimp,tr to ^tr- It turns out that Eq.([T4]) 
holds within a few percents also for large Aep(> 1). Note, 
that ttj,. iFi{v) ^ afFi{v) and the index I enumerates all 
scattering bosons - phonons, spin fluctuations, etc. For 
comparison, the quasi-particle scattering rate 7(0;, T) is 
given by 



7(0;, T) = 27r / dva^F{v){2nB{v) 



-I- npi^ + oj) + np{i' — uj)} + 7'' 



(15) 



where up is the Fermi distribution function. For com- 
pleteness we give also the expression for the quasi-particle 
effective mass to*(w) 



^ I (12) 

Att - tUJmtr{uj)/moo Wtr(w)' 

which is an useful representation for systems with sin- 
gle band electron-boson scattering which is justified in 
HTSC cuprates. However, this procedure is inadequate 
for interpreting optical data in multi-band systems such 
as new high-temperature superconductors Fe-based pnic- 
tides since even in absence of the inelastic intra- and 
inter-band scattering the effective optic relaxation rate 
may be strongly frequency dependent [65j . (The useful- 
ness of introducing the optic relaxation ujtriuj) will be 
discussed below and in Appendix B.) Here, i — a,b enu- 
merates the plane axis, Wp, jf^{u!,T) and mop{uj) are 
the electronic plasma frequency, the transport (optical) 
scattering rate and the optical mass, respectively. Very 
frequently it is analyzed the quantity 7j^(a;,r) given by 

M 



ltA^,T) = — 7i^(cj,T) - 



mtr{uj) 



uiRea{uj) ' 



which is determined from the half-width of the Drude- 
like expression for cr(aj) and is independent of w^. In the 
weak coupling limit Agj, < 1, the formula for conductivity 
given in Appendix B Eqs. (|125m28l) can be written in 
the form of Ea. ((T2|) where 7(,. reads [21]- 22] 



, Jo 



d,.al,Fi{^)[2{l 



00 

1 + - dQafFiin) 



„ r /I UJ + Q, ,,1 U}—fl,^ 

The term ^""p is due to the impurity scattering. By 
comparing Ea.([T^ and Ea.([T^. it is seen that 7^^ and 
7 are different quantities, i.e. 7^^ ^ 7, where the for- 
mer describes the relaxation of Bose particles (excited 
electron-hole pairs) while the latter one describes the re- 
laxation of Fermi particles. This difference persists also 
at T = K where one has (due to simplicity we omit in 
the following summation over / ) pcj] 



27r 

UJ 



dy{uj - v)al{y)F{u) (17) 



(13) and 



7(0;) = 27r / dva''{v)F{v). 
Jo 



(18) 



In the case of EPI with the constant electronic density 
of states, the above equations give that j^p{uj) — const 
for 



cj > uj'^i^^ while Jep.tri^) (^s well as "/Ip^tr) is mono- 



tonic growing for w > 



.max 

-'ph ' WllClC UJpf^ 

mal phonon frequency So, the growing of 7^^ tri^) (^^nd 
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FIG. 4. (top) Scattering rates 7(aj,r), 7j,,(cj,r) and 7*^ 
- from top to bottom, for the Eliashberg function in (bot- 
tom). From [SJl. (bottom) Eliashberg spectral function 
a'lp{uj)Fph(uj) obtained from tunnelling experiments on break 
junctions [27|. Inset shows 7^^ with (full line) and without 
(dashed line) inter-band transitions 0]. 

7*p t^) for uj > uj'pf^^ is ubiquitous and natural for EPI 
scattering and has nothing to do with some exotic scat- 
tering mechanism. This behavior is clearly seen by com- 
paring 7(0;, T), 7tr(w, T) and 7^^ which are calculated for 
the EPI spectral function alp{uj)Fph{Lo) extracted from 
tunnelling experiments in YBCO (with w^f^ ~ 80 meV) 
[13 - see Fig. m 

The results shown in Fig. |4] clearly demonstrate the 
physical difference between two scattering rates 7gp and 
7ep,*r (and 74^)- It is also seen that 7tr(w, T) is even more 



linear function of oj than ^^^{uj,T). From these calcula- 
tions one concludes that the quasi-linearity of j^^(uj,T) 
(and 7j^) is not in contradiction with the EPI scatter- 
ing mechanism but it is in fact a natural consequence 
of EPI. We stress that such behavior of 7gp and 7ep,tr 
(and 7*p j^), shown in Fig. |4l is in fact not exceptional 
for HTSC cuprates but it is generic for many metallic 
systems, for instance 3D metallic oxides, low tempera- 
ture superconductors such as Al, Pb, etc. - see more in 
[1], [1] and References therein. 

Let us discuss briefly the experimental results for R{uj) 
and 7(^(a;,T) and compare these with theoretical pre- 
dictions obtained by using a single band model with 
a1p{u})Fph{(jj) extracted from the tunnelling data with 
the EPI coupling constant Agp > 2 [23|. In the case of 
YBCO the agreement between measured and calculated 
R{lS) is very good up to energies lo < 6000 cm~^ which 
confirms the importance of EPI in scattering processes. 
For higher energies, where a mead-infrared peak appears, 
it is necessary to account for inter-band transitions [31 . 
In optimally doped Bi2Sr2CaCu20G {Bi - 2212) ^ 
the experimental results for 74^(0;, T) are explained the- 
oretically by assuming that the EPI spectral function 
a1p{u!)F{u!) ~ Fph{ui), where Fph{ui) is the phononic den- 
sity of states in BISCO while a1p{w) ^ oj^'^, Xep — 1.9 
and 7ij„p « 320 cm~^ - see Fig. [Sl^top). The agreement 
is rather good. At the same time the fit of ^Ij.{uj,T) by 
the marginal Fermi liquid phenomenology fails as it is 
evident in Fig. [5Ubottom). 

Now we will comment the so called pronounced lin- 
ear behavior of 7(^(a;,r) (and 7j^(a;,T)) which was one 
of the main arguments for numerous inadequate conclu- 
sions regarding the scattering and pairing bosons and 
EPI. We stress again that the measured quantity is re- 
flectivity i?(a;) and derived ones are cr(a;), 74^(0;, T) and 
TOtr(w), which are very sensitive to the value of the dielec- 
tric constant Ecx) ■ This sensitivity is clearly demonstrated 
in Fig. [S] for Bi-2212 where it is seen that 7(^(a;,T) (and 
7t^(w, T)) for = 1 is linear up to much higher w than 
in the case £00 > 1- 

However, in some experiments [63|, [1^ the extracted 
7t^(w,r) (and 74^(0;, T)) is linear up to very high uj sa 
1500 cm~^ . This means that the ion background and 
inter-band transitions (contained in £00) are not properly 
taken into account since too small £00 (?J 1) is assumed. 
The recent elipsometric measurements on YBCO [69^ give 
the value £00 « 4 — 6, which gives much less spectacular 
linearity in the relaxation rates 7tj.(w, T) (and 7tj.(a;, T)) 
than it was the case immediately after the discovery of 
HTSC cuprates, where much smaller £00 was assumed. 

Furthermore, we would like to comment on two points 
concerning cr, 7^^, 7 and their interrelations. First, the 
parametrization of cr{uj) with the generalized Drude for- 
mula in Eq. and its relation to the transport scatter- 
ing rate 7(^(a;,T) and the transport mass mt,.(a;,T) is 
useful if we deal with electron-boson scattering in a sin- 
gle band problem. In [25|, [65*] it is shown that cr(w) of 
a two-band model with only elastic impurity scattering 
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FIG. 5. (top) Experimental transport scattering rate 7*^ 
(solid lines) for BISCO and the theoretical curve by using 
Eq. p25|l and transport mass m*r with a'^F{uj) due to EPI 
which is described in text (dashed lines), (bottom) Compar- 
ison with the marginal Fermi liquid theory - dashed lines. 
From [|. 



can be represented by the generalized (extended) Drude 
formula with uj and T dependence of effective parame- 
ters {uj.T), rall^{uj,T) despite the fact that the in- 
elastic electron-boson scattering is absent. To this end 
we stress that the single-band approach is justified for a 




1500 



FIG. 6. Dependence of 7j^(a;,r) on e^o in Bi2Sr2CaCu20s 
for different temperatures: Soo ~ 4 (solid lines) and £00 = 1 
(dashed lines). On the horizontal is lo in units cm~^ . From 



number of HTSC cuprates such as LSCO, BISCO etc. 
Second, at the beginning we said that 7j^(a;,T) and 
7(cj, T) are physically different quantities and it holds 
7j^(aj, T) 7^ j{uj,T). In order to give the physical picture 
and qualitative explanation for this difference we assume 
that al^F{v) k. a'^F{v). In that case the renormalized 
quasi-particle frequency uj{uj) = Z{ll!)uj — uj — S(cli) and 
the transport one (I'tr(w) - defined in Eg (fT2|) . are related 
and at T = they are given by [1^ , (2^ 



1 



(19) 



(For the definition of Z{uj) see Appendix A.) It gives the 
relation between 7f^(w) and 7(1^), mtr{uj) and m*(w), 
respectively 



1 



rfw'7(w') 



dw'2w'm*(w')- 



(20) 



(21) 



The physical meaning of Eg. dTSl) is the following: in op- 
tical measurements one photon with the energy uj is ab- 
sorbed and two excited particles (electron and hole) are 
created above and below the Fermi surface. If the elec- 
tron has energy uj' and the hole uj — uj\ then they relax as 
quasi-particles with the renormalized frequency uj. Since 
ui' takes values < w' < a; then the optical relaxation 
uJtr{i^) is the energy-averaged u}{uj) according to Ea.([TO|. 
The factor 2 is due to the two quasi-particles - electron 
and hole. At finite T, the generalization reads [23|j [1^ 

ujtri^) — — / duj'[l — npi^') — npi^iJ — uj')]2uj{uj'). 
w Jo 

(22) 

2. Inversion of the optical data and a'^^{uj)F{uj) 
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In principle, the transport spectral function 
afj.{tjj)F{uj) can be extracted from a{uj) (or 7j^(cij)) 
only a.tT ^0 K, which follows from Eq.([T7l) 



87r2 duj- 



\ujRe 



a{uj) 



(23) 



or equivalently a'^^{ijj)F{ijj) = {l/2'K)d'^{ijj^^^{uj))/duj'^. 
However, real measurements are performed at finite T 
(at T > Tc which is rather high in HTSC cuprates) 
and the inversion procedure is an ill-posed problem since 
a'^^{uj)F{uj) is the de-convolution of the inhomogeneous 
Fredholm integral equation of the first kind with the tem- 
perature dependent Kernel K2{ij-',i',T) - see Ea. ([Ti| . It 
is known that an ill-posed mathematical problem is very 
sensitive to input since experimental data contain less in- 
formation than one needs. This procedure can cause first, 
that the fine structure of a'^j.{u!)F{u!) gets blurred (most 
peaks are washed out) in the extraction procedures and 
second, the extracted a^^(uj)F{uj) is temperature depen- 
dent even when the true afj.{uj)F{uj) is T independent. 
This artificial T-dependence is especially pronounced in 
HTSC cuprates because Tc{^ 100 K) is very high. In the 
context of HTSC cuprates, this problem was first stud- 
ied in [2lj, [i^l where this picture is confirmed by the 
following results: (1) the extracted shape of aj^(aj)F(aj) 
in Y Ba2Cu307-x (and other cuprates) is not unique and 
it is temperature dependent, i.e. at higher T > Tc the 
peak structure is smeared and only a single peak (slightly 
shifted to higher w) is present. For instance, the exper- 
imental data of i?(w) in YBCO were reproduced by two 
different spectral functions a'fj.{uj)F{uj), one with single 
peak and the other one with three peaks structure as it is 
shown in Fig. [3 where all spectral functions give almost 
identical R{uj). The similar situation is realized in opti- 
mally doped BISCO as it is seen in Fig. |5] where again 
different functions a'^{LL>)F{LL>) reproduce very well curves 
for R{lu) and cr('j-'). However, it is important to stress 
that the obtained width of the extracted a^^{Lu)F{Lu) in 
both compounds coincide with the width of the phonon 
density of states Fph{uj) HH, [H], [zT 
energy bound for 



(2) The upper 
^j.{u})F{u}) is extracted in [ 



?)Th 

m 



[2^ and 

it coincides approximately with the maximal phonon fre- 
quency in cuprates ^ 80 meV as it is seen in Figs. 

These results demonstrate the importance of EPI in 
cuprates [2l|, [l^. We point out that the width of 
af^{uj)F(uj) which is extracted from the optical measure- 



ments [21|, [22| coincides with the width of the quasi- 
particle spectral function a'^{uj)F{uj) obtained in tun- 
nelling and ARPES spectra (which we shall discuss be- 
low), i.e. both functions are spread over the energy in- 
terval < oj < uj'^h^i^ 80 meV). Since in cuprates this 
interval coincides with the width in the phononic den- 
sity of states F{ljj) and since the maxima of q^{uj)F{lxj) 
and F{ijj) almost coincide, this is further evidence for the 
importance of EPI. 

To this end, we would like to comment two aspects 
which appear from time to time in the literature. First, in 
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FIG. 7. Experimental (solid lines) and calculated (dashed 
lines) of R{ijj) in optimally doped YBCO XI] at T=100, 200, 
300 K (from top to bottom). Inset: the two (solid and dashed 
lines) reconstructed q^^{uj)F{uj) at T=100 K. The phonon 
density of states F{uS) - dotted line in the inset. From !21]. 




FIG. 8. Experimental (solid line) and calculated (dashed line) 
of Riu)) in optimally doped BISCO [zl at T=100 K. Inset: 
the reconstructed a^r{uji)F{LLi) - solid line. The phonon density 
of states F{ijj) - dotted line. From [2j|. 



some reports [1J|, [Tj|, [T^ it is assumed that a1^{uj)F{uj) 
of cuprates can be extracted also in the superconducting 
state by using Eq. (^51 . However, Eq. (^51) holds exclu- 
sively in the normal state (at T = 0) since criuj) can be 
described by the generalized (extended) Drude formula in 
Eq. (fT2|) only in the normal state. Such an approach does 
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not hold in the superconducting state since the dynami- 
cal conductivity depends not only on the electron-boson 
scattering but also on coherence factors and on the mo- 
mentum and energy dependent order parameter A(k, w). 
In such a case it is unjustified to extract at^{u!)F{u!) from 
Eq. (|23p. Second, if (and cr[uj)) in cuprates are due 
to some other bosonic scattering which is pronounced up 
to much higher energies lUc ^ '^p]^^, this should be seen 
in the extracted spectral function a'j^{uj)F{uj). Such an 
assumption is made, for instance, in the phenomenologi- 
cal spin-fluctuation approach Q where it is assumed that 
a'^{u})F{uj) g^j:lmx{uj) where Imx(a;) is extended up to 
the large energy cutoff ujc ~ 400 meV. This assumption 
is in conflict with the above theoretical and experimental 
analysis which shows that solely EPI can describe -R(w) 
very well and that the contribution from higher energies 
w ^ '^^ZlT ™'^st be small and therefore irrelevant for 
pairing ]2l| , [g^] , [z3| • The finding of the importance of 
EPI is also confirmed by tunnelling measurements - see 
discussion in Subsection D. Despite the experimentally 
established facts, that the energy width of the extracted 
a'jj,{iu)F{u!) coincides with the phononic range, which fa- 
vor EPI, some reports appeared recently claiming that 
SFI dominates and that af^{uj)F{uj) ~ gl^Imxi'-o) where 
Imxil^) = J c?^fcx(k, w) [l3|, [l^. This claim is based on 
reanalyzing of some IR measurements [3] , [l^ and the 
transport spectral function Q;j^(a;)i^(a;) is extracted in 
by using the maximum entropy method in solving 
the Fredholm equation. However, in order to exclude 
negative values in the extracted Q;(^(a;)F(a;), which is an 
artefact and due to the chosen method, in [l3| it is as- 
sumed that af^{uj)F{uj) has a rather large tail at large 
energies - up to 400 meV. It turns out that even such an 
assumption in extracting a'}^(uj)F(uj) does not reproduce 
the experimental curve Imx{uj) [73j in some important 
respects: (i) the relative heights of the two peaks in the 
extracted spectral function cej^{Lo)F{u!) at lower temper- 
atures are opposite to the experimental curve Imx{(^) 
ff^ - see Fig. 1 in [l3|- {ii) the strong temperature de- 
pendence of the extracted atj.{ui)F{ui) found in [l3|, [l^ 
is in fact not an intrinsic property of the spectral func- 
tion but it is an artefact due to the high sensitivity of the 
extraction procedure on temperature. As it is already ex- 
plained before, this is due to the ill-posed problem of solv- 
ing the Fredholm integral equation of the first kind with 
strong T-dependent kernel. Third, in fact the extracted 
spectral weight Ib{i^) = cttri^)P{^) H much 
smaller values at larger frequencies (w > 100 meV) than 
it is the case for the measured Imx(w), i.e. (Ib{<^ > 100 
TOey)//_B(a;,nax)}< Imxi^ > 100 mey)//mx(wmax) - 
see Fig. 1 in [1J|. Fourth, the recent magnetic neu- 
tron scattering measurements on optimally-doped large- 
volume crystals Bi2Sr2CaCu20s+s [13 (where the ab- 
solute value of /TOx(q, w) is measured) are not only 
questioning the theoretical interpretation of magnetism 
in HTSC cuprates in terms of itinerant magnetism but 
they also oppose the finding in [l3|-[l3- This experi- 
ment shows that the local spin susceptibility Imxi^) = 



J2q I^xil' ^) is temperature independent across the su- 
perconducting transition Tc — 91 K, i.e there is only 
a minimal change in Imx{ui) between 10 K and 100 
K. This T-independence of Irn y juj) st rong ly opposes 



the (above discussed) results in [isj, [ij], 15 1, where the 
fit of optic measurements give strong T-dependence of 
Imxiuj). 

Fifth, the transport coupling constant Xtr extracted in 
is too large, i.e. Xtr > 3 contrary to the previous 
findings that A*^ < 1.5 - 2 [2l|, [H, [t^. Since in HTSC 
one has A > Xtr this would probably give A « 6 what 
is not confirmed by other experiments. Sixth, the inter- 
pretation of at^{u})F{u}) in LSCO and BISCO solely in 
terms of Imxi^j) is in contradiction with the magnetic 
neutron scattering in the optimally doped and slightly 
underdoped YBCO [18]. The latter was discussed in Sub- 
section A, where it is shown that Imx(Q,a;) is small in 
the normal state and its magnitude is even below the ex- 
perimental noise. This means that if the assumption that 
atr{i^)F{u!) « g^j:lrax{uj) were correct then the contribu- 
tion to Imx(a;) from the momenta < k <^ Q would be 
dominant which is very detrimental for d-wave supercon- 
ductivity. 

Finally, we point out that very similar (to cuprates) 
properties, of <t{uj), R{u}) (and p{T) and electronic Ra- 
man spectra) were observed in 3D isotropic metallic ox- 
ides Lao,5Sro,5Co03 and CaQ ^SrQ c^RuOs which are non- 
superconducting [TJI and in Bai^xK^BiO^ which super- 
conducts below Tc ~ 30 if at a; = 0.4. This means that in 
all of these materials the scattering mechanism might be 
of similar origin. Since in these compounds there are no 
signs of antiferromagnetic fluctuations (which are present 
in cuprates), then EPI plays important role also in other 
oxides. 

3. Restricted optical sum-rule 

The restricted optical sum-rule was studied intensively 
in HTSC cuprates. It shows some peculiarities not 
present in low-temperature superconductors. It turns out 
that the restricted spectral weight W{flc,T) is strongly 
temperature dependent in the normal and superconduct- 
ing state, that was interpreted either to be due to EPI 
24 1, [2^ or to some non-phononic mechanisms [tH]. In 
the following we demonstrate that the temperature de- 
pendence of Wiflc, T) = ly (0) - fiT'^ in the normal state 
can be explained in a natural way by the T-dependence 
of the EPI transport relaxation rate 74^(0;, T) [H], [H]. 
Since the problem of the restricted sum- rule has attracted 
much interest it will be considered here in some details. 
In fact, there are two kinds of sum rules related to cr(a;). 
The flrst one is the total sum rule which in the normal 
state reads 



'pi 



nne 
2m 



(24) 



while in the superconducting state it is given by the 
Tinkham-Ferrell-Glover (TFG) sum-rule 
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CTl {Uj)duj = — 2- 

oAj^ 



fff = (25) 



+0 



Here, n is the total electron density, e is the electron 
charge, m is the bare electron mass and Xl is the London 
penetration depth. The first (singular) term c^/8A^ in 
Eq. psp is due to the superconducting condensate which 
contributes <jf condi^) — (c^/4^i)<5(w). The total sum 
rule represents the fundamental property of matter - the 
conservation of the electron number. In order to calculate 
it one should use the total Hamiltonian -fftot = + ^mt 
where all electrons, electronic bands and their interac- 
tions Hint (Coulomb, EPI, with impurities, etc.) are ac- 
counted for. Here, is the kinetic energy of bare elec- 
trons 



vertex correction H(ia;„) (at the Matsubara frequency 
ujn = TTT{2n + 1)) is given by [H], [H] 



H(*u;) = 2^ 



dep 
dp 



G(p, iUJ„ + iU!m)G{p, iUJm), 



where G(p,«aj„) — (ia;„ — — E(p,jaj„))~-^ is the elec- 
tron Green's function. In the case when the inter-band 
gap Eg is the largest scale in the problem, i.e. when 
Wb < < Eg, in this region one has approximately 
ImH(u;) « and the limit fie — 5" oo in Eg. ([29]) is jus- 
tified. In that case one has H(0) « J^" (lmll{uj)/uj)duj 
which gives the approximate formula for W{Clc,T) 



5:/rf3,^i(.)|-^.(x)^5:^ct, 



p,cr 



2m^ 



(26) 



The partial sum rule is related to the energetics solely 
in the conduction (valence) band which is described by 
the Hamiltonian of the conduction (valence) band elec- 
trons 



(27) 



Hy contains the band-energy with the dispersion ep 
(5p = ~ f-) s-nd the effective Coulomb interaction of 
the valence electrons T4,c- In this case the partial sum- 
rule in the normal state reads [t^ (for general form of 



ne 
2V 



T'V.pjH^ 
Tod 



(28) 



where the number operator fiy^p — c|jg.Cpo-; 1/TOp = 
d'^e^/dpl. is the momentum dependent reciprocal mass 
and V is volume. In practice, measurements are per- 
formed up to finite frequency and the integration over 
Lo goes up to some cutoff frequency f2c (of the order of 
the band plasma frequency). In this case the restricted 
sum-rule has the form 



W{n,,T)^ / <„(w)(ia; 
Jo 



/toH(w) 



du. 



(29) 



where K''' is the diamagnetic Kernel given by Eg. ([50)) 
below and H(w) is the paramagnetic (current-current) 
response function. In the perturbation theory without 



9p2 



(30) 



where ' 



•'V,p, 



is the quasi-particle distribution func- 



tion in the interacting system. Note, that W{VLc,T) is 
cutoff dependent while K'^ in Eg. ipU)) does not depend 
on the cutoff energy J7c. So, one should be careful not 
to over-interpret the experimental results in cuprates by 
this formula. In that respect the best way is to calculate 
W{Vlc,T) by using the exact result in Eg. ip^ which ap- 
parently depends on ilc- However, Eg. ([50]) is useful for 
appropriately chosen ilc, since it allows us to obtain semi- 
quantitative results. In most papers related to the re- 
stricted sum-rule in HTSC cuprates it was assumed, due 
to simplicity, the tight-binding model with nearest neigh- 
bors (n.n.) with the energy ep = —2t{cospxa -\- cospyo) 
which gives 1/TOp — —2ta^ cosp^a. It is straightforward 
to show that in this case Eg. ([50)) is reduced to a simpler 
one 



2V 



-Ty), 



(31) 



where {Ty)H^ = X]p ^p("'w)^ft, the average kinetic en- 
ergy of the band electrons, a is the Cu — Cu lattice dis- 
tance and V is the volume of the system. In this approxi- 
mation W{Vlc, T) is a direct measure of the average band 
(kinetic) energy. In the superconducting state the partial 
band sum-rule reads 



2V 



'Juj)duj 



-Ty) 



(32) 
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In order to introduce the reader to (the complexity of) 
the problem of the T-dependence of W{Qc, T), let us con- 
sider the electronic system in the normal state and in 
absence of quasi-particle interaction. In that case one 
has Up = /p (/p is the Fermi distribution function) and 
Wn{^c,T) increases with the decrease temperature, i.e. 
Wn{n„T) = Wn{Q) - P^T^ where /3f, ~ \/Wb. To this 
end, let us mention in advance that the experimental 
value /Jgxp is much larger than /J^, i.e. /J^^p ^ Pb thus 
telling us that the simple Sommerfeld-like smearing of /p 
by the temperature effects cannot explain quantitatively 
the T-dependence of W^(ric, T). We stress that the smear- 
ing of /p by temperature lowers the spectral weight com- 
pared to that aiT = Q K, i.e. Wn{nc,T) < Wn{^c,0). 
In that respect it is not surprising that there is a low- 
ering of Ws{i^c,T) in the BCS superconducting state, 
Wf^^in^T < Tc) < WniflcT < Tc) since at low 
temperatures /p is smeared mainly due to the super- 
conducting gap, i.e. /p = [1 - i^p/Ep)thiEp/2T)]/2, 

IJ. The maximal decrease of 



Ep = ^Cp + A^ ^p - ep 
Wsi^cT) is at T = 0. 

Let us enumerate and discuss the main experimental 
results for W{Qc,T) in WTSC cwpiates: 1. in the norma/ 
state (T > Tc) of most cuprates, one has Wn{^c,T) = 
W„(0) - /3,,r2 with /33,p > /3b, i.e. VF„(17„T) is in- 
creasing by decreasing T, even at T below T* - the tem- 
perature for the opening of the pseudogap. The increase 
of Wn{^c,T) from room temperature down to Tc is no 
more than 5 %. 2. In the superconducting state (T < Tc) 
of some unde rdop ed and optimally doped Bi-2212 com- 
pounds [ig (and underdoped Bi-2212 films [Z9|) 
there is an effective increase of Ws{^c, T) with respect to 
that in the normal state, i.e. Ws{^c,T) > Wn{^c,T) for 
T < Tc- This is a non-BCS behavior which is shown in 
Fig. ini Note, that in the tight binding model the effec- 
tive band (kinetic) energy {Ty) is negative ((T„) < 0) and 
in the standard BCS case Eq.(l32]) gives that Ws{T < Tc) 
decreases due to the increase of (T„). Therefore the ex- 
perimental increase of Ws{T < Tc) by decreasing T is 
called the non-BCS behavior. The latter means a lower- 
ing of the kinetic energy (T„) which is frequently inter- 
preted to be due either to strong correlations or to a Bose- 
Einstein condensation (BEC) of the preformed tig htly 
bound Cooper pair-bosons, for instance bipolarons '8C|. 
It is known that in the latter case the kinetic energy of 
bosons is decreased below the BEC critical temperature 
Tc- In [m it is speculated that the latter case might be 
realized in in underdoped cuprates. 

However, in some optimally doped and in most over- 
doped cuprates, there is a decrease of Ws{^c,T) at 
T <Tc {W,inc,T) < Wni^cT)) which is the BCS-like 
behavior |82| as it is seen in Fig. [TU] 

We stress that the non-BCS behavior of Ws{flc,T) 
for underdoped (and in some optimally doped) systems 
was obtained by assuming that fie ~ (1 ~ 1-2) eV. 
However, in Ref. these results were questioned 

and the conventional BCS-like behavior was observed 
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FIG. 9. Measured spectral weight Ws{ilc,T){r^ Ai+n in 
figures) for ujc ~ 1.25eV in two underdoped Bi2212 (with 
Tc^88 K and n = 66 K). From [7|. 



{Ws{nc, T) < Wni^c, T)) in the optimally doped YBCO 
and slightly underdoped Bi-2212 by using larger cutoff 
energy ftc — 1.5 eV. This discussion demonstrates how 
risky is to make definite conclusions on some fundamen- 
tal physics based on the parameter (such as the cutoff 
energy flc] dependent analysis. Although the results ob- 
tained in 6£?] looks very trustfully, it is fair to say that 
the issue of the reduced spectral weight in the super- 
conducting state of cuprates is still unsettled and un- 
der dispute. In overdoped Bi-2212 films, the BCS-like 
behavior Ws{^c,T) < Wn{^c,T) was observed, while 
in LSCO it was found that Ws{^c,T) w const, i.e. 
Wsifl^T <Tc)^Wninc,Tc). 

The first question is - how to explain the strong tem- 
perature dependence of W^(r2c,r) in the normal state? 
In [i^l, [11] W{T) is explained solely in the framework 
of the EPI physics where the EPI relaxation 7ep(T) 
plays the main role in the T-dependence of W{ilc,T). 
The main theoretical results of [l^l, [l^ are the follow- 
ing: the calculations of W{T) based on exact Eg. ([30)1 
give that for ftc ^ (the Debye energy) the differ- 
ence in spectral weights of the normal and supercon- 
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FIG. 10. (top) Spectral weight W„(r2c,r) of the overdoped 
Bi2212 for flc = leV. Closed symbols - normal state. Open 
symbols - superconducting state, (bottom) Change of the 
kinetic energy AEtin = Ekin,s — Eki„,N in meV per Cu site 
vs the charge p per Cu with respect to the optimal value Popt- 
From 1831. 



ductmg state is small, i.e. Wn{^c,T) « Ws{^c,T) since 
Wniflc,T)-W,{nc,T) - A^/nl. (2) In the case of large 
r2c the calculations based on Eg. ([30)) gives 
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In the case of EPI, one has 7 = 7ep(^) + 7 



/ imp 



(33) 



where 



7^p(T) = dza2(z)F(z)coth(z/2T). It turns out that 
for a'^{ijj)F{ijj) shown in Fig. HI one obtains: (i) 7ep(T) ~ 
in the temperature interval 100 K <T < 200 K as it 
is seen in Fig. [lT]for the T-dependence of WiQc, T) [H, 
[1^; (ii) the second term in Eg. (1551) is much larger than 
the last one (the Sommerfeld-like term). For the EPI 
coupling constant Aep.tr = 1.5 one obtains rather good 



FIG. 11. Spectral weight W(Q.c,T) in the normal state for 
Einstein phonons with Q.e = 200 K (full triangles) and Q.e = 
400 K (open circles, left axis). Dashed lines is asymptotic. 
From [2a|. 



agreement between the theory in [23|) |2a| and experi- 
ments [6§|, [T^I, [ill. At lower temperatures, 7ep(r) de- 
viates from the behavior and this deviation depends 
on the structure of the spectrum in a'^{uj)F{uj). It is 
seen in Fig. [TT] that for a softer Einstein spectrum (with 
U.E = 200 K), W{^c,T) lies above the curve with the 
asymptotic behavior, while the curve with a harder 
phononic spectrum (with rig = 400 K) lies below it. 
This result means that different behavior of W{^lc, T) in 
the superconducting state of cuprates for different dop- 
ing might be simply related to different contributions of 
low and high frequency phonons. We stress that such a 
behavior of W{nc,T) was observed in experiments [69j . 
fr^ . [78| . To summarize, the above analysis demon- 
strates that the theory based on EPI is able to explain in 
a satisfactory way the strange temperature behavior of 
W{^lc, T) above and below Tc in systems at and near op- 
timal doping and that there is no need to invoke exotic 
scattering mechanisms. In that respect we would like 
to stress that at present we still do not have fully mi- 
croscopic theory which comprises strong correlations and 
EPI and which is able to predict the complex behavior 
of W{T) as a function of T, doping etc. 

4. a'^{uj)F(uj) and the in-plane resistivity Pab{T) 
The temperature dependence of the in-plane resistivity 
Pab{T) in cuprates is a direct consequence of the quasi-2Z? 
motion of quasi-particles and of the inelastic scattering 
which they experience. At present, there is no consensus 
on the origin of the linear temperature dependence of the 
in-plane resistivity /f)jjj(T) in the normal state. Our in- 
tention is not to discuss this problem, but only to demon- 
strate that the EPI spectral function a^(w)i^(w), which 
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FIG. 12. (top) Calculated resistivity p{T) for the EPI spectral 
function ajr(ci;)F(cj) in ^8^. (bottom) Measured resistivity in 
a(x)- and b(y)-crystal direction of YBCO Isil] and calculated 
Bloch-Griineisen curve (points) for X^^ = 1, [Qa ]. 



is obtained from tunnelling experiments in cuprates (see 
Subsection D), is able to explain the temperature de- 
pendence of Pab{T) in optimally doped YBCO. In the 
Boltzmann theory Pab{T) is given by 



Pab{T) 



47r 



(34) 



T 



duj- 



sinh^(a;/2T) 



MF{^). (35) 



where Pi^p is the impurity scattering. Since p{T) = 
l/a{u) = 0,T) and having in mind that the dynamical 
conductivity a{uj,T) in YBCO (at and near the opti- 
mal doping) is satisfactory explained by the EPI scatter- 
ing, then it is to expect that Pab{T) is also dominated 
by EPI in some temperature region T > T^. This is 
indeed confirmed in the optimally doped YBCO, where 
Pirap is chosen appropriately and the spectral function 
<^Aoj)F{uj) is taken from the tunnelling experiments in 
[27| . The very good agreement with the experimen- 
tal results [is'] is shown in Fig. [T^l We stress that 
in the case of EPI there is always a temperature re- 
gion where "ftAT) ^ T for T > aQo, a < 1 de- 
pending on the shape of atj.{uj)F{uj) (for the simple De- 
bye spectrum a ~ 0.2). In the linear regime one has 

P{T) - P^rnp + STT^XepMkBT/fuvl) = p„„p + p'T. 

There is an experimental constraint on Xtr, i.e. Xtr ~ 
0.25ujpi{eV)p' {p,^ cm/K), which imposes a limit on it. 
For instance, for cj^; « (2 — 3) eV 74] and p' « 0.6 in 



the oriented YBCO films and p' « 0.3 — 0.4 in single 
crystals of BSCO, one obtains Xtr ~ 0.6 — 1.4. In case of 
YBCO sin gle crystals, there is a pronounced anisotropy 
in p„fc(T) [8l which gives p'^{T) = 0.6 pilcm/K and 
p'y{T) — 0.25 p,flcm/K. The function Xtrit^pi) is shown 
in Fig ll3l where the plasma frequency Upi can be calcu- 
lated by LDA-DFT and also extracted from the width 
('^ Wpj) of the Drude peak at small frequencies, where 
^pi = V^^pi ■ We stress that the rather good agreement 
of theoretical and experimental results for PabiT) should 
be not over-interpreted in the sense that the above rather 
simple electron-phonon approach can explain the resis- 
tivity in HTSC for various doping. In fact this in highly 
underdoped systems where PabiT) is very different from 
the behavior in Fig. [121 In this case one should certainly 
take into account polaronic effects strong corre- 

lations, etc. and in these cases the simplified Migdal- 
Eliashberg theory is not sufficient. The above analy- 
sis on the resistivity in optimally doped system demon- 
strates only, that if in Eq. ([35]) one uses the EPI spectral 
function a^(cj)i<"(cj) obtained from the tunnelling experi- 
ments (and optics) one obtains the correct T-dependence 
of p^ f,{T) for some optimally doped cuprates, which is an 
additional evidence for the importance of EPI. 

Concerning the linear (in T) resistivity we would like 
to point out that there is some evidence that linear resis- 
tivity is observed in HTSC cuprates sometimes at tem- 
peratures T < 0.2 9d [stI, In that respect, it was 
shown by [89] that in two-dimensional systems with a 
broad interval of phonon spectra the quasi-linear behav- 
ior of Pab{T) is realized even at T < 0.2 Qd- The quasi- 
linear behavior of the resistivity at T ^ 0.2 O^i has been 
observed in Bi^fSrn o7Prn nn'^)9Cu0 fi l86j inLSCO and 1- 
layer Bi-2201 |7|, [88|, (^1, [9l|, UTan systems with 
rather small Tc ~ 10 K. We would like to emphasize here 
that some of these observations are contradictory. For ex- 
ample, the results obtained in the group of Vedeneev [9]| . 
[9^ show that some samples demonstrate the quasi-linear 
behavior of the resistivity up to T = 10 if but some 
others with approximately the same Tc have the usual 
Bloch-Griineisen type behavior characteristic for EPI. In 
that respect it is very unlikely that the linear resistivity 
up to T = 10 if can be simply explained in the standard 
way by interactions of electrons with some known bosons 
either by phonons or spin- fluctuations (magnons). The 
question, why in some cuprates the linear resistivity is 
observed up to T = 10 if is still a mystery and its expla- 
nation is a challenge for all kinds of the electron-boson 
scattering, not only for EPI. In that respect it is interest- 
ing to mention that the existence of the forward scatter- 
ing peak in EPI (with the width qc <C kp), which is due 
to strong correlations, may give rise to the linear behav- 
ior of p(T) down to very low temperatures T ^ 9£)/30 
[1, [sl, [il] - see more in Part II, Section VII. D item 
6.. We shall argue in Subsection D that if one interprets 
the tunnelling experiments in systems near optimal dop- 
ing [231 -[111 in the framework of the Eliashberg theory 



one obtains large EPI coupling constant Ae 



3.5 



20 




Plasma frequency (eV) 



FIG. 13. Transport EPI spectral function coupling constant in 
YBCO as a function of plasma frequency ujp as derived from 
the experimental slope of resistivity p'{T). \x for p'^{T) = 
0.6^iricm/A' and for p^(T) = 0.25iiQcm/K 84]. Squares 
are LDA values 18511. 



which implies that Xtr ^ (A/3). This nieans that EPI is 
reduced much more in transport properties than in the 
self-energy, which is due to some reasons that shall be 
discussed in Part II. Such a large reduction of Xtr cannot 
be obtained within the LDA-DFT band structure calcu- 
lations which means that Xep and Xtr contain renormal- 
ization which do not enter in the LDA — DFT theory. 
In Part II we shall argue that the strong suppression of 
Xtr n iay have its origin in strong electronic correlations 
[Slj . (Oq and in the long-range Madelung energy [l], [J- 
4. Femtosecond time-resolved optical spectroscopy 
The femtosecond time-resolved optical spectroscopy 
(FTROS) has been developed in the last couple of years 
and applied to HTSC cuprates. In this method a fem- 
tosecond (1 fs = 10^^^ sec) laser pump excites in mate- 
rials electron-hole pairs via interband transitions. These 
hot carriers release their energy via electron-electron 
(with the relaxation time Tee) and electron-phonon scat- 
tering reaching states near the Fermi energy within 
10 — 100 fs - see The typical energy density of 

the laser pump pulses with the wavelength A « 810 nm 
{hLj = 1.5 eV) was around F ~ l^J/cm'^ (the excitation 
fluence F) which produces approximately 3 x 10^" car- 
riers per pulse (by assuming that each photon produces 
Huj/ A carriers, A is the superconducting gap). By mea- 
suring photoinduced changes of the reflectivity in time, 
i.e. AR{t)/ Rq, one can extract information on the relax- 
ation dynamics of the low-laying electronic excitations. 
Since AR{t) relax to equilibrium the fit with exponential 
functions is used 



AR(t) 
Rf) 



fit) 



Ae 



Be 



(36) 



where f{t) = H[t)[l - exp{-t/rej] {H{t) is the 



Heavyside function) describes the finite rise-time. The 
parameters A, B depends on the fiuence F. This 
method was used in studying the superconucting phase 
of La2-xSr^CuOi, with x = 0.1, 0.15 and Te ^ 30 K 
and 38 K respectively [l^. In that case one has A ^ 
for T<Tc and A = ior T > Tc, while the signal B was 
present also at T > Tc. It turns out that the signal A 
is related to the quasi-particle recombination across the 
superconducting gap A(r) and has a relaxation time of 
the order ta > W ps a,t T — 4.5 K. At the so called 
threshold fluence (Ft = 4.2 ± 1.7 fiJ/cm'^ for x = 0.1 
and Ft ~ 5.8 ± 2.3 /cm? for x = 0.15) the vaporiza- 
tion (destroying) of the superconducting phase occurs, 
where the parameter A saturates. This vaporization pro- 
cess takes place at the time scala « 0.8 ps. The 
external fluence is distributed in the sample over the 
excitation volume which is proportional to the optical 
penetration depth Aop(~ 150 nm at A « 810 nm) of 
the pump. The energy densities stored in the excita- 
tion volume at the vaporization threshold for a; = 0.1 
and a: = 0.15 are Up = Fx/Xop = 2.0 ± 0.8 K/Cu and 
2.6 ± 1.0 K/Cu, respectively. The important fact is that 
Up is much larger than the superconducting condensa- 
tion energy which is Ucond ~ 0.12 K/Cu for x = 0.1 
and Ucond ~ 0.3 K/Cu for x = 0.15, i.e. Up > Ucond- 
This means that the energy difference Up — Ucond must be 
stored elsewhere on the time scale r^. The only present 
reservoir which can absorb the difference in energy are 
the bosonic baths of phonons and spin fluctuations. The 
energy required to heat the spin reservoir from T — 4.5 
K to Tc is Usf = ff''Css{T)dT. The measured spe- 
cific heat Csf{T) in La2CuOi (2^ gives very small value 
Usf ~ 0.01 K. In the case of the phonon reservoir on 
obtains Upu = j^" Cph{T)dT = 9 K/Cu for x = 0.1 and 
28 K/Cu for x = 0.15, where Cph is the phonon spe- 
cific heat. Since Ugf <^ Up — Ucond the spin reservoir 
cannot absorb the rest energy Up — Ucond- The situa- 
tion is opposite with phonons since Upu ^ Up — Ucond 
and phonon can absorb the rest energy in the excitation 
volume. The complete vaporization dynamics can be de- 
scribed in the framework of the Rothwarf- Taylor model 
which describes approaching of electrons and phonons to 
quasi-equilibrium on the time scale of 1 ps [99| . We shall 
not go into details but only summarize by quoting the 
conclusion in (ooj . that only phonon-mediated vaporiza- 
tion is consistent with the experiments, thus ruling out 
spin-mediated quasi-particle recombination and pairing 
in HTSC cuprates. The FTROS method teU us that at 
least for non-equilibrium processes EPI is more impor- 
tant than SFI. It gives also some opportunities for ob- 
taining the strength of EPI but at present there is no 
reliable analysis. 

In conclusion^ optics and resistivity measurements in 
the normal state of cuprates give evidence that EPI is 
important while the spin fluctuation scattering is weaker 
than it is believed. However, some important questions 
related to the transport properties remain to be an- 
swered: (i) what are the values of Xtr and uipi ; (ii) what is 
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the reason that Xtr ^ A is reahzed in cuprates; (iii) what 
is the role of Coulomb scattering in 17(0;) and p{T). Later 
on we shall argue that ARPES measurements in cuprates 
give evidence for an appreciable contribution of Coulomb 
scattering at higher frequencies, where 7(w) ~ 7o + 
for u! > wfj^ax with Ac ^ 1. One should stress, that de- 
spite the fact that EPI is suppressed in transport prop- 
erties it is sufficiently strong in the quasi-particle self- 
energy, as it comes out from tunnelling measurements 
discussed below. 



C. ARPES and the EPI self-energy 

The angle-resolved photoemission spectroscopy 
(ARPES) is nowadays one of leading spectroscopy 
methods in the solid state physics [l^l- In some favor- 
able conditions it provides direct information on the 
one-electron removal spectrum in a complex many body 
system. The method involves shining light (photons) 
with energies between Ei = b — 1000 eV on samples and 
by detecting momentum (k)- and energy(w)-distribution 
of the outgoing electrons. The resolution of ARPES 
has been significantly increased in the last decade 
with the energy resolution of AE « 1 — 2 meV (for 
photon energies ~ 20 eV) and angular resolution of 
Ae < 0.2°. On the other side the ARPES method is 
surface sensitive technique, since the average escape 
depth {lesc) of the outgoing electrons is of the order 
of lesc ~ 10 A, depending on the energy of incoming 
photons. Therefore, very good surfaces are needed in 
order that the results to be representative for bulk 
samples. The most reliable studies were done on the 
bilayer Bi2Sr2CaCu20s [Bi — 2212) and its single 
layer counterpart Bi2Sr2CuOe {Bi2201), since these 
materials contain weakly coupled BiO planes with the 
longest inter-plane separation in the cuprates. This 
results in a natural cleavage plane making these materi- 
als superior to others in ARPES experiments. After a 
drastic improvement of sample quality in other families 
of HTSC materials, the ARPES technique has became 
an important method in theoretical considerations. The 
ARPES can indirectly give information on the momen- 
tum and energy dependence of the pairing potential. 
Furthermore, the electronic spectrum of the (above 
mentioned) cuprates is highly quasi-2D which allows 
rather unambiguous determination of the initial state 
momentum from the measured final state momentum, 
since the component parallel to the surface is conserved 
in photoemission. In this case, the ARPES probes 
(under some favorable conditions) directly the single 
particle spectral function A(k, w). In the following we 
discuss mainly those ARPES experiments which give 
evidence for the importance of the EPI in cuprates - see 
more in p2| . 

ARPES measures a nonlinear response function of the 
electron system and it is usually analyzed in the so-called 
three-step model, where the total photoemission inten- 



sity Itot (k, w) w /i • /2 • I3 is the product of three inde- 
pendent terms: (1) Ii - describes optical excitation of 
the electron in the bulk; (2) I2 - describes the scatter- 
ing probability of the travelling electrons; (2) I3 - the 
transmission probability through the surface potential 
barrier. The central quantity in the three-step model 
is /i(k, cj) and it turns out that for k = ky it can be 
written in the form Ii{]i,uj) ~ /o(k, u)/((u;)A(k, w) [l2| 
with /o(k, u) ~| {ipj: I pA ji/ij) p and the quasi-particle 
spectral function A(k, cj) = — ImG'(k, a;)/7r 

1 /mi;(k,w) 
^ ^ "^[w-e(k) -i?eS(k,cj)]2 -f /mS2(k,a;)' 

(37) 

Here, {ip j: \ p ■ A \ip^) is the dipole matrix element which 
depends on k, polarization and energy E.^ of the incom- 
ing photons. The knowledge of the matrix element is of 
a great importance and its calculation from first prin- 
ciples was done in |lOOj |. /(w) is the Fermi function, 
G and S = ReH + ilmT, are the quasi-particle Green's 
function and the self-energy, respectively. We summarize 
and comment here some important ARPES results which 
were obtained in the last several years and which confirm 
the existence of the Fermi surface and importance of EPI 
in the quasi-particle scattering [l^ . 
ARPES in the normal state 

(Nl) There is well defined Fermi surface in the metal- 
lic state of optimally and near optimally doped cuprates 
with the topology predicted by the LDA-DFT. However, 
the bands are narrower than LDA-DFT predicts which 
points to a strong quasi-particle renormalization. (N2) 
The spectral lines are broad with |ImE(k, w) |~ w (or 
~ T for T > Lu) which tells us that the quasi-particle 
liquid is a non-canonical Fermi liquid for larger values of 
T, oj. (N3) There is a bilayer band splitting in Bi - 2212 
(at least in the over-doped state) what is also predicted 
by LDA-DFT. In the case when the coherent hopping 
t^ between two layers in the bilayer dominates, then the 
anti-bonding and bonding bands ^j^'^ — zt t^ with 
t^ = [i-'-(cos2 — cos^ ky) + ...] have been observed. It 
is worth to mention that the previous experiments did 
not show this splitting what was one of the reasons for 
various speculations on some exotic electronic scattering 
and non-Fermi liquid scenarios. (N4) In the under-doped 
cuprates and at temperatures Tc < T < T* there is a d- 
wave like pseudogap Apg(k) ~ Apgfi{coskx — cos ky) in 
the quasi-particle spectrum where Apg^ increases by low- 
ering doping. We stress that the pseudogap phenomenon 
is not well understood at present and since we are inter- 
ested in systems near optimal doping where the pseudo- 
gap phenomena are absent or much less pronounced we 
shall not discuss this problem here. Its origin can be due 
to a precursor superconductivity or due to a competing 
order, such as spin- or charge-density wave, strong corre- 
lations, etc. (N5) The ARPES self-energy gives evidence 
that EPI interaction is rather strong. The arguments for 
latter statement are the following: (i) at T > Tc there are 
kinks in the quasi-particle dispersion a;(^]j.) in the nodal 
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FIG. 14. (top) Quasi-particle dispersion of Bi2212, Bj2201 
and LSCO along the nodal direction, plotted vs the momen- 
tum k for (a) — (c) different doping, and (d) — (e) different T; 
black arrows indicate the kink energy; the red arrow indicates 
the energy of the q = (vr, 0) oxygen stretching phonon mode; 
inset of (e)- T-dependent E' for optimally doped _Bi2212; (/) - 
doping dependence of the effective coupling constant A' along 
(0,0) - (7r,7r) for the different HTSC oxides. From Ref. 
(bottom) Quasi-particle dispersion E{k) in the normal state 
(al, bl, c), at 107 K and 115 K, along various directions (j) 
around the anti-nodal point. The kink at i? = 40meV is 
shown by the horizontal arrow. (a2 and b2) is E{k) in the su- 
perconducting state at 10 K with the shifted kink to 70meV. 
(d) kink positions as a function of in the anti-nodal region. 
From Ref. 101]. 



direction (along the (0, 0) 



teristic phonon energy ld'^!'^^ 



(tt, tt) line) at the charac- 
^ph (60 - 70) meV 0, see 
Fig. I14[ and near the anti-nodal point (tt, 0) at 40 meV 
[Toil - see Fig. M 

(ii) This kink structure is observed in a variety 
of the hole-doped cuprates such as LSCO, Bi — 
2212, Si2201, rZ2201 {Tl2Ba2CuOfi), Na - CCOC 
{Ca2-xNaxCu02Cl2)- These kinks exist also above Tc 
what excludes the scenario with the magnetic resonance 
peak in JmXg(Q,a;). Moreover, since the tunnelling and 
magnetic neutron scattering measurements give small 
SFI coupling constant gsf < 0.2 eV then the kinks can 



FIG. 15. (a) Effective real self-energy for the non- 
superconducting La2-xSrxCuOi, x = 0.03. Extracted 
alff{u))F{uj) is in the inset, (b) Top: the total MDC width - 
open circles. Bottom: the EPI contribution shows saturation, 
impurity contribution - dotte d black line. The residual part 
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not be due to SFI. {Hi) The position of the nodal kink 
is practically doping independent which points towards 
phonons as the scattering and pairing boson. (N6) The 
quasi-particles (holes) at and near the nodal-point kjv 
couple practically to a rather broad spectrum of phonons 
since at least three group of phonons were extracted 
in the bosonic spectral function q^F(\lisi ^ui) fro m t he 



ARPES effective self-energy in La2~xSrxCuOi \1Q1 

Fig.im 

This result is in a qualitative a gree ment with numer- 
ous tunnelling measurements 27|-[32l| which apparently 
demonstrate that the very broad spectrum of phonons 
couples with holes without preferring any particular 
phonons - see discussion below. (N7) Recent ARPES 
measurements in B2212 6lj show very different slope 
dijj/dS^-^ of the quasi-particle energy uj{£,-^) at small | 
1^ LUph and at large energies | £y |^ ujph - see 
Fig. [161 The theoretical analysis [l03| of these results 
gives the total coupling constant = + -^f w 3, 
and for the EPI coupling A^ ~ 2 while the Coulomb 

coupling (SFI is a part of it) is Af « 1 jl03{ - see 
Fig. [ini (Note, that the upper index Z in the coupling 
constants means the quasi-particle renormalization in the 
normal part of the self-energy.) To this end let us men- 
tion some confusion which is related to the value of the 
EPI co uplin g cons tant extracted from ARPES. Namely, 
in [12, UQI], [Toil the EPI self-energy was obtained by 
subtracting the high energy slope of the quasi-particle 
spectrum i^(^fc) at a; ~ 0.3 eV . The latter is appar- 
ently due to the Coulomb interaction. Although the po- 
sition of the low-energy kink is not affected by this pro- 
cedure (if ijJ^^^ ^ i^c), this subtraction procedure gives 
in fact an effective EPI self-energy Egjy(k, w) and the 
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(as we alre ady ar gued ) most probably of the phononic 
origin [l^], |104| . [l05j |. while the second kink at higher 
energy uj2 ~ <^c ~ 350 meV is due to the Coulom b int er- 
action. However, the important results in Ref. 'IOG] is 
that the slopes of i?eE°'^P(k, w) at low (w < uj^^^) and 

high energies (Wp^^'' < w < Wc) are very different. The 
low-energy and high-energy slope near the nodal point are 
shown in Fig. [16] schematically (thin lines). From Fig. [TBI 
it is obvious that EPI prevails at low energies uj < u/^]'^^ ■ 
More precisely digitalization of i?eS°''P(k, w) in the in- 



terval w 



high 
ph 



< uj < 0.4 eV gives the Coulomb couphng 

< 



Af ~ 1.1 while the same procedure at 20 meV ~ cj'"™ 



U! < LU 



high 
ph 



^ph 

50 — 70meV gives the total coupling con- 



0.0 Q2 



oKeV) 



stant (A2 =)A^ = Af + A^ « 3.2 and the EPI couphng 



the 



FIG. 16. Fig.4b from [lM|: ReT,{uj) measured in 612212 
(thin line) and model f?eE(a;) (bold line) obtained In [lOd ]. 
The three thin lines (Ai, A2, A3) are the slopes of _ReE(i:i;) In 
different energy regions - see the text. 



effective coupling constant A^^ j (k) only. We demon- 
strate below that the A^ f (k) is smaller than the real 

EPI coupling constant A^(k). The total self-energy is 
T.{k,uj) = E'=P(k,tj) -I- S]'=(k,cj) where is the contri- 
bution due to the Coulomb interaction. At very low en- 
ergies LU <ti LOc one has usually E'^(k, w) = — Af(k)w, 
where Wc(^ 1 eV) is the characteristic Coulomb ener- 
gies and Aj, is the Coulomb coupling constant. The 
quasi-particle spectrum a;(k) is determined from the con- 
dition w - f(k) - i?e[S^P(k,w) + I]'=(k,w)] = where 
f (k) is the bare band structure energy. At low energies 



UJ < LU 



ph 



ujc it can be rewritten in the form 



^-r"(k)-i?eS:^/k,c.)=0, 



with ^'■'"(k) = [1 + Af (k)]-ie(k) and 



i?eS:^y(k,c.) 



Since at very low energies w ^ 
i?eE^P(k,tj) = 



-Agp(k)a; 



-Af 



gyy(k)a;, then the real coupling constant is related 



to the effective one by Agp(k) = [1 



^fp,e//(k)- At higher energies uj'^^'' 
effects are suppressed and I]'^^(k, w) stops growing, one 
has i?eS(k,w) « i?eE'=P(k,w) - Af (k)u;. The measured 
i?eS®''P(k, w) at T = 10 iiT near and shghtly away from 
the nodal point in the optimally doped Bi-2212 with 
Tc^91 K [I06l | is shown in Fig. [Hj 

It is seen that i?eS]™P(k, w) has two kinks - the first 



one at low energy oji 



high 
ph 



constant X^p{= Kp,high) ~ 2.1 > 2X^j,^^ffy 
EPI coupling is at least twice larger than the effective 
EPI coupling constant obtained in th e previous analy- 
sis of ARPES results [loi|, [loi|. This estimation 
tells us that at (and near) the nodal point, the EPI in- 
teraction dominates in the quasi-particle scattering at 
low energies since A^(w 2.1) 2A^ > 2Afj, while at 
large energies (compared to LUph), the Coulomb interac- 
tion with Af « 1.1 dominates. We point out that EPI 
near the anti-nodal point can be even larger than in the 
nodal point, mostly due to the higher density of states 
near the anti-nodal point. (N8) Recent ARPES spectra 
in the optimally dop ed Bi — 2212 near the nodal and 
anti-nodal point [l05l [ show a low energy isotope effect in 
Re'S'^^P (k, Lu) , which can be well described in t he fr ame- 
work of the Migdal-Eliashberg theory for EPI jl07l |. At 
higher energies w > ujph the obtained in [10^ very pro- 
nounced isotope effect can not b e explained by the simple 
Migdal-Eliashberg theory (l07j |. However, there are con- 
troversies with the strength of the high-energy isotope 
effec t since it was not confirmed in other measurements 
[l08j |. |l09l | - see the discussion in Section F.2 related to 
the isotope effects in HTSC cuprates. (N9) ARPES ex- 
periments on Ca2Cu02Cl2 give strong evidence for the 
formation of small polarons in undoped cuprates which 
are due to phonons and strong EPI, while by doping 
quas i-particles appear and there are no small polarons 
[llOj |. Namely, in ,11(| a broad peak around —0.8 eV is 
observed at the top of the band (k = (7r/2, 7r/2)) with the 
dispersion similar to that predicted by the i — J model - 
see Fig. [T71 

However, the peak in Fig. [17] (left) is of Gaussian shape 
and can be described only by coupling to bosons, i.e. this 
peak is a boson side band - see more in Q and refer- 
ences therein. The theory based on the t-J model (in the 
antiferromagnetic state of the undoped compound) by 
including coupling to several (half-breathi ng, a pical oxy- 
gen, low- lying) phonons, which is given in [lll| . explains 
successfully this broad peak of the boson side band by 
the formation of small polarons due to the EPI coupling 
(Aep ~ 1.2). Note that this value of Aep is for the po- 
50 — 70 meV which is laron at the bottom of the band while in the case where 



i?el]r,,(k,c.) 



1 + Af(k) 



(38) 



(39) 



ph 1 "-"^^"^ 
and ReYfJ'j j{k.,uj) 



has 



KinXtp^effi^) > 

< oj < LUr the EPI 
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FIG. 17. (left-b) The ARPES spectrum of undoped 
Ca2Cu02Cl2 at k = (7r/2, 7r/2). Gaussian shape - solid line, 
Lorentzian shape - dashed line, (right-c) Dispersion of the po- 
laronic band - A and of the quasi-particle band - B along the 
nodal direction. Horizontal lines are the chemical potentials 
for a large number of samples. From [llO|] . 



the Fermi surface exists (in doped systems) this coupling 
is even larger due to the larg er density of states at the 
Fermi surface [m|. In [m| it was stressed that even 
when the electron-magnon interaction is stronger than 
EPI the polarons in the undoped systems are formed due 
to EPI. The latter mechanism involves excitation of many 
phonons at the lattice site (where the hole is seating), 
while it is possible to excite only one magnon at the 
given site. (NIO) Recent soft x-ray ARPES measure- 
ments on the electron doped HTSC Ndi,^c^CeQ,i^CuOi 
[ll2l] . and Sm(2-x)Ce^CuOi {x = 1, 0.15, 0.18), 
-^rfi.85Ceo.i5CM04, Eui,%^CeQ,izCuOi jll3| show kink 
at energies 50 — 70 meV in the quasi-particle dispersion 
relation along both nodal and antinodal directions as it 
is shown in Fig. [iHl 

It is seen from this figure that the effective EPI cou- 
pling constant Aep,e//(< Aep) is isotropic and Xep,eff ~ 
0.8 — 1. It seems that the kink in the electron-doped 
cuprates is due solely to EPI and in that respect the sit- 
uation is similar to the one in the hole-doped cuprates. 
ARPES results in the superconducting state 
(SI) There is an anisotropic superconducting gap in 
most HTSC compounds which is predominately d- 
wave like, i.e. A(k) « Ao(cos fc^; — cos fc^) with 2Aq/Tc ~ 
5 — 6. (S2) The particle-hole coherence in the super- 
conducting state which is expected for the BCS-like the- 
ory of superconductivity has been observed first in 1^ 
and confirmed with better resolution in JJLSj, where the 
particle-hole mixing is clearly seen in the electron and 
hole quasi-particle dispersion. To remaind the reader, the 
excited Bogoliubov-Valatin quasi-particles (ak,±) with 



energies i?^* 



lAkl are a mixture of electron 
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FIG. 18. NCCO electron-doped: (a) /mE(aj) measured in 
the nodal point. Curves are offsets by 50 meV for clarity. 
The change of the slope in the last bottom curve is at the 
phonon energy, (b) /mE(a;) for the antinodal direction with 
30 meV offset, (c) Experimental phonon dispersion of the 
bond stretching modes, (d) Estimated Ajj^y. from /mS(a;). 
From [HI- 



cik.- = MkC-k^ + ^^kcj(.^ where the coherence factors Mk, 
Vk are given by |Mk|^ = 1 - \vk\^ = (1 + ^k/-E^k)/2. Note, 
that |uk|^ + \v-k\^ = 1 what is exactly observed, together 
with d-wave pairing A(k) = Ao(cosA:a; — cosfcj,), in ex- 
periments JT5J. This is very important result since it 
proves that the pairing in HTSC cuprates is of the BCS 
type and not exotic one as was speculated long time af- 
ter the discovery of HTSC cuprates. (S3) The kink at 
(60 — 70) meV in the quasi-particle energy around the 
nodal point is not shifted (in energy) while the antin- 
odal kink at ui^p^'' ~ 40 meV is shifted (in energy) in 
the superconducting state by Ao(= (25 — 'iQ)meV), i.e. 

^^ph ^^ph + ^0 = (65 - 7Q)meV To remind the 
reader, in the standard Eliashberg theory the kink in the 
normal state at w = should be shifted in the super- 
conducting state to ujph -f Aq at all points at the Fermi 
surface. This puzzling result (that the quasi-particle en- 
ergy around the nodal point is not shifted in the super- 
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conducting state) might be a smoking gun result since 
it makes an additional constraint on the quasi-particle 
interaction in cupr ates . Until now there is only one plau- 
sible explanation of this non-shift puzzle which is 
based on an assumption of the forward scattering peak 
(FSP) in EPI - see more in Part 11. The FSP in EPI means 
that electrons scatter into a narrow region {q < Qc <^ kp) 
around the initial point in the fc-space, so that at the most 
part of the Fermi surface there is practically no mixing of 
states with different signs of the order parameter A(k). 
In that case the EPI bosonic spectral function (which is 
defined in Appendix A) a'^F(k,k',n) « F {ip , ip' , fl) 
{if is the angle on the Fermi surface) has a pronounced 
forward scattering peak (at dip = ip — ip' = 0) due 
to strong correlations - see Part II. Its width S(p^ is 
narrow, i.e Sip^ <^ 2t: and the angle integration goes 
over the region 5ip^ around the point ip. In that case 
the kink is shifted (approximately) by t he l ocal gap 
A{ip) = AniaxCOs2(/3 - for more details see |ll6l |. As the 
consequence, the anti-nodal kink is shifted by the max- 
imal gap, i.e. \A{ipj^j^ « tt/2)\ = A^ax while the nodal 
gap is practically unshifted since l Ajtp ^M ~ 7r/4)| w 0. 
(S4) The recent ARPES spectra |ll7l | on an undoped 
single crystalline 4-layered cuprate with apical fluorine 
(F), Ba2Ca3Cu40sF2 (F0234) give rather convincing ev- 
idence against the SFI mechanism of pairing - see Fig. [T21 

First, F0234 is not a Mott insulator - as expected from 
valence charge counting which puts Cu valence as 2+, 
but it is a superconductor with Tc = 60 K. Moreover, 
the ARPES data [117{ reveal at least two metallic Fermi- 
surface sheets with corresponding volumes equally below 
and above half-filling - see Fig. [20l 

Second, one of the Fermi-surfaces is due to the electron- 
like (N) band (with 20 ± 6% electron-doping) and the 
other one due to the hole-like (P) band (with 20 ± 8 
% hole-doping) and their splitting along the nodal di- 
rection is significant and cannot be explained by the 
LDA-DFT calculations [118]. This electron and hole self- 
doping of inner and outer layers is in an appreciable con- 
trast to other multilayered cuprates where there is only 
hole self-doping. For instance, in HgBa2Ca„Cu„+i02n+2 
{n = 2,3) and (Cu,C)Ba2Ca„Cu„+i03„+2 {n = 2,3,4), 
the inner Cu02 layers are less hole-doped than outer lay- 
ers. It turns out, unexpectedly, that the superconducting 
gap on the N-band Fermi-surface is significantly larger 
than on the P-one, where in Ba2Ca^Cu40fiF2 the ra- 
tio is anomalous (A^r/Ap) « 2 and A^r is an order 
of magnitude larger than in the electron-doped cuprate 
Nd2~xCexCu04^. Third, the N-band Fermi-surface is 
rather far from the antinodal point at (tt, 0). This is very 
important result which means that the antiferromagnetic 
spin fluctuations with the AF wave- vector Q = (tt, tt), as 
well as the van Hove singularity, are not dominant in the 
pairing in the TV-band. To remind the reader, the SFI 
scenario assumes that the pairing is due to spin fluctu- 
ations with the wave-vector Q (and near it) which con- 
nects two anti-nodal points which are near the van Hove 
singularity at the hole-surface (at (tt, 0) and (0, tt)) giv- 
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FIG. 19. Crystal structure of Ba2CaiCuiOs,{05Fx^5)2. 
There are four Cu02 layers in a unit cell with the outer having 
apical F atoms. CRL - charge reservoir layer; SCL/OP - su- 
perconducting layer/out er pl ane; SCL/IP - superconducting 
layer/inner plane. From [llTl ]. 



ing rise to large density of states. This is apparently 
not the case for the TV-band Fermi-surface - see Fig. [20l 
The ARPES data give further that there is a kink at 
^ 85 meV in the quasi-particle dispersion of both bands, 
while the kink in the iV-band is stronger than that in the 
P-band. This result, together with the anomalous ratio 
(Aat/Ap) 2, strongly disfavors SFI as a pairing mech- 
anism. (S5) Despite the presence of significant elastic 
quasi-particle scattering in a number of samples of opti- 
mally doped Bi-2212, there are dramatic sharpening of 
the spectral function near the anti-nod al po int (7r,0) at 
T < Tc (in the superconducting state) |ll9j |. This effect 
can be explained by assuming that the small g-scattering 
(the forward scattering peak) dominates in the e lasti c 
impurity scattering as it is pointed in [5l|, [o^, |l20j |. 
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FIG. 20. (a) Fermi surface (FS) contours from two samples 
of F0234. - electron- like; P - hole-like. Bold arrow is (tt, tt) 
scattering vector. Angle 9 defines the horizontal axis in (b). 
(b) Leading gap edge along k-space angle from the two FS 
contours. From [llTll. 



[l2l|. As a result, one finds that the impurity scatter- 
ing rate in the superconducting state is almost zero, i.e. 
7»mp(k,t^) = 7«(k,t^) + 7a(k,a;) w for | w |< Aq for 
any kind of pairing (s- p- c?-wave etc.) since the normal 
(7„) and the anomalous (7^) scattering rates compen- 
sate each other. This collapse of th e elastic scattering 
rate is elaborated in details in |ll9l | and it is a conse- 
quence of the Anderson-like theorem for unconventional 
superconductors which is du e to the d ominance of the 
small g-scattering [5l|, [ol], |l2Cll |- jl2l| . In such a case 
d-wave pairing is weakly unaffected by impurities and 
there is small reduction in Tc [121], [122]. The physics 
behind this result is rather simple. The small g-scattering 
(usually called forward scattering) means that electrons 
scatter into a small region in the k-space, so that at the 
most part of the Fermi surface there is no mixing of states 
with different signs of the order parameter A(k). In such 
a way the detrimental effect of impurities on d-wave pair- 
ing is significantly reduced. This result points to the im- 
portance of strong correlations in the renormalization of 
the impurity scattering too - see discussion in Part II. 

In conclusion, in order to explain ARPES results in 
cuprates it is necessary to take into account: (1) the 
electron-phonon interaction (EPI) since it dominates in 
the quasi-particle scattering in the energy region impor- 
tant for pairing; (2) effects of elastic nonmagnetic im- 
purities with the forward scattering peak (FSP) due to 
strong correlations; (3) the Coulomb interaction which 
dominates at higher energies uj > ujph- In this respect, 
the presence of ARPES kinks and the knee-like shape of 
the T-dependence of the spectral width are important 
constraints on the scattering and pairing mechanism in 
HTSC cuprates. 



D. Tunnelling spectroscopy and spectral function 

By measuring current-voltage I — V characteristics 
in NIS (normal metal-insulator-superconductor) tun- 
nelling junctions with large tunnelling barrier one ob- 
tains from tunnelling conductance GNsiV) = dl/dV the 
so called tunnelling density of states in superconductors 
Nt{uj)- Moreover, by measuring of GNsiV) at voltages 
eV > A it is possible to determine the Eliashberg spectral 
function a'^F{ui) and finally to confirm the phonon mech- 
anism of pairing in LTSC materials. Four tunnelling 
techniques were used in the study of HTSC cuprates: (1) 
vacuum tunnelling by using the STM technique - scan- 
ning tunnelling microscope; (2) point-contact tunnelling; 
(3) break-junction tunnelling; (4) planar-junction tun- 
nelling. Each of these techniques has some advantages 
although in principle the most potential one is the STM 
techniq ue si nce it measures superconducting properties 
locally |l23j |. Since tunnelling measurements probe a 
surface region of the order of superconducting coherence 
length , then this kind of measurements in HTSC ma- 
terials with small coherence length (^^^^ ^ 20 A in the 
a — h plane and 1 — 3 A along the c — axis) de- 

pends strongly on the surface quality and sample prepa- 
ration. Nowadays, many of the material problems in 
HTSC cuprates are understood and as a result consis- 
tent picture of tunnelling features is starting to emerge. 

From tunnelling experiments one obtains the (energy 
dependent) gap function A(ti;) in the superconducting 
state. Since we have already discussed this problem in 
we will only briefiy mention some important result, 
that in most cases GNsiV) has F-shape in all families of 
HTSC hole and electron doped cuprates. The F-shape 
is characteristic for d — wave pairing with gapless spec- 
trum, which is also confirmed in the interference exper- 
iments on hole and electron doped cuprates [i^. Some 
experiments give an [/-shape of Gns{V) which resembles 
s-wave pairing. This controversy is explained to be the 
property of the tunnelling matrix element which filters 
out states with the maximal gap. 

Here we are interested in the bosonic spectral function 
a^Fluj) of HTSC cuprates near optimal doping which 
can be extracted by using tunnelling spectroscopy. We 
inform the reader in advance, that the shape and the 
energy width of a^i^(w), which are extracted from the 
second derivative d^I/dV"^ at voltages above the super- 
conducting gap, in most HTSC cuprates resembles the 
phonon density of states Fph{uj). This result is strong 
evidence for the importance of EPI in the pairing po- 
tential of HTSC cuprates. For instance, plenty of break- 
junctions made from Bi — 2212 single crystals 27] show 
that the peaks (and shoulders) in ~d^I/dV^ (or dips- 
negative peaks in d^I/dV^) coincide with the peaks (and 
shoulders) in the phonon density of states Fph{uj) mea- 
sured by neutron scattering - see Fig. [5TJ 

The tunnelling spectra in Bi-2212 break junctions p?} , 
which are shown in Fig. [^TJ indicates that the spectral 
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FIG. 21. (a) Second derivative of I{V) for a Bi2212 break 
junction in various magnetic fields (from — 20 T). The 
structure of dips (minima) in d?I /dV^ can be compared with 
the phonon density of states F{uj); (b) the spectral functions 
F{ijj) in various magnetic fields. From [27| ]. 

function a^Fiui) is independent of magnetic field, which 
is in contradiction with the theoretical prediction based 
on the SFI pairing mechanism where this function should 
be sensitive to magnetic field. The reported broadening 
of the peaks in a'^F{uj) are partly due to the gapless spec- 
trum of d — wave pairing in HTSC cuprates. Addition- 
ally, the tunnelling density of states Nt{(^) at very low 
T and for w > A show a pronounced gap structure. It 
was found that 2A/Tc = 6.2 - 6.5, where = 74 - 85 if 
and A is some average value of the gap. In order to ob- 
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FIG. 22. The spectral functions a^F{uj) and the calculated 
density of states at OK (upper solid line) obtained from the 
conductance measurements the _Bi(2212) — Au planar junc- 
tions. GPDS - generalized phonon density of states. From 

tain a'^F^Lu) the inverse procedure was used by assuming 
s — wave superconductivity and the effective Coulomb 
parameter fi* « 0.1 [27j. The obtained a'^F{uj) gives 
large EPI coupling constant Xep ~ 2.3. Although this 
analysis [27*1 was done in terms of s — wave pairing, it 
mimics qualitatively the case of d — wave pairing, since 
one expects that d—wave pairing does not change signif- 
icantly the global structure of (Pl/dV^ at eV > A albeit 
introducing a broadening in it - see the physical meaning 
in Appendix A. We point out that the results obtained 
in j27| were reproducible on more than 30 junctions. In 
that respect very important results on slightly overdoped 
Bz— 2212 — Ga^s and on Bi — 2212— Au planar tunnelling 
junctions are obtained in [2^ - see Fig. [22] 

These results show very similar features to those ob- 
tained in (27j on break-junctions. It is worth mentioning 
that several groups [2§|, [l^l, HH have obtained simi- 
lar results for the shape of the spectral function a'^F(Ld) 
from the I — V measurements on various HTSC cuprates 
- see the comparison in Fig. [531 These facts leave no 
much doubts about the importance of the EPI in pair- 
ing mechanism of HTSC cuprates. 

In that respect, tunnelling measurements on slightly 
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FIG. 23. The spectral functions F(uj) from measurements 
of various groups: Vedeneev et al. [27|, G onnelli et al.fsH. 
Miyakawa et al. [2^, Shimada et al. [2^. The generalized 
density of states GPDS for Bi2212 is at the bottom. From 
[H. 
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FIG. 24. The spectral functions o?F{ui) from the tunnelling 
conductance of Bi2Sr2CaCu20s for the positive and the neg- 
ative bias voltages, and the averaged one 28]. The averaged 
one is divided into 13 components. The origin of the ordinate 
is 2, 1,0 and -0.5 from the top down. From 
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overdoped Bi2Sr2CaCu20% [28|, (32| give impressive 
results. The Eliashberg spectral function a^F(aj) of 
this compound was extracted from the measurements 
of d^I/dV^ and by solving the inverse problem - see 
Appendix A. The extracted a^Fiuj) has several peaks 
in broad energy region up to 80 meV as it is seen in 
Figs. [^211^ which coincide rather well with the peaks in 
the phonon density of states Fph{uj) - more precisely the 
generalized phonon density of states GPDS{uj) defined 
in Appendix A. In [sl] numerous peaks, from PI — P13, 
in a^Fiiu) are discerned as shown in Fig. 1241 which corre- 
spond to various groups of phonon modes - laying in (and 
around) these peaks. Moreover, in [2^, [sll the coupling 
constants for these modes are extracted as well as their 
contribution (ATc) to Tc as it is seen in Fig. [55] Note, 
that due to the nonlinearity of the problem the sum of 
ATc is not equal to Tc- 

The next remarkable result is that the ex- 
tracted EPI coupling constant is very large, i.e. 
Aep(= 2jdwa^F{uj)/uj) = A, w 3.5 - see Fig. [551 
It is obvious from Figs. (|24II25I) that almost all phonon 
modes contribute to Agp and Tc, which means that on 
the average each particular phonon mode is not too 
strongly coupled to electrons since < 1.3, thus keeping 



For a better understanding of the the EPI coupling in 
these systems we show in Fig. [26|the total and partial 
density of phononic states. Let us discuss the content of 
the Table in Fig. [25]in more details where it is shown the 
strength of the EPI coupling and the relative contribution 
of different phononic modes to Tc. In Fig. [25] it is seen 
that lower frequency modes from PI — P3, corresponding 
to Cu, Sr and Ca vibrations, are rather strongly coupled 
to electrons (with A^ ^ 1) which give appreciable contri- 
butions to Tc. It is also seen in Fig. [25] that the coupling 
constants Ai of the high-energy phonons (P9 — P13 with 
a; > 70 meV) have Ai <IC 1 and give moderate contri- 
bution to Tc - around 10 %. These results give sohd 
evidence for the importance of the low-energy modes re- 
lated to the change of the Madelung energy in the ionic- 
metallic structure of HTSC cuprates - the idea advocated 
in , and discussed in Part II. If confirmed in other 
HTSC families, these results are in favor of the moder- 
ate oxygen isotope effect in cuprates near the optimal 
doping since the oxygen modes are higher energy modes 
and gives smaller contribution to Tc. We stress that each 
peak PI — P13 in a'^F{ui) corresponds to many modes. 
In order to obtain information on the structure of vibra- 
tions which are strongly involved in pairing, we show in 
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FIG. 25. Table I - Phonon frequency lj, EPI coupling constant 
Xi of the peaks PI — P13 and contribution ATc to Tc of each 
peak in a^_F(ij;)-shown in Fie: l24l obtained from the tunnelling 
conductance of Bi2Sr2CaCu20s- ATc is the decrease in Tc 
when the peak in a'^F(aj) is eliminated. From [3^ . 

Figs. [?71l^ the structure of these vibrations at special 
points in the Brillouin zone. It is seen in Fig. [57]that the 
low-frequency phonons PI — P2 are dominated by Cu, 
Sr, Ca vibrations. 

Further, based on the Table in Fig. [55] one can conclude 
that the P3 modes are stronger coupled to electrons than 
P4 ones, although the density of state for the P4 modes is 
larger. The reason for such an anomalous behavior might 
be due to symmetries of the corresponding phonons as it 
is seen in Fig. [211 Namely, to the P3 peak contribute 
axial vibrations of 0(1) in the Cu2 plane which are odd 
under inversion, while in the P4 peak these modes are 
even. The in-plane modes of Ca and 0(1) are present in 
P3 which are in-phase and out-of-phase modes, while in 
P4 they are all out-of-phase modes. For more informa- 
tion on other modes P5 — P13 see [32|. We stress that 
the Eliashberg equations based on the extracted a'^F{uj) 
of the slightly overdoped Bi2Sr2CaCu20% with the ra- 
tio (2A/Tc) ~ 6.5 describe rather well numerous optical, 
transport and thermodynamic properties [32j . However, 
in underdoped systems with (2A/Tc) ~ 10, where the 
pseudogap phenomena are pronounced, there are seri- 
ous disagreements between experiments and the Migdal- 
Eliashberg theory [32|. We would like to stress that the 
contribution of the high frequency modes (mostly the 
oxygen modes) to a^Fiui) may be underestimated in tun- 
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FIG. 26. The phonon density of states F{ijj) (PDOS) of 
Bi2Sr2CaC'u20s compared with the generalized density of 
states (GPDOS) [l24| . Atomic vibrations: Ol - O in the 
Cu02 plane; 02 - apical O; 03 - O in the BiO plane. From 

El- 

nelling measurements due to their sensitivity to surface 
contamination and defects. Namely, the tunnelling cur- 
rent probes a superconductor to a depth of order of the 
quasi-particle mean-free path liui) — vfJ~^{^)- Since 
7~^(w) decreases with increasing ut the mean free path 
can be rather small and the effects of the high energy 
phonons are sensitive to the surface contamination. 

Similar conclusion regarding the structure of the EPI 
spectral function Q!^P(w) in HTSC cuprates comes out 
from tunnelling measurements on Andreev junctions 
(Z ^ 1 - low barrier) and Giaver junctions {Z ^ 1 - 
high barrier) in La2-xSrxCu04 and YBCO compounds 
[l25j |. where the extracted a^P(aj) is in good accordance 
with the phonon density of states Fphiuj) - see Fig. 

Note, that the BTK parameter Z is related to the 
transmission and reflection coefficients for the normal 
metal {1 + Z'^)~^ and Z'^{1 + Z'^y^ respectively. 

Although most of the peaks in a'^F{uj) of HTSC 
cuprates coincide with the peaks in the phonon den- 
sity of states it is legitimate to put the question - can 
the magnetic resonance in the superconducting state give 
significant contribution to a^P(a;)? In that respect the 
inelastic magnetic neutron scattering measuremen ts of 
the magnetic resonance as a function of doping |l26l | 
give that the resonance energy Er scales with Tc, i.e. 
Er — {5 — 6)Tc as shown in Figi30l This means that 
if one of the peaks in a^P(w) is due to the magnetic 
resonance at lo — Er, then it should shift strongly with 
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FIG. 27. Atomic polarization vectors and their frequencies (in 
meV) at special points in the Brillouin zone. The larger closed 
circles in the lattice are O-ions. F — X is along the Cu-O- 
Cu direction. Arrows indicate displacements. The modes in 
square and round brackets are the transverse and longitudinal 
optical modes respectively, (top) - modes of the PI peak, 
(bottom) - modes of the P2 peak. From [111, [11]. 



FIG. 28. Atomic polarization vectors and their frequencies 
(in meV) at special points in the Brillouin zone. The larger 
closed circles in the lattice are O-ions. F — X is along the Cu- 
O-Cu direction. Arrows indicate displacements. The modes 
in square and round brackets are the transverse and longi- 
tudinal optical modes, respectively. (top) - modes of the P3 
peak, (bottom) - modes of the P4 peak. From [2^, [3^ . 



doping as it is observed in |126l |. This is contrary to 
phonon peaks (energies) whose positions are practically 
doping independent. To this end, recent tunnelling ex- 
periments on Bi-2212 (ssj show clear doping independence 
of a^Fluj) as it is seen in Fig. This remarkable 

result is an additional and strong evidence in favor of 
EPI and against the SFI mechanism of pairing in HTSC 
cuprates which is based on the magnetic resonance peak 
in the s uperconducting state. In that respect the analysis 
in [l27l | of tunneling spectra of electron-doped cuprates 
Pro.88Ceo.i2CM04 with Tc — 24: K shows the existence 
of the bosonic mode at = 16 meV which is signif- 
icantly larger than the magnetic-resonance mode with 
ujr — (10 — 11) meV . This result excludes the magnetic- 



resonance mode as an important factor which modifies 
superconductivity. 

The presence of pronounced phononic structures (and 
the importance of EPI) in the liV) characteristics was 
quite recently demonstrated by the tunnelling measure- 
ments on the very good Lai.85<S'ro.i5Cu04 films prepared 
by the molecular beam epitaxy on the [001]-symmetric 
SrTiO^ bi-crystal substrates fsll- They give unique evi- 
dence for eleven peaks in the (negative) second derivative, 
i.e. —d^I/dV^. Furthermore, these peaks coincides with 
the peaks in the intensities of the phonon Raman scatter- 
ing data measured at 30 K in single crystals of LSCO 
with 20 % of Sr These results are shown in Fig. 

[32l In spite of the lack of a quantitative analysis of the 
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FIG. 29. (a) (fl/dV'^ of a Giaver-like contact in La 
2-xSrxCu04 - note the large structure below 50meV; (b) 
d^I/dV^ of an Andreev- and Giaver-like contact compared to 
the peaks in the phonon density of states. From [l25| |. 



data in the framework of the Eliashberg equations the 
results in Ref. [sij are important evidence that phonons 
are relevant pairing bosons in HTSC cuprates. 

It is interesting th at in the c-axis vacuum tunnelling 
STM measurements (l28j the fine structure in cfl/dV"^ 
at eV > A was not seen below Tc, while the pseudogap 
structure is observed at temperatures near and above Tc- 
This result could mean that the STM tunnelling is likely 
dominated by the nontrivial structure of the tunnelling 
matrix element (along the c-axis) , wh ich is derived from 
the band structure calculations Il29l. However, recent 
STM experiments on Bi — 2212 [S^l give information 
on possible nature of the bosonic mode which couples 
with electrons. In [s^ the local conductance dI/dV{r, E) 
is measured where it is found that d^I/dV^{r,E) has 
peak at E{r) = A(r) -^^^(r) where dI/dV{r,E) has the 
maximal slope - see Fig. ISSTa). 

It turns out that the corresponding average phonon 
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energy fl depends on the oxygen mass, i.e. 51 — i.^^ 
with CtiQ — 52 meV and fiis ~ 48 meV - as it is seen 
in Fig. IMTb). This result is interpreted in [sl] as an ev- 
idence that oxygen phonons are strongly involved in the 
quasi-particle scattering. A possible explanation is put 
forward in ^32i] by assuming that this isotope effect is due 
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FIG. 30. Doping dependence of the energy Er of the magnetic 
resonance peak at tt, tt in YBCO and Bi2212 measur ed at low 
temperatures by inelastic neutron scattering. From [l26| |. 









0«. Oou Oo« Om 



















Up. T/i»ir^ 


IK f\ ll 






00. T\^0 






-< • 


OVO, T,- 


■*X : \ 





40 4a so 
2«V. m»V 



FIG. 31. Second derivative of I{V) for a Bi2212 tunnelling 
junctions for various doping: UD-underdoped; OD-optimally 
doped; OVD-overdoped system. The structure of minima in 
can be compared with the phonon density of states 
F{ui). The full and vertical lines mark the positions of the 
magnetic resonance energy Er ~ 5ATc for various doping 
taken from Fig l30l Red tiny arrows mark positions of the 
magnetic resonance Er in various doped systems. Dotted ver- 
tical lines mark various phonon modes. From 'sS]. 
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phonon mode can be singled out in the spectral function 
a^F(uj) as being the only one which dominates in pair- 
ing mechanism. This important result means that the 
high Tc is not attributable to a particular phonon mode 
in the EPI mechanism but all phonon modes contribute 
to Xep. Having in mind that the phonon spectrum in 
HTSC cuprates is very broad (up to 80 raeV ), then the 
large EPI constant (Aep > 2) obtained in tunnelling ex- 
periments is not surprising at all. Note, that similar con- 
clusion holds for some other oxide superconductors such 
as Bag eKQABiO^ with Tc = 30 K where the peaks in 
the bosonic spectral function /extracted from tunnelling 
measurements) c oinci d e wit h the peaks in the phononic 
density of states jl30l |. [l3l| ). 

E. Phonon spectra and EPI 



FIG. 32. Second derivative data d^I(y)/dV^ of the tun- 
nelling spectra on thin films of Lai.ssSro.isCuOi are shown 
along with phonon Raman scattering data on single crystals 
of LSCO with 20 % Sr. The polarization of the incident and 
scattered light in the Raman spectra is parallel to the Cu02 
planes. From [33 |. 
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FIG. 33. (left) Typical conductance dI/dV{r,E). The ubiq- 
uitous feature at eV > A(gap) with maximal slopes, which 
give peaks in d^I/dV^{r,E) are indicated by arrows, (right- 
b) The histograms of all values of f2(r) for samples with O^^ 
- right curve and with O^* - left curve. From [39l |. 



to the Big phonon which interacts with the anti- nodal 
quasi-particles. However, this result requires a reanaly- 
sis since the energy of the bosonic mode in fact coincides 
with the dip and not with the peak of (Pl/dV'^{r, E) - as 
reported in [39j . 

The important message of most tunnelling experiments 
in HTSC cuprates near and at the optimal doping is 
that there is strong evidence for the importance of EPI 
in the quasi-particle scattering and that no particular 



Although experiments related to phonon spectra and 
their renormalization by EPI, such as inelastic neutron, 
inelastic x-ray and Raman scattering, do not give the 
spectral function a'^F{uj) they nevertheless can give use- 
ful, but indirect, information on the strength of EPI for 
some particular phonons. We stress in advance that the 
interpretation of experimental results in HTSC cuprates 
by the theory of EPI for weakly correlated electrons is 
risky since in strongly correlated systems, such as HTSC 
cuprates, the phonon renormalization due t o EPI is dif- 
ferent than in weakly corre lated metals [132] . Since these 
questions are reviewed in |l32j | we shall briefly enumer- 
ate the main points: (^) in strongly correlated systems 
the EPI coupling for a number of phononic modes can 
be significantly larger than the LDA-DFT and Hartree- 
Fock methods predict. This is du e to m any-body effects 
not contained in LDA-DFT [111], ^M- The lack of the 



LDA-DFT calculations in obtaining phonon line-widths 
is clearly demonstrated, fo r ins tance in experiments on 
L2_2;Sra;Cu04 - see review jl34| and References therein, 
where the bond-stretching phonons at q = (0.3,0,0) 
are softer and much broader than the LDA-DFT cal- 
culations predict. (Note, the wave vector q is in units 
(27r/a, 27r/&, 2tt/c) - for instance in these units (tt, tt) cor- 
responds to (0.5,0.5).) (2) The calculation of phonon 
spectra is in principle very difficult problem since be- 
sides the complexity of structural properties in a given 
material one should take into account appropriately the 
long-range Coulomb interaction of electrons as well as 
strong short-range repulsion. Our intention is not to dis- 
cuss this complexity here - for that see for instance in 
(4^ , but we only stress some important points which will 
help to understand problems with which is confronted 
the theory of phonons in cuprates. 

The phonon Green's function D{q,Lu) depends on 
the phonon self-energy H(q, uj) which takes into ac- 
count all the enumerated properties (note £'^^(q, w) — 
_D(^^(q, w) — H(q, w)). In cases when the EPI coupling 
constant 17 (k, k') is a function on the transfer momentum 
q = k — k' only, then H(g) {q — (q, ia;„)) depends on the 
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quasi-particle charge susceptibility xd^) = P{l)l^e{<i) 



(40) 



and P{q) is the irreducible electronic polarization given 
by 



P{q)^-Y,G{p + q)T,{p,q)G{p). 



(41) 



The screening due to the long-range Coulomb interac- 
tion is contained in the electronic dielectric function 
Seiq) while the" screening" due to (strong) correlations 
is described by the charge vertex function Tdp^q)- Due 
to complexity of the physics of strong correlations the 
phonon dynamics was studied in the t — J mod e l bu t 
with out the long-range Coulomb interaction |l32l |. |l33l |. 
|l35j |. in which case one has Eg = 1 and xdl) = 
However, in studying the phonon spectra in cuprates it 
is expected that this deficiency is partly compensated by 
choosing the bare phonon frequency ljq (q) (contai ned in 
Pn^ {c\,uj)) to correspond to undoped compounds |l32l |. 
[l35j . It is a matter of future investigations to incorpo- 
rate all relevant interactions in order to obtain a fully 
microscopic and reliable theory of phonons in cuprates. 
Additionally, the electron-phonon interaction (with the 
bare coupling constant (/q) is dominated by the change 
of the energy of the Zhang-Rice singlet - see more Part II 
Sect ion V I I, an d in that Eq. pOl) for \l{q) is adequate one 
[I, [ili], [ill]. Since the charge fluctuations in HTSC 
cuprates are strongly suppressed (no doubly occupancy 
of the Cu 3d^ state) due to strong correlations and since 
the suppressed value of Xcit) can not be obtained by the 
band-structure calculations this means that LDA-DFT 
underestimates the EPI coupling constant significantly. 
In the following we discuss this important result briefiy. 




£ 70 
a. 

65 



■DFT (GiusSnoot al.jj 45 



1 0.2 03 4 0.5 
h(r.l.u.) 



-DFT (Bohnen«lal.> 
, ExpQiima;!! 



0.0 0. 1 0.2 3 0.4 0.S 
h(r,l.u.) 



16 
t4 
12 



CIO 



q-(hOO) ^ 



DFT (BotindfiOt aJ.J 
Expetiimni 



0.0 



0.2 0.4 
h(r.l.u.) 



FIG. 34. Comparison of DFT calculations with experimen- 
tal results of inelastic x-ray scattering; (a) phonon energies 
in La\,s5Srii,\r^Cu04, and (b) in Y Ba2CuzOT, (c) phonon 
line widths in Lai.asSroAhCuOA,. DFT calc ulations ilQ] gives 
much smaller width than experiments [l4C| |. From [l39l ]. 



1. Inelastic neutron and x-ray scattering - phonon 
softening and width due to EPI 

The softening and broadening of numerous phonon 
modes has been observed in the normal state of cuprates 
giving evidence for pronounced EPI effects and for inad- 
equacy of the LDA-DFT calculations in treating strong 
correlati ons an d su ppression of the charge susceptibility 
3, 0, [l32j |. Il35i . There are several relevant reviews 



on this subject [6| , [l34j . |l32j |. 136l ] and here we dis- 
cuss briefiy two important examples which demonstrate 
the inefficiency of the LDA-DFT band structure calcu- 
lations to treat quantitatively EPI in HTSC cuprates. 
For instance, the Cu-0 bond-stretching phonon mode 
shows a substantial softening at q^b = (0.3, 0, 0) b y dop - 
ing of Lai.85S'ro.i5CM04 and ySasCMsOr [l33|. pf - 
called the half-breathing ph onon ,, and a large broadening 
by 5 meV at 15 % doping [l40| as it is seen in Fig. IMl 
While the soft ening can be partly described by the LDA- 
DFT method [l37|, the latter theory predicts an order 
of magnitude smaller broadening than the experimental 



one. This failure of LDA-DFT is due to incorrect treat- 
ment of the ffects of strong correlations on the charge 
susceptibility xdl) and due to absence of many body 
effects which can increase the coupling constant - see 
more in Part II. The neutron scattering measurements in 
Lai,s5Sro,i5Cu04 give evidence for large (30 %) soften- 
ing of the of with Ai symmetry with the energy a; ~ 60 
meV, which is theoretically predicted in [l38|, and for 
the large line- width about 17 meV which also suggests 
strong EPI. These apex modes are favorable for d-wave 
pairing since their coupling constants are peaked at small 
momentum q Having in mind the above results, there 
is not surprising at all that the recent calculations of the 
EPI coupling constant Aep in the framework of LDA-DFT 
give very small EPI coupling constant Aep ~ 0.3 [l6|, [17| . 
The critical analysis of the LDA-DFT resuhs in HTSC 
cuprat es is done in ^ and recently additionally argued 
in , |l39| by their disagreement with inelastic neutron 
and x-ray scattering measurements - as it is shown in 
Fig. El 
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In Part II we shall discuss some theoretical approaches 
related to EPI in strongly correlated systems but we 
shall not study the phonon renormalization. The lat- 
ter problem was studied mor e detailed and thoroughly in 
the review article Q, |l34l |. Here, we point out only 
three important results. First, there is an apprecia- 
ble difference in the phonon renormalization in strongly 
and weakly correlated systems. Namely, the change of 
phonon frequencies in the presence of conduction elec- 
trons is proportional to the squared coupling constant 
|(7q| and charge susceptibility Xc) i-G- Suj{q) ~ |(?q| ReXcj 
while the line width is given by ^ ISql^ l-^'^ix^j. All 

these quantities can be calculated in LDA-DFT and as we 
discussed above where for some modes one obtains that 
p(LDA) p(™P) However, it turns out that in strongly 

correlated systems doped by holes (with the concentra- 
tion (5^1) the charge fluctuations are s uppre ssed in 
which case the following sum-rule holds @, [l35j 



ttN 



\ImxA(l,^)\d^^2Sil-6)N, (42) 



theory gives that for coupling to some modes the cou- 
pling constant \gq\ in HTSC cuprates can be significantly 
larger than LDA-DFT predicts Q , which i s due to some 
many-body effects not present in the latter [l33j | . In Part 
II it will be argued that for some phonon modes one has 

2 

ISql^ ^ 9(1^^^ ■ For instance, for the half-breathing 

mode, one has \gq\^ ~ 3 Iffql^^,^ 0, |l33l | - see Part II. 
These two results point to an inadequacy of LDA-DFT in 
calculations of EPI effects in HTSC cuprates. Third, the 
phonon self-energy (H(g)) and quasi-particle self-energy 
S(A:) are differentl y r enormalized by strong correlations 



i, 0, 0, 1271, which is the reason that H(g) is 
much more suppressed than S(fc) - sec Part II. The ef- 
fects of the charge vertex on Il{q) and S(fc) are differently 
manifested. Namely, the vertex function enters quadrat- 
ically in S(A:) and the presence of the forward scattering 
peak in the charge vertex strongly affects the EPI cou- 
pling constant (jq in S(A;) 



- - E l5q • 7c(k, qf D{q)g{k + q), 



(44) 



while in the LDA-DFT method one has 

-^E / |/mx,(q,a.)|^^^d^ = (l-<5)7V. 



(43) 



q#0 



The inequality F 



(LDA) 
cj(q) 



< F 



(cxp) 

'(q) 



(for some phonon 



modes) together with Eqs. (14211^ means that for low 
doping 5 <^ 1 the LDA strongly underestimates the EPI 
coupling constant in the large portion of the Brillouin 

zone, i.e. one has g[^^^^ <C g^^^^ ■ The large soften- 
ing and the large line width of the half-breathing mode 
at q — (0.5,0) but very moderate effects for the breath- 
ing mode at q = (0.5,0.5) is explained in the framework 
of the one slave-boson (SB) theory (for U — oo) in 135|, 
where Xc(qi'^) (i-e. Tc{p,q) = Tc{p,q)) is calculated in 
leading 0{1/N) order. We stress that there is another 
method for studying strong correlations - the X-method, 
where the controllable 1/iV expansion is performed in 
terms of Hubbard operators and where the charge ver- 
tex Fc(p,g) is calculated [3], [EH, It turns out 
that in the adiabatic limit (za;„ = 0) the vertex functions 
rc(PF,q) in two methods have some important differ- 

(X) 

ences. For instance, Fc (pf,?) (in the X-method) is 
peaked at g = - the so called forward scattering peak, 
while F*;"^^-* {pF,q) has maximum at finite g ^ [l4l| - 
see Part II, Section VI. E. The properties of Fc (pf,?) 
are confirmed by numerical Monte Carlo calculations in 
the finitc-U Hubbard model |l68l |. where it is found that 
FSP exists for all U, but it is especially pronounc ed in 
the limit U ^ t. These results are also confirmed in 171 
where the calculations are performed in the four slave- 
bosons technique - see more Part II Section VILE. Having 
in mind this difference it would be useful to have calcula- 
tions of XcIQ' ^) in the framework of the X-method which 
are unfortunately not done yet. Second, the many body 



where g{k){= G{k)/Q) is the quasi-particle Green's func- 
tion, 7^(k, g) = Fc(k, (7)/Q is the quasi-particle ver- 
tex and Q(~ 8) is the Hubbard quasi-particle spectral 
weight - see Part II, Section VII. In the adiabatic limit 
|q| > q^ = ujph/vF one has 7c(k, g) w 7c(k,q) and for 
9 3> qc{~ 6-Tr/a) the charge vertex is strongly suppressed 
(7f,(k, q) ^ 1) making the effective EPI coupling (which 
also enters the pairing potential) small at large (transfer) 
momenta q. This has strong repercussion on the pairing 
due to EPI since for small doping it makes the d-wave 
pairing coupling constant to be of the order of the s-one 
(Arf ~ Xs). Then in the presence of the residual Coulomb 
interaction EPI gives rise to d-wave pairing. On the other 
side the charge vertex Fc(k, g) enters Tl{q) linearly and 
it is additionally integrated over the quasi-particle mo- 
mentum k- see Eq. ([4T|) . Therefore, one expects that 
the effects of the forward scattering peak on H(g) are 
less pronounced than on S(fc). Nevertheless, the peak 
of Fc(k, g) at g = may be (partly) responsible that 
the maximal experimental softening and broadening of 
the stretching (half-breading) mode in L ai .HSi SroA^C'uOi 
and Y Ba2Cu^O-j is at q^j, = (0.3, 0, 0) 134^ and not at 
q/ifc — (0.5, 0) for which g^^ reaches maximum. This 
means that the charge vertex function pushes the renor- 
malized EPI coupling constant to smaller q. It would be 
very interesting to have calculations for other phonons by 
including the vertex function obtained by the X-method 
- see Section. VII 



2. The phonon Raman scattering 

The phonon Raman scattering gives an indir ec t ey 
iden ce for importance of EPI in cuprates Il42l . Il4i 
We enumerate some of them 



144| . We enumerate some of them - see more in [3| and 
References therein, (i) There is a pronounced asymmetric 
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line-shape (of the Fano resonance) in the metalhc state. 
For instance, in Y Ba2C'u30^ two Raman modes at 115 
cm~^ [Ba dominated mode) and at 340 cm^^ [O domi- 
nated mode in the Cu02 planes) show pronounced asym- 
metry which is absent in Y Ba2Cu^0e. This asymmetry 
means that there is an appreciable interaction of Raman 
active phonons with continuum states (quasi-particles). 
(ii) The phonon frequencies for some Aig and Big are 
strongly renormalized in the superconducting state, be- 
tween (6—10) %, poin ting again to the importance of EPI 
[144| - see also in Q , [23f . To this point we mention that 
there is a remarkable correlation between the electronic 
Raman cross-section S'exp(w) and the optical conductiv- 
ity in the a — b plane it(w), i.e. Scxpi^) ^ <^ah{^) 3- 
In previous Subsections it is argued that EPI with the 
very broad spectral function a'^F{uj) (0 < < 80 meV) 
explains in a natural way the uj and T dependence of 
Gabi'^)- This means that the electronic Raman spectra 
in cuprates can be explained by EPI in conjunction with 
strong correlations. This conclusion is su pported by cal- 
culations of the Raman cross-section (l45l] which take into 
account EPI with a^i^(w) extracted from tunnelling mea- 
surements in Y Ba2CuzOQ+x and Bi2Sr2CaCu20g,+x 3, 
[131- nil. Quite similar properties (to cuprates) of the 
electronic Raman scattering, as well as of 17(0;), R{i^) 
and p{T), were observed in experiments [t^ on isotropic 
3D metallic oxides LaQ ^Sro^^CoO^ and Cao.^Sro^^RuOs 
where there are no signs of antiferromagnetic fluctua- 
tions. This means that low-dimensionality and antifer- 
romagnetic spin fluctuations cannot be a prerequisite for 
anomalous scattering of quasi-particles and EPI must be 
inevitably taken into account since it is present in all 
these compounds. 



F. Isotope effect in Tc and E(k,aj) 

The isotope effect axc in th^ critical temperature Tc 
was one of the very important proof for the EPI pair- 
ing in low-temperature superconductors (LTSC). As a 
curiosity the isotope effect in LTSC systems was mea- 
sured almost exclusively in mono-atomic systems and in 
few polyatomic systems; the hydrogen isotope effect in 
PdH , the Mo and 5*6 isotope shift of T^ in MogS'eg, and 
the isotope effect in Nb^Sn and MgB2- We point out 
that very small (ckTc ~ in Zr and Ru) and even nega- 
tive (in PdH) isotope effect in some polyatomic systems 
of LTSC materials are compatible with the EPI pairing 
mechanism but in the presence of substantial Coulomb 
interaction or lattice anharmonicity. The isotope effect 
aTc cannot be considered as the smoking gun effect since 
it is sensitive to numerous influences. For instance, in 
MgB2 it is with certainty proved that the pairing is due 
to EPI and strongly dominated by the boron vibrations, 
but the boron isotope effect is significantly reduced, i.e. 
axc ~ 0.3 and the origin for this smaller value is still un- 
explained. The situation in HTSC cuprates is much more 
complicated because they are strongly correlated systems 



and contain many- atoms in unit cell. Additionally, the 
situation is complicated with the presence of intrinsic 
and extrinsic inhomogeneities, low dimensionality which 
can mask the isotope effects. On the other hand new 
techniques such as ARPES, STM, ^iSR allow studies of 
the isotope effects in quasi-particle self-energies, i.e. as, 
which will be discussed below. 



L Isotope effect Tc 

This problem will be discussed only briefly since more 
extensive discussion can be found in Q. It is well 
known that in the pure EPI pairing mechanism, the to- 



tal isotope coefficient a is given by = X^i p ^ 



(p) 



-J2z,pd^nTc/d\nM^''\ where M^''' is the mass of the 
i-th element in the p-th crystallographic position. We 
stress that the total isotope effect is not measured in 
HTSC cuprates but only some partial ones. Note, 
that in the case when the screened Coulomb interac- 
tion is negligible, i.e. /i* = 0, the theory predicts 
axc = 1/2. From this formula one can deduce that 
the relative change of Tc, 5Tc/Tc, for heavier elements 
should be rather small - for instance it is 0.02 for 
i35Ba ^138 0.03 for ^^Cu (j^ ^^^^^ g.O? for 

^^^La This means that measurements of 

ai for heavier elements are at/or beyond the ability of 
present day experimental techniques. Therefore most 
isotope effect measurements were done by substituting 
light atoms ^^O by ^^O only. It turns out that in 
most optimally doped HTSC cuprates ao is rather small. 
For instance ao ~ 0.02 — 0.05 in Y Ba2Cu^Oj with 
Tc.max ~ 91 K, but it is appreciable in Lai,g,zSrQ,i^CuOn 
with Tcmax ~ 35 X where ao ~ 0.1 — 0.2. In 
Bi2Sr2CaCu20s with Tc_niax « 76 iiT one has ao ~ 
0.03 — 0.05 while ao ~ 0.03 and even negative (—0.013) 
in Bi2Sr2Ca2Cu20io with Tc^max ~ 110 if. The ex- 
periments on Tl2Can-iBaCun02n+i {n — 2,3) with 
Tc.max ~ 121 K are still unreliable and ao is unknown. 
In the electron-doped [N di-xCex)2CuOi with Tc_max ~ 
24 K one has ao < 0.05 while in the underdoped materi- 
als ao increases. The largest ao is obtained even in the 
optimally doped compounds like in systems with substi- 
tution, such as Lai.sBSro i^Cui-xMxOi, M = Fe, Co, 
where ao ~ 1-3 for x « 0.4 %. In La2-xMxCuO^ there 
is a Cu isotope effect which is of the order of the oxy- 
gen one, i.e. acu ~ cao giving acu + cto ~ 0.25 — 0.35 
for optimally doped systems {x = 0.15). In case when 
X = 0.125 with Tc < C max one has acu ~ 0.8 — 1 
with acu + oto « 1.8 [l46j |. The appreciation of copper 
isotope effect in La2-xMxCuOi tells us that vibrations 
other than oxygen ions are important in giving high Tc. 
In that sense one should have in mind the tunnelling ex- 
periments discussed above, which tell us that all phonons 
contribute to the Eliashberg pairing function Q!^i^(k,w) 
and according to these results the oxygen modes give 
moderate contribution to Tc [32]. Hence the small oxy- 



r(p) 
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gen isotope effect in optimally doped cuprates, if it 
is an intrinsic property at all (due to pronounced local 
inhomogeneities of samples and quasi-two dimensionality 
of the system), does not exclude the EPI mechanism of 
pairing. 

2. Isotope effect qe in the self-energy 

The fine structure of the quasi-particle self-energy 
S(k,CLi), such as kinks and slopes, can be resolved in 
ARPES measurements and in some respect in STM. It 
turns out that there is isotope effect in the self -e nerg y 
in th e optimally doped Bi — 2212 samples |105L |108j . 
[lOO]. In the first paper on this subject [l05j | it is re- 
ported a red shift Sujk.ra —(10— 15) meV of the nodal 
kink at ujk.io — 70 meV for the ^^O — ^ O^^ substitution. 
In (1051] it is reported that the isotope shift of the self- 
energy (5S = Ei6 — Si8 ~ 10 meV is very pronounced at 
large energies oj = 100 — 300 meV. Concerning the latter 
result, there is a d i spute since it is not confirmed exper- 
imentally in [lOSl j. [l09j . How e ver, t he isotope effect in 
i?eS(k, w) at low energies |l08j |. |l09j | is well described in 
the framework of the Migdal-Eliashberg theory for EPI 
[l07| which is in accordance with the recent A RPE S mea- 
surements with low-energy photons ~ 7 eV [l47l |. The 
latter allowed very good precision in measuring the iso- 
tope effect in the nodal point of Bi-2212 with T^^ = 92.1 
K and T^^ = 91.1 K They observed shift in the 

maximum of i?el](kjv, w) - at uik^ro ~ 70 meV (it corre- 
sponds to the half-breathing or to the breathing phonon) 
by Sujk^Yo ~ 3.4 ± 0.5 meV as shown in Fig. 1351 

By analyzing the shift in Iml](kjv,i^) - shown in Fig. 
[551 one finds similar result for Sujk^jo ~ 3.2 ± 0.6 meV. 
The similar shift was obtained in STM measurements [s^ 
which is shown in Fig. l33r b) and can have its origin in dif- 
ferent phonons. We would like to stress two points: (?) in 
compounds with Tc ^ 100 X the oxygen isotope effect in 
Tc is moderate, i.e. ai^J < 0.1, 147]. If we consider this 
value to be intrinsic then even in this case it is not in con- 
flict with the tunnelling experiments [s^] since the latter 
give evidence that vibrations of heavier ions contribute 
significantly to Tc - see the discussion in Subsection D 
on the tu nnel ling spectroscopy; (ii) in ARPES measure- 
ments of jl47l ] the effective EPI coupling constant Xep,eff 
> 0.6 is extracted, while the theory in Subsection C gives 
that the real coupling constant is larger, i.e. Xep > 1.2. 
This value is significantly larger than the LDA-DFT the- 
ory predicts Xep.LDA < 0.3 [3], [13] ■ This again points 
that the LDA-DFT method does not pick up the many- 
body effects due to strong correlations- see Part II. 

IV. SUMMARY OF PART I 

The analysis of experimental data in HTSC cuprates 
which are related to optics, tunnelling and ARPES mea- 
surements near and at the optimal doping give evidence 
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FIG. 35. (a) Effective ReT, for five samples for O^^ (blue) and 
O^* (red) along the nodal direction, (b) Effective JmE deter- 
mined from MDC full widths. An impurity term is subtracted 
at a; = 0. From ^147.] . 

for the large electron-phonon interaction (EPI) with the 
coupling constant 1 < X^p < 3.5. We stress that 
this analysis is done in the framework of the Migdal- 
Eliashberg theory for EPI which is a reliable approach for 
systems near the optimal doping. The spectral function 
a^Fluj), averaged over the Fermi surface, is extracted 
from various tunnelling measurements on bulk materials 
and tin films. It contains peaks at the same energies as 
the phonon density of states Fph{u>). So obtained spectral 
function when inserted in the Eliashberg equations pro- 
vides sufficient strength for obtaining high critical tem- 
perature Tc ~ 100 K. These facts are a solid proof for 
the important role of EPI in the normal state scattering 
and pairing mechanism of cuprates. Such a large (exper- 
imental) value of the EPI coupling constant and the ro- 
bustness of the d-wave superconductivity in the presence 
of impurities imply that the EPI and impurity scatter- 
ing amplitude must be strongly momentum dependent. 
The IR optical reflectivity data provide additional sup- 
port for the importance of EPI since by using the spec- 
tral function (extracted from tunnelling measurements) 
one can quantitatively explain frequency dependence of 
the dynamical conductivity, optical relaxation rate and 
optical mass. These findings related to EPI are addi- 
tionally supported by ARPES measurements on BISCO 
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compounds. The ARPES kinks, the phononic features 
and isotope effect in the quasi-particle self-energy in the 
nodal and anti- nodal points at low energies (w <C Wc) per- 
sist in the normal and superconducting state. They are 
much more in favor of EPI than for the spin fluctuation 
(SFI) scattering mechanism. The transport EPI coupling 
constant in HTSC cuprates is much smaller than Agp, 
i.e. Xtr ~ Aep/3, which points to some peculiar scatter- 
ing mechanism not met in low-temperature superconduc- 
tors. The different renormalization of the quasi-particle 
and transport self-energies by the Coulomb interaction 
(strong correlations) hints to the importance of the small- 
momentum scattering in EPI. This will be discussed in 
Part II. 

The ineffectiveness of SFI to solely provide pairing 
mechanism in cuprates comes out also from magnetic 
neutron scattering on YBCO, BISCO. As a result, the 
imaginary part of the susceptibility is drastically reduced 
in the low energy region by going from slightly under- 
doped toward optimally doped systems, while Tf. is prac- 
tically unchanged. This implies that the real SFI cou- 
pling constant As/(~ g^j:) is small since the experimental 

value g^™^'' < 0.2 eV is much smaller than the assumed 

theoretical value g^^j^^ « (0.7 — 1.5) eV. 

Inelastic neutron and x-ray scattering measurements in 
cuprates show that the broadening of some phonon lines 
is by an order of magnitude larger than the LDA-DFA 
methods predict. Since the phonon line- widths depend 
on the EPI coupling and the charge susceptibility it is ev- 
ident that calculations of both quantities are beyond the 
range of applicability of LDA-DFT. As a consequence, 
LDA-DFT overestimates electronic screening and thus 
underestimates the EPI coupling, since many-body ef- 
fects due to strong correlations are not contained in this 
mean- field type theory. However, in spite of the promis- 
ing and encouraging experimental results about the dom- 
inance of EPI in cuprates the theory is still confronted 
with difficulties in explaining sufficiently large coupling 
constant in the d-channel. At present there is no such a 
satisfactory microscopic theory although some concepts, 
such as the the dominant EPI scattering at small trans- 
fer momenta, are understood at least qualitatively. These 
set of problems and questions will be discussed in Part 
IL 

Part II 

Theory of EPI in HTSC 

The experimental results in Part I give evidence that 
the electron-phonon interaction (EPI) in HTSC cuprates 
is strong and in order to be conform with d-wave pair- 
ing EPI must be peaked at small transfer momenta. A 
number of other experiments in HTSC cuprates give ev- 
idence that these are strongly correlated systems with 
large on-site Coulomb repulsion of electrons on the Cu- 



ions. However, at present there is no satisfactory mi- 
croscopic theory of pairing in HTSC cuprates which is 
capable to calculate Tc and the order parameter. This 
is due to mathematical difficulties in obtaining a solu- 
tion of the formally exact ab initio many-body equations 
which take into account two important ingredients - EPI 
and strong correlations Q. In Section V we discuss 
first the ab initio many body theory of superconduc- 
tivity in order to point places which are most difficult 
to be solved. Since the superconductivity is low energy 
phenomenon (also in HTSC cuprates) one can simplify 
the structure of the ab initio equations in the low-energy 
sector (the Migdal-Eliashberg theory) where the high- 
energy processes are incorporated in the (so called) ideal 
band structure (non-local) potential V7b5(x, y) and ver- 
tex functions T. Due to its complexity this program of 
calculation of V/_bs(x, y), F, EPI couplings gep{x,y) is 
not realized in HTSC superconductors. However, one 
pragmatical way out is to calculate g^p in the framework 
of the LDA-DFT method which is at present stage un- 
able to treat strong correlations in a satisfactory manner. 
Some achievements and results of LDA-DFT which are 
related to HTSC cuprates are discussed in Section VI. 

In the case of very complicated systems, such as the 
HTSC cuprates, the standard (pragmatical) procedure 
in physics is to formulate a minimal theoretical model 
- sometimes called toy model, which includes minimal 
set of important ingredients necessary for qualitative and 
and semi-quantitative study of a phenomenon. As the 
consequence of the experimental results, the minimal the- 
oretical model for cuprates must comprise two important 
ingredients: 1. EPI and 2. strong correlations. In Sec- 
tion VII we shall introduce such a minimal theoretical 
t — J model which includes EPI and discuss the renor- 
malization of EPI by strong correlations. In recent years 
the interest in this set of problems is increased and nu- 
merous numerical calculations were done mostly on small 
clusters with n x n atoms (n < 8). We shall not discuss 
this subject which is fortunately covered in the recent 
comprehensive review [6|. The analytical approaches in 
studying the renormalization of EPI by strong correla- 
tions, which are based on a controllable and systematic 
theory, are rather scarce. We shall discuss such a system- 
atic and controllable theory in the framework of the t— J 
model with EPI, which is formulated and solved in terms 
of Hubbard operators. The theory of this (toy) model 
predicts some interesting effects which might be impor- 
tant for understanding the physics of HTSC cuprates. 
It predicts that the high-energy processes (due to the 
suppression of doubly occupancy for U 3> Wb) give rise 
to a nonlocal contribution to the band structure poten- 
tial (self-energy I](x, y,a; = 0)) as well as to EPI. This 
non-locality in EPI is responsible for the peak in the ef- 
fective pairing potential (Vev.effiQ, ^)) at small trans- 
fer momenta qc <^ kp) 0, HH, [o^. The latter 
property allows that the (strong) EPI is conform with 
d-wave pairing in HTSC cuprates. Furthermore, the pe- 
culiar structural properties of HTSC cuprates and cor- 
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responding electronic quasi-two-dimensionality give an 
additional non-locality in EPI. The latter is due to the 
weakly screened Madelung energy for a number of lat- 
tice vibrations along the c-axis. Since at present there is 
no quantitative theory for the latter effect we tackle this 
problem here only briefly. The next task for the future 
studies of the physics of HTSC cuprates is to incorpo- 
rate these structural properties in the minimal theoreti- 
cal model problem. 

Finally, by writing this chapter our intention is not to 
overview the theoretical studies of EPI in HTSC cuprates 
- an impossible task, but more (i) to elucidate the de- 
scending way from the (old) well-defined ab initio micro- 
scopic theory of superconductivity to the one of the min- 
imal model which treats the interplay of EPI and strong 
correlations, (ii) to encourage the reader to further de- 
velop the theory of HTSC cuprates. 



V. MICROSCOPIC THEORY OF 
SUPERCONDUCTIVITY 

A. Ab iniitio many-body theory 

The many-body theory of superconductivity is based 
on the fully microscopic electron-ion Hamiltonia n fo r 
electrons and ions in the crystal - see for instance 148|, 
[l49j |. It comprises mutually interacting electrons which 
also interact with ions and ionic vibrations. In or- 



der to pass continually to the problem of interplay of 
EPI and strong correlations and also to explain why 
the LDA-DFT method is inadequate for HTSC cuprates 
we discuss this problem here with restricted details - 
more extended discussion can be found in i, [Hi. In 
order to describe superconductivity the Nambu-spinor 

4'\r) = (V'|(r) ■0|(r)) is introduced which operates in 

the electron- hole space (V'(r) = ( V'^r))^') where '0-|~(r), 

'0|(r) are annihilation and creation operators for spin 
up, respectively etc. The microscopic Hamiltonian which 
in principle describes normal and superconducting states 
of the system under consideration contains three parts 
H = He -\- Hi + He-i- The electronic Hamiltonian He, 
which describes the kinetic energy and the Coulomb in- 
teractions of electrons is given by 



He 



d'^rip (r)f3eo(p)-0(r) 



+ ^ I d^r(fr'i;\r)T3i;{r)Ve{r~r')^j\r')T3^j{r'), (45) 



1 
2 

where eo(p) = fP' /2ni is the kinetic energy of electron 
and Vc{y — y') = e? j \ r — r' | is the electron-electron 
Coulomb interaction. Note, in the electron-hole space 
the pseudo-spin (Nambu) matrices t^, i = 0,1,2,3 are 
Pauli matrices. Since we shall discuss only the electronic 
properties the explicit form of the lattice Hamiltonian Hi 



3; [l49j | is omitted here. The electron-ion Hamiltonian 
describes the interaction of electrons with the equilibrium 
lattice and with its vibrations, respectively 



He-i 



d^rVe- 



Kr-R°)V'\r)f3^(r) 



d-^r^{r)^ (r)f3V'(r) 



(46) 



Here, T4_i(r — R°) is the electron-ion potential where Ze 
is the ionic charge and in the all-electron calculations it 
is given by V;_,(r - RO) = - ZeV | r - R° |. The sec- 
ond term which is proportional to the lattice distortion 

operator $(r) = - Y.„,a Uan"^ aVe-i{r - R°) + ^anh{r) 

(because of convenience it includes also the EPI cou- 
pling VctT4_i) describes the interaction of electrons with 
harmonic (^ Uan) (or anharmonic ~ $an/i(r)) lattice vi- 
brations. 

The Dyson's equations for the electron and phonon 
Green's functions G_(l,2) = -(r^/'(l)V^^(2)), D{1 - 2) ^ 
-(T$(l)$(2)) are G-\l,2) = G(7^(l,2) - E(l,2) and 
l)-i(l,2) = Do^il,2) - n(l,2), where the Gq^{1,2) = 
[(-g77 - eo(Pi) + IJ-)to - Ueff{l)T3]5{l - 2) is the bare 
inverse electronic Green's function. Here, 1 = (ri,Ti), 
where ri is the imaginary time in the Matsubara tech- 
nique and the effective one-body potential Me//(1) = 
Ve-i{l) + Vh + +($(1)), Vh is the Hartree potential. 
The electron and phonon self-energies E(l, 2) and 11(1, 2) 
take into account many-body dynamics of the interacting 
system. The electronic self-energy E(l,2) = Sc(l,2) + 
Sep(l, 2) is obtained in the form 

t{l, 2) = -Veff{l, i)f3G(l, 2)fe//(2, 2; I), (47) 

where integration (summation) over the bar indices is 
understood. The effective retarded potential Vef / (1, 1) in 
Eq. (l47|) contains the screened (by the electron dielectric 
function ee(l,2)) Coulomb and EPI interactions 



Vefsil, 2) = Ve{\-\)e-e^{\, 2)+e-i(l, 1)^(1, ~2)e-e\% 2). 

(48) 

The inverse electronic dielectric permeability is 



Hl,2) 



,5(1 - 2) + K(l - l)P(l,2)e-i(2,2) is 



defined via the irreducible electronic polarization op- 
erator P(l,2) = -5p{f3G(l, 2)f e//(2, 3; 2)G(3, 1+)}. 
The vertex junction re//(l,2;3) = — (5G(1, 2)/(5ue//(3) 
in Eq. (jT7|) is the solution of the complicated (and 
practically unsolvable) integro-differential func- 
tional equation fe//(l,2;3) = f3(5(l - 2)5(1 - 3) 
{5S(1, 2)/<5G(i, 2)}G(1, 3)f e//(3, 4; 3)G(4, 2). Note, 
that the effective vertex function re//(l;2;3), which 
takes into account all renormalizations going beyond the 
simple Coulomb (RPA) screening, is the functional of 
both the electronic and phononic Green's functions G 
and D, thus making at present the ab initio microscopic 
equations practically unsolvable. 
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B. Low-energy Migdal-Eliashberg theory 

If the vertex function Tef / would be known we would 
have closed set of equations for Green's functions which 
describe dynamics of interacting electrons and lattice vi- 
brations (phonons) in the normal and superconducting 
state. However, this is a formidable task and at present 
far from any practical realization. Fortunately, we are 
mostly interested in low-energy phenomena (with en- 
ergies \ujn \ ,^<^ujc and for momenta k = kp + Sk in 
the shall Sk Skc near the Fermi momentum kp] Wc 
and Skc are some cutoffs which disappear after calcu- 
lations) which allows us further simplification of equa- 
tions 1]. Therefore, the strate gy is to integrate high- 
energy processes - see more in [l49j, and we sketch it 
only briefly. Namely, the Green's function (^(k, a;„) — 
[iujn — (k^/2TO — /i)-r3 — S(k, w„)]~^ can be formally writ- 
ten in the form G(k,w„) = G"°'"(k,w„) -f G'''3^{k,ujn), 
where G''°"'(k,w„) = G'(k,w„)e(wc - |a;„|)e((5/cc - Sk) 
is the low-energy Green's function and G''*^''(k, cj„) = 
G{'k,u}n)'d{\oJn\—^c)'d{Sk — Skc) the high-energy one and 
analogously D = D'°'"(k, w„) + i:)'**f'*(k, w„). By intro- 
ducing the small parameter of the theory s ^ (w/wc) ^ 
(Sk/Skc) 1 one has in leading order G'°'"(k,a;„) ~ s"^, 
G''*9''(k,a;„) < 1 and Z)'°^(k,w„) - s", D^*9''(k, a;„) - 
s^. Note, that the coupling constants {Vei,Wei, Vu etc.) 
are of the order s° = 1. 

This procedure of separating low-energy and high- 
energy processes lies also behind the adiabatic approxi- 
mation since in most materials the characteristic phonon 
(Debye) energy of lattice vibrations lojj is much smaller 
than the characteristic electronic Fermi energy Ep 
{lod < Ep). In the small s(< 1) hmit the Migdal the- 
ory Q keeps in the total self-energy I] linear terms in the 
phonon propagator D [D) only. In that case the effective 
vertex function can be written in the form Fe/ / (1, 2; 3) = 
fc(l, 2; 3) (Jf ep(l, 2; 3) Q, where the Coulomb charge 

vertex fc(l, 2; 3) = f 3<5(1-2)5(1-3)-H5Ec(1, 2)/Sueff{i) 
contains correlations due to the Coulomb interaction only 
but does not contain EPI and phonon propagator D 
explicitly. The part (5fep(l,2;3) = Step{l,2) / Suef f{2>) 
contains all linear terms with respect to EPI. Note, 
that in these diagrams enters the dressed Green's func- 
tion which contains implicitly EPI up to infinite order. 
By careful inspection of all (explicit) contributions to 
STep{^^ 2; 3) which is linear in D one can express the self- 
energy in terms of the charge (Coulomb) vertex rc(l, 2; 3) 
only. As a result of this approximation, the part of the 
self-energy due to Coulomb interaction is given by 

Sc(l, 2) = ~V^%1, l)f3G(l, 2)fe(2, 2; T), (49) 

where V/=(l,2) = Vc(l, 2)£-i(2, 2) is the screened 
Coulomb interaction. The part which is due to EPI has 
the following form 

Sep(l, 2) = -Kp(l, 2)f,(l, 3; 1)G(3, 4)f e(4, 2; 2), (50) 



where 14^(1,2) = £^1(1, l)l)(l, 2)£-i(2, 2) is the 
screened EPI potential. Note, that I]ep(l,2) depends 
now quadratically on the charge vertex Fc, which is due 
to the adiabatic theorem. 

It is well known that the Coulomb self-energy Sc(l, 2) 
is the most complicating part of the electronic dynamics 
but since we are interested in low-energy physics when 
s ^ 1 then the term Ec(l, 2) can be further simplified by 
separating it in two parts 

S,(1,2)-E(")(1,2) + eW(1,2). (51) 

The term Sc''^(l,2) is due to high-energy processes con- 
tained in the product G"f'»(i^ 2)f ^*f'»(2, 2; 1) (for in- 
stance due to the large Hubbard U in strongly corre- 
lated systems) and Sc''(l,2) is due to low energy pro- 
cesses. The leading part part of Ec''^(l,2) is 1, i.e. 
tf\l,2) - fiO, while t^\l,2) is small of order 1, i.e. 
Sc ''(1, 2) ^ s^. For further purposes we define the quan- 
tity Vq 

Vb(l,2) = {K-.,(1) + T^h(1)}t3<5(1-2) + S('')(1,2). (52) 

where Ve-i, Vh are also of order s°. After the Fourier 
transform with respect to time (and for small |aj„| <C uJq) 
si'''' is given by 

S(")(X1,X2,C^„) ~ Eio^(xi,X2,0)f3 + (Sio^)'(xi,X2,0)-i^, 

(53) 

As we said E^q"* ~ while (e|;q'')' ■ ujn ^ s^ because 
From Eq.dMl) it is seen that the part t^c\l, 2) 
contains the low-energy Green's function G'°'"(l,2) and 
this skeleton diagram is of order . The similar analysis 
based on Eg. (15(11) for Eep(l, 2) gives that the leading order 
is which describes the low-energy part of EPI. After 
the separations of terms (of and orders) the Dyson 
equation in the low- energy region has the form 

[ia;„Zc(x,x) - Ho{^,^) - ^c\^,^,(^n) - Sep(x, x,cj„)] 

xG''""(x,y,w„) = '5(x-y)ro, (54) 

where x means integration J d^x over the crystal vol- 
ume. The Coulomb renormalization function Zc(x, y) — 

(5(x — y) — (Sq|!^)'(x, y, 0) and the single-particle Hamil- 
tonian Hq (x, y) collect formally all high-energy processes 
which are unaffected by superconductivity (which is low- 
energy process) where 

7?o(x,y) = {(-^V^ - ^)J(x - y) + Vo^")(x,y, 0)}f3 

(55) 

with 

Fj'')(x,y,0) = {K^,(x) + T/H(x)}(5(x-y)-Hl]i^)(x,y,0). 

(56) 
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One can further absorb Zc(x, y) into the renormalized 
Green's function 

Gr{x,y,iu,,) = zy2(x,x)G"°-(x,y,,.„)zy2(y,y), 

(57) 

the renormalized vertex function rren(l,2;3) = 
Zc ^^^TcZc ^^"^j the renormalized self-energies 

and introduce the ideal band- structure Hamiltonian 
ho{x,y) = Ztr^/^(x,x)i7o(i,y)Zc"^^^(y,y) given by 

/io(x, y) = {(-^V^ - ^l)6{^ - y) + Vibs{^: y)}T3. 

(58) 

Here, 

ViBsi^.y) - Z-^'\^,^)vt\x,y)Z-''\y,y) (59) 

is the ideal hand structure potential (sometimes called the 
excitation potential) and apparently nonlocal quantity, 
which is contrary to the standard local potential V^(x) 
in the LDA-DFT theories - see Section VI. The static po- 
tential V7b5(x, y) is of order s° and includes high-energy 
processes 

Finally, we obtain the matrix Dyson equation for the 
renormalized Green's function Gr(x,y,aj„) which is the 
basis for the (strong-coupling) Migdal-Eliashberg theory 
in the low-energy region 

[iw„(5(x - x) - /io(x, x) - Si'|(x, x,w„) - Sep,r(x, x,w„)] 

X G^(x,y,tj„) =(5(x-y)fo, (60) 

where Ec'r and Sep,r have the same form as Eqs.(|49t- 
1501) but with renormalized Green's and vertex functions 
Gr, fr instead of G, f. We stress that Eq.® holds 
in the low-energy region only. In the superconducting 
state the set of Eliashberg equations in Eq. (pO)) are 
writen explicitely in Appendix A, where it is seen that 
the superconducting properties depend on the Eliashberg 
spectral function a\^F{ijj). The latter function is defined 
also in Appendix A, Ea. (jl04p . and it depends on material 
properties of the system. 

The important ingredients of the low-energy Migdal- 
Eliashberg theory is the ideal band-structure Hamilto- 
nian /io(x, y) - given by Ea.(|58p. which contains many- 
body (excitation) ideal band-structure non-local periodic 
crystal potential yfss(x, y). The Hamiltonian /io(x, y) 
determines the ideal energy spectrum e(k) of the conduc- 
tion electrons and the wave function p (x) through 

^o(x, y)-0j,k(y) h(k) - Ail'i/'i.k W, (61) 

where /i is the chemical potential. We stress that the 
Hamiltonian ho{x,y) also governs transport properties 
of metals in low-energy region. 

After solving Ea. ([6T|) the next step is to expand all 
renormalized Green's function, self- energies, vertices and 



the renormalized EPI constant (symbolically gep.r = 
gepfl^ren/se) in the basis of "i/jj p(x) and write down the 
Eliashberg equations in this basis. We shall not elaborate 
further this program and we refer the reader to the rele- 
vant literature [ 148i] , |l49j |. We point out, that even such 
simplified program of the low-energy Migdal-Eliashberg 
theory was never fully realized in low-temperature super- 
conductors, because the non-local potential ViBs{^,y) 
(enters the ideal band-structure Hamiltonian ho{x,y)) 
and the renormalized vertex function (entering the EPI 
coupling constant gep,r) which include electronic corre- 
lations are difficult to calculate especially in strongly 
correlated metals. Therefore, it is not surprising at all, 
that the situation is even worse in HTS materials which 
are strongly correlated systems with complex structural 
and material properties. Due to these difficulties the cal- 
culations of the electronic band structure and the EPI 
coupling is usually done in the framework of LDA-DFT 
where the many-body excitation potential ViBs{^,y) is 
replaced by some (usually local) potential Vlda{^) which 
in fact determines ground state properties of the crystal. 
In the next section we briefly describe: (z) the LDA-DFT 
procedure in calculating the EPI coupling constant and 
(a) some results related to HTSC cuprates. We shall also 
discuss why this approximation is a serious drawback if 
applied to HTS materials. 

VI. LDA-DFT CALCULATIONS OF THE EPI 
COUPLING 

We point out again two results which are important 
for the future microscopic theory of pairing in HTSC 
cuprates. First, numerous experiment (discussed in Part 
I) give evidence that the EPI coupling constant which 
enters the normal part of the quasi-particle self-energy 
X^p = Xs + Xd + •■• is rather large, i.e. 1 < < 3.5. In 
order to be conform with d-wave pairing the effective EPI 
potential must be rather nonlocal (peaked at small trans- 
fer momenta q) which implies that the s-wave and d-wave 
coupling constants are of the same order, i.e. ~ As. 
Second, the theory based on the minimal t — J model, 
which will be discussed in Section VII, gives that strong 
electronic correlations produce a peak at small transfer 
momenta in the effective EPI pairing potential thus giv- 
ing rise to Xd ~ A^. This is a striking property which 
allows that EPI is conform with d-wave pairing. How- 
ever, the theory is seriously confronted with the problem 
of calculation of the coupling constants A^. It turns out, 

that at present it is an illusory task to calculate A^ and 
Xd since it is extremely difficult (if possible at all) to 
incorporate the peculiar structural properties of HTSC 
cuprates (layered structure, ionic-metallic system, etc.) 
and strong correlations effects in a consistent and reli- 
able microscopic theory. As it is stressed several times, 
the LDA-DFT methods miss some important many body 
effects (especially in the band structure potential) and 
therefore fail to describe correctly screening properties 
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of HTSC cuprates and the strength of EPI. However, the 
LDA-DFT methods are capable to incorporate structural 
properties much better than the simplified minimal t — J 
(toy) model. Here, we discuss briefly some achievements 
of the improved LDA-DFT calculations which are able to 
take partially into account some nonlocal effects in the 
EPI. The latter are mainly due to the almost ionic struc- 
ture along the c-axis which is reflected in the very small 
c-axis plasma frequency ujc ^ ^ab)- 

The main task of the LDA-DFT theory in obtaining 
the EPI coupling is to calculate the change of the ground 
state (self-consistent) potential SVg{r) /SRa and the EPI 
coupling constant {matrix element) g^^^Ck^k') (see its 
definition below) what is the most difficult part of calcu- 
lations. Since in LDA-DFT there is no EPI coupling, the 
recipe is that the calculated "EPI" coupling is inserted 
into the Eliashberg equations. By knowing (7^^"^(k, k') 
one can define the total (A) and partial (Aq.^) i^P/ cou- 
pling constants for the i^-th mode, respectively |l50l | 



A 



A. 



PI. 



7rA^(0)Wq,^' 



(62) 



where p = 3k is the number of phonon branches (k is the 
number of atoms in the unit cell) and N{0) is the density 
of states at the Fermi energy (per spin and unit cell). 
The phonon line-width 7^ ^ is defined in the Migdal- 
Eliashberg theory by 
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e"(q)-5a,i/'(k,k-q) 



(63) 

Here, e"(q) is the phonon polarization vectors, is the 
Fermi function. Since the ideal energy spectrum from 
Eq. (j6ip k — -^',k ~ M s-nd the corresponding eigen- 
functions ?/'ki unknown then instead of these one sets 
in Eq. ^ the LDA - DFT eigenvalues for the I - th 

band Cl.t^'^^ and V'Lt^^^- In the LDA-DFT method the 

EPI coupling constant {matrix element) g^^r^^ is de- 
fined which by the change of the ground-state potential 
6Vg{r)/SR^ 
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E^l^^r). (64) 



The index n means summation over the lattice sites, 
a = x,y^ z and the wave function ^A^t^"^^ are the so- 
lutions of the Kohn-Sham equation - see Q. In the 
past various approximations within the LDA — DFT 
method have been used in calculating SVg{r)/SRa and 
A while here we comment some of them only: (i) In 
most calculations in LTS systems and in HTSC cuprates 
the rigid-ion {RI) approximation was used as well as 
its further simplifications which inevitable (due to its 



shortcomings and obtained small A) deserves to be com- 
mented. The RI approximation is based on the very 
specific assumption that the ground-state (crystal) po- 
tential Vg{r) can be considered as a sum of ionic po- 
tentials Vg{r) = '^n^gi^ ~ where the ion poten- 
tial Vg{r — R„) and the electron density Pe(r) are carried 
rigidly with the ion at R„ during the ion displacement 



(R„ 



1). In the RI approximation the change 



of Vg (r) is given by 



SVg{r)^J2^^Vg{r-R'i^)u. 



SVg{r) 



V„F3(r-R° 



(65) 
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which means that RI does not take into account changes 
of the electron density during the ion displacements. In 
numerous calculations applied to HTSC cuprates the 
rigid-ion model is even further simplified by using the 
rigid muffin-tin approximation {RAIT A) (or similar ver- 
sion with t he ri gid- atomic sphere) - see discussions in 
[l50t . [mi, [113. The RMTA assumes that the ground- 
state potential and the electron density follow ion dis- 
placements rigidly inside the Wigner-Seitz cell but out- 
side it Vg (r) is not changed because of the assumed very 
good metallic screening (for instance in simple metals) 



V„l/g(r-R„) = { 



S/aVg{r — Rn), r in cell n 
0, outside 



(67) 



This means that the dominant EPI scattering is due to 
the nearby atoms only and that the scattering potential 
is isotropic. All nonlocal effects related to the interac- 
tion of electrons with ions far away, are neglected in the 
RMTA. In this case 5^_^'^(k, k') is calculated by the 
wave function centered at the given ion R° which can 
be expanded inside the muffin-tin sphere (outside it the 
potential is assumed to be constant) in the angular mo- 
mentum basis {l,m}, i.e. 
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a„(fc,RO)ft(|r-R° !)!,„((/), 0) 

(68) 

(the angles cj), 9 are related to the vector f = (r — R°)/ | 
r — R° I). The radial function p;(| r — R° |) is zero out- 
side the muffin-tin sphere. In that case the EPI matrix 
element is given by g^_^^'^'^(k, k') - (Yi™ | f \ Yv^') 
and because f is a vector the selection rule implies that 
only terms with = I' — I = ±1 contribute to the 
EPI coupling constant in the RMTA. This result is 
an immediate consequence of the assumed locality of 
the EPI potential in RMTA. However, since nonlo- 
cal effects^ such as the long-range Madelung-like interac- 
tion, are important in HTSC cuprates then additional 
terms contribute also to the coupling constant gam i-e. 
g„,„(k, k') = g«r^^(k, k') + sS^n '°'(k, k'), where a part 
{^9a°n''°'^) of the nonlocal contribution to 52°"'°'^ is rep- 
resented schematically 



Sg. 



nonloc I 



k,k')'-(15™|(R^-R^)a|i1'. 



(69) 
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From Ea. ([69|) comes out the selection rule Al = Ir — I = 
for the nonlocal part of the E — P interaction. We 
stress that the Al — (nonlocal) terms are omitted 
in the RMTA approach and therefore it is not surpris- 
ing that this approximation works satisfactory in ele- 
mental (isotropic) metals only. The latter are charac- 
terized by the large density of states at the Fermi sur- 
face which makes electronic screening very efficient. This 
gives rise to a local EPI where an electron feels poten- 
tial changes of the nearby atom only. One can claim 
with certainty, that the RMTA method is not suitable 
for HTSC cuprates which are highly anisotropic sys- 
tems with pronounced ionic character of binding and 
pronounced strong electronic correlations. The RMTA 
method applied to HTSC cuprates misses just this im- 
portant part - the long-range part EPI due to the change 
of the long range Madelung energy in the almost ionic 
structure of HTSC cuprates. For instance, the first cal- 
culations done in 1153] which are based on the RMTA 

/? K IT' A 

give very small EPI coupling constant A ^ 0.1 in 

YBCO, which is in apparent contradiction with the ex- 
perimental finding that Aep is large - discussed in Part 
I. 

Howe ver, these nonlocal effects are taken into account 
in [l50j | by using the frozen-in phonon (FIP) method 
in evaluating of Xep in La2-xMxCu04. In this method 
some symmetric phonons are considered and the band 
structure is calculated for the system with the super-cell 
which is determined by the periodicity of the phonon 
displacement. By comparing the unperturbed and per- 
turbed energies the corresponding EPI coupling X^, (for 
the considered phonon z/-th mode) is found. More pre- 
cisely speaking, in this approach the matrix elements of 
6Vg{r)/6RQ ^ are determined from the finite difference of 
the ground-state potential 

AK,,q,.(r) = VgiRl^ + Ar« ,(i)) - VgiR^^) 



where L, k enumerate elementary lattice cells 
and atoms in the unit cell, respectively. The 
frozen-in atomic displacements of the phonon 
Arqj,(L) of the v-th mode is given by At'^^{L) — 

Auq,^(;i/2M^Wq,i,)i/2i?e[e^_^(q)e"i-^] where Auq^^ is 
the dimensionless phonon amplitude and the phonon 
polarization (eigen) vector eK,iy(q) fulfills the condition 
^^e* j^(q) • eK.i/'(q) = S,y^^>. Based on this approach 
various symmetric Ag (and some B^g) modes of 
La2-xMxCuOi were studied [l50| . where it was found 
that the large matrix elements are due to unusually 
long-range Madelung-like, especial for the c-axis phonon 
modes. The obtained large \ep ~ 1.37 is the consequence 
of the following three main facts: (i) The electronic 
spectrum in HTSC cuprates is highly anisotropic, i.e. 
it is quasi-two-dimensional. This is an important fact 



for pairing because if the conduction electrons would be 

uniformly spread over the whole unit cell then due to the 

rather low electron density [n ^ Kfi^cm"^) the density 

of states on the Cu and O in-plane atoms would be an 

order of magnitude smaller than the real value. This 

would further give an order of magnitude smaller EPI 

coupling constant Agp. Note, that the calculated density 

of states on the (heavy) Cu and (light) O in-plane 

atoms, iV^"(0) ^ 0.54 states/ eV and AfO-«(0) 0.35 

states /eV^ are of same order of magnitude as in some 

LTS materials. For instance, in NhC where Tc « 11 

K one has on (the heavy) Nb atom N^^{fi) ^ 0.58 

states/eV and on (the light) C atom N^\lS) ^ 0.25 

states /eV. So, the quasi-two-dimensional character of 

the spectrum is crucial in obtaining appreciable density 

of states on the light O atoms in the Cu02 planes, (ii) 

In HTSC cuprates there is strong Cu — O hybridization 

leading to good in-plane metallic properties. This large 

covalency in the plane is due to the (fortunately) small 

energy separation of the electron levels on Cu and 

0^,1 atoms which comes out from the band structure 
I h 

calculations [154| . i.e. A =| ecu — ^o^y \^ 3 eV . The 
latter value gives rise to strong covalent mixing (the 
hybridization parameter tpd) of the Cufi^2_y2 s-nd Op^ ^ 
states, i.e. tpd = —1.85 eV. It is interesting, that 
the small value of A is not due to the ionic structure 
(crystal field effect) of the system but it is mainly due 
to the natural falling of the Cud^2_y2 states across the 
transition-metal series. So, the natural closeness of the 
energy levels of the Cud^2^y2 ^^id Op^^ states is this 
distinctive feature of HTSC cuprates which basically 
allows achievement of high Tc- {Hi) The ionic structure 
of HTSC cuprates which is very pronounced along the 
c-axis is responsible for the weak electronic screening 
along this axis and according to that for the significant 
contribution of the nonlocal (long-range) Madelung-like 
interaction to EPI. It turns out that because of the 
ionicity of the structure the La and Oz axial modes are 
strongly coupled with charge carriers in the Cu02 planes 
despite the fact that the local density of states on these 
atoms is very small |150( . i.e. ^^^"(0) = 0.13 states /eV 
and A^*^^ (0) = 0.016 states /eV. For comparison, on 
planar atoms Cu and Oxy one has A^'^"(0) = 0.54 
states /eV and N'^-y{0) = 0.35 states /eV. These 
calculations show that the lanthanum mode (with 
cjq iy — 202 cm^^) at the q = (0, 0.27r/c) zone boundary 
(fully symmetric Z-point) has ten times larger coupling 
constant \^1{FIP) — 4.8 than it is predicted in the 
RMT approximation X^l{RMT) = 0.48. The similar 
increase holds for the average coupling constant, where 

>ilverage{FlP) = 1-0 but X^%^^^^g^{RMT) = 0.1. 

Note, that for the q « La-mode one obtains 
A^"(F/P) = 4.54 compared to X^^iRMT) ^ 0.12. Sim- 
ilar results hold for the axial apex-oxygen q — (0, 0.27r/c) 
mode (Oz) with Wq.^ = 396 cm,~^ where the large 
(compared to the RMT method) coupling constant 
is obtained: Aq^, — 14 and A^J^^^g^^gg = 4.7, while 
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for 



cm 



axial apex-oxygen modes with Wq.^ = 415 



one has A, 



= 11.7. 



After 



'u,average " — averaging 

over all calculated modes it was estimated A — 1.37 
and Ld\og w 400 K. By assuming that /i* = 0.1 one 
obtains T^. — 4:9 K by using Allen-Dynes formula for 
Tc w 0.83a;iogexp{-1.04(l -I- A)/[A - fi*{l + 0.62A)]} 
with wiog — 2 J du!duja'^{uj)F{uj)hiuj/ Xuj. For /i* — 0.15 
and 0.2 one obtains Tc — 41 and 32 K, respectively. 
We stress that the rather large Aep (and Tc) are due to 
the nonlocal (long range) effects of the metallic-ionic 
structure of HTSC cuprates and no n- mu ffin-tin correc- 
tions in EPI, as was first proposed in [l55l |. However, we 
would like to st ress that the optimistic results for Aep 
obtained in (l50j | is in fact based on the calculation of the 
EPI coupling for some wave vectors q with symmetric 
vibration patterns and in fact the obtained Agp is an 
extrapolated value. The all-q calculations of Aep,q which 
take into account long-range effects is a real challenge 
for the LDA-DFT calculations and are still awaiting. 

Finally, it is worth to mention important calculations 
of the EPI coupling constant in the framework of the 
linear-response full-potential linear- muffi n-tin- orbital 
method (LRFP — LMTO) invented in |l5a | and applied 
to the doped HTSC cuprate (Cai -^,Sr ^)\ jiCuO-2 for 
X ~ 0.7 and y - 0.1 with = 110 K |157t . Namely, 



these calculations give strong evidence that structural 
properties of HTSC cuprates already alone make the 
dominance of small-q scattering in EPI, which effect 
is additionally increased by strong corr elations. In 
order to analyze this compound in |l57l | calculations 
are performed for CaCu02 doped by holes in an 
uniform, neutralizing back-ground charge. The momen- 
tum (q — (q||,gj_)) dependent EPI coupling constant 
(summed over all phonon branches v) in different L chan- 
nels (s,p, d..) is calculated by using a standard expression 



Ai(q) - Af^ri(k + q) |5k,q,.|'>L(k) x 5{£,^+^)5{i^) 

(71) 

Here, is the quasi-particle energy, ffk.q.z^ is the EPI 
coupling constant (matrix element) with the v-th branch, 
Yl (k) is the L-channel wave function and the normaliza- 
tion factor M cx A''^"'^(0) with the partial density of states 
is Nl{0) oc X^k (^)'^(^k)- The total coupling constant 
in the L-channel is an average of AL(q) over the whole 
2D Brillouin zone (over qy), i.e. XL{q±) = (A_L(qii )) py 
We stress three important results of Ref. jl57j . First, 
the s- and d-coupling constants, As(q), Ad(q), are peaked 
at small transfer mom enta q ^(7r/3a, 0, 0) as it is shown 
in Fig. 3 of Ref. jl57j . This result is mainly caused by 
the nestin g pro perties of the Fermi surface shown in Fig. 
1 of Ref. |157| . Second, the q-dependence of the inte- 
grated EPI matrix elements |gLq| =AL(q)/xi(q) (with 
xKq) « Eki^i(k + q)>^L(k)<5(ek+q)<5(Ck)) for L = s,d 
is similar to that of Ai(q), i.e. these are peaked at 
small transfer momenta q <C 2kp. Both these results 
mean that the structural properties of HTSC cuprates 



imply the dominance of small-q EPI scattering. Third, 
the calculations give similar values for Xs{q± — 0) and 
^d{q± = 0), i.e. Xs = 0.4 7 for s — wave and A^ = 0.36 
for d — wave pairing (l57| . This finding (that A^ ~ As) 
is due to the dominance of the small g-scattering in EPI. 
It means that the nonlocal effects (long range forces) in 
EPI of HTSC cuprates are very important. This result 
together with the finding of the dominance of small-g 
scattering in EPI due to strong correlations [5l|, (96| . 
(OTj mean that strong correlations and structural prop- 
erties of HTSC cuprates make EPI conform with d-wave 
pairing, either as its main cause or as its supporter. We 
stress that the obtained coupling constant Ad = 0.36 is 
rather small to give d-wave pairing with large Tc and on 
the first glance this result is against the EPI mechanism 
of pairing in cuprates. However, it is argued throughout 
this paper, that the LDA methods applied to strongly 
correlated systems overestimate the screening effects and 
underestimate the coupling constant and therefore their 
quantitative predictions are not reliable. 



VII. EPI AND STRONG CORRELATIONS IN 
CUPRATES 

A. Minimal model Hamiltonian 

The minimal microscopic model for normal and su- 
perconducting state of HTSC cuprates must include at 
least three orbitals: one (i2,2_j,2-orbital of the Cu-ion and 
two p-orbitals {px,y) of the 0-ion since they participate 
in transport properties of these materials - see more in 
3] and References therein. The electronic part of the 
Hamiltonian (of the minimal mode l) is H = Hq + Hint 
- usually called the Emery model [158|, where the one 
particle tight-binding Hamiltonian Hq describes kinetic 
energy in the three-band (orbital) model 



E 



ifj',c.pPjaaP3'f}<y (72) 



Here t^^^ enumerate Cu- and 0-sites, respectively) 
is the hopping integral between the Paia = x, y)- and d- 
states and t^j'ap between pa- and p/3-states, while e° and 
epc, are the bare d- and p- local energy levels and fi is 
the chemical potential. This tight-binding Hamiltonian 
is written in the electronic notation where the charge- 
transfer energy A^p o = ^3 ~ ^p > by assuming that 
there is one Sd^^-y^ electron on the copper {Cu"^^) while 
electrons in the p-levels of the ions occupy filled 

bands. Hq contains the main ingredients coming from 
the comparison with the LDA — DFT band structure 
calculations. The LDA— DFT results are reproduced by 
assuming t^P <^ tpd (and e^^, — ^p) where the good fit to 
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the LDA — DFT band structure is found for A^p^o = <^d^ 
6° « 3.2 eV &nAV"^{= tpd) = (V3/2)(pdcr), {pda) = -1.8 
eV. The total LDA band- width Wb ^ (4v^) | |= 9 
eV [T59t . 

The electron interaction is described by iJi„( 



where Ud and ?7p are the on-site Coulomb repulsion en- 
ergies at Cu and O sites, respectively while Vc and 
describe the long-range part of the Coulomb interaction 
of electrons (holes) and EPI, respectively. Note, the 
Hubbard repulsion Ud on the Cu-ion is different from its 
bare atomic value Udo{~ 16 eV for C u) due to various 
kind of screening effects in solids |l6Clj | . In turns out that 
in most transition metal oxides one has Ud <SC Udo- This 
problem is thoroughly studied in (l60j | and applied to 
HTSC cuprates. An estimation from numerical cluster 
calculations [161] gives Ud = 9 — 11 eV and C/p = 4 — 6 
eV but because nfUd ^ n^Up the on-site repulsion on 
oxygen is usually neglected at the first stage of the anal- 
ysis. 

Note, that in the case of large Ud{^ tpd, ^dp,o) the hole 
notation is usually used where in the parent compound 
(and for \tpd\ ^ ^dp,o) one has (nf) — 1, i.e. one hole in 
the 3(i-shall (in the ^d^2_y2 state) in the ground state. 
In the limit of large Ud oo the doubly occupancy on 
the Cu atoms is forbidden and only two copper states are 
possible: Cm^"*" described by the quantum state dj^ | 0) 
and Cu^+ described by | 0). In this (hole) notation the 
oxygen p-level is fully occupied by electrons, i.e. there is 
no holes ( (np = 0) in the occupied oxygen 2p-shall of 
O^^ . In this notation the vacuum state | 0„) (not the 
ground state) of the Hamiltonian H for large Ud corre- 
sponds to the closed shell configuration Cu^~^0^~. In 
the hole notation the hole p-level e^^ lies higher than the 
hole d- level e°^, i.e. Apdfl = e^^ — e",^ > (note in the 
electron picture it is opposite) and Ud means repulsion of 
two holes (in the 3dx2-y2 orbital ) with opposite spins - 
3d® configuration of the Cu'^^ ion. Note, that e^i^ = — Ep, 

^dh ~ ~^d ^^^'^ tpd^h — —tpd- In the following the index 
h in tpd,h is omitted. The reason for 6°,^ > ejij^ is partly 
in different energies for the hole sitting on the oxygen 
and copper, respectively [l59j |. From this model one can 
derive in the limit U — » cxj t h e t — J model for the 2D lat- 
tice in the Cu02 plane jl62j | , jl63| , where each lattice site 
corresponds to a Cw-atom. In the presence of one hole in 
the 3(i-shall then in the undoped (no oxygen holes) HTSC 
cuprate each lattice site is occupied by one hole. By dop- 
ing the systems with holes the additional holes go onto 
O-sites. Furthermore, due to the strong Cu-0 covalent 
binding the energetics of the system implies th at an O- 
hole forms a Zhang- Rice singlet with a Cu-hole (l62j |. In 
the t—J model the Zhang-Rice singlet is described by an 
empty site. Since in the t— J model the doubly occupancy 
is forbidden one introduces annihilation (Hubbard) oper- 



describes creation of an hole (in the 3d-shall of the Cu- 
atoms) on the j-th site if this site is previously empty 
(thus excluding doubly occupancy) , i.e. ni^„ -{-Ui^^cr < 1 
must be fulfilled on each lattice site. In this picture the 
doped hole concentration 5 means at the same time the 
concentration of the oxygen holes and of the Zhang- Rice 
singlets. 

The bosonic-like operators X'^^'^^ — X^^'^ X'j''^'^ , ai ^ 
02 create a spin fluctuation for (Ti ^ 02 at the i — th site 
and the spin operator is given by S = X^^°{(T)gj^a2^i'^'^ 
where summation over bar indices is understood. The 
operator X^"^ has the meaning of the hole number 
on the i-th site. It is useful to introduce the operator 
= X^'^X['^ at a given lattice site which is the number 
of Zhang-Rice singlets on the i-th site, i.e. if Xf^ \ 0) = 
1 I 0) the i-th site is occupied by the Zhang- Rice singlet, 
while for X?" | 1) = | 1) there is no Zhang-Rice singlet 
on the i-th site (i.e. this site is occupied only by one 
3d^ hole on the Cu site). The latter is due to the local 
constraint 



X\ 
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xr = 1, 



(74) 



which forbids doubly occupancy of the i-th site by holes. 
By projecting out doubly occupied (high energy) states 
the t—J model reads 



+ ^ Jij{St ■ Sj - - fiifij ) + H3. 



(75) 



The first term (^ e°) describes an effective local energy 
of the hole (or the Zhang- Rice singlet), the second one 
(^ tij) describes hopping of the holes, the third one 
(^ Jij) is the Heisenber g-lik e exchange energy between 



two holes. The theory |162l | predicts that len ^ \t 



This property is very important in the study of EPI. H3 
contains three-sites term which is usually omitted believ- 
ing it is not important. For charge fluctuation processes 
it is plausible to omit it, while for spin- fluctuation pro- 
cesses it is questionable approximation. If one introduces 
the enumeration a,/3,7,A = 0,t,i than the Hubbard op- 
erators satisfy the following algebra 



7/3 



(76) 



where Sij is the Kronecker symbol. The 
(anti) commutation relations in Eq. (j76l) are more 
complicated than the canonical Fermi and Bose 
(anti) commutation relations which complicates the 
mathematical structure of the theory. Note, that the 
Hubbard operators possess also projection properties 



with xr^x7^ 



yXf"^. To escape these complications 



ator of the composite fermion X 



aO 



4o-(l-rtj,-o-) which 



some novel techniques have been used, such as the 
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slave boson technique. In this technique X?'^ = ficrb], 
^(Tio-a _ jj^^fia^ are represented in terms of the fermion 
(spinon) operator fiu which annihilates the spin of the 
hole and the boson (holon) operator b\ which creates 
the Zhang-Rise singlet. 

In the minimal theoretical model the electron-phonon 
interaction (EPI) contains in principle two leading terms 
F^p = Hl^+H1°\ the "ionic" one and the "cova- 

lent" one (7?^°^). The "ionic" term describes the change 
of the energy of the hole (or the Zhang-Rice singlet) at 
the i-th site due to lattice vibrations and it reads j3| , [5l| , 

M 



(77) 



where the "displacement" operator — '^L^lf-i^^i — 
R-Lac + - - £(1^° - R-D] (which as in Section V 
includes the bare coupling constant) describes the change 
of the hole (or Zhang-Rice singlet) energy e" ^ by dis- 
placing atoms in the lattice by the vector ulk- In the 
harmonic approximation the EPI potential is given by 
= I].9j(q; A)exp{iqRJ[6q,A + ^-q^aI where 6q,A and 
6q are annihilation and creation operators of phonons 
with the polarization A, respectively. This term describes 
in principle the following processes: (1) the change of the 
O-hole and Cw-hole bare energies e"^, in the three- 
band model due to lattice vibrations; (2) the change of 
the long-range Madelung energy (which is due to the ion- 
icity of the structure) by lattice vibrations along the c- 
axis; (3) the change of the Cu — O hopping parameter tpd 
in the presence of vibrations, etc. Here, L and k enumer- 
ate unit lattice vectors and atoms in the unit cell, respec- 
tively. Until now the phonon operator is calculated 
in the harmonic approximation for the EPI interaction 
of holes with some specific phononic modes, such as the 
breathing and half- breathing ones Q, |l33l |. The theory 
which includes also all other (than oxygen) vibrations in 
$i is still awaiting. 

It is interesting to make comparison of the EPI cou- 
pling constants in the t—J model and in the Hartree-Fock 
(HF) approximation (which is the analogous of the LDA- 
DFT method) of the three-band model in Eas. (|72j|73p 
when the problem is projected on the single band. For 
instance, the coupling constant with the half-breathing 
mode at the zone boundary in HF approximation (which 
mimics LDA-DFT) is given by 
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while the coupling constant in the t-J model gjii,'^i= 
de°/dRcu~o) is given by 
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where p ~ 0.96 - see Q, |l33l | and references therein. It 
is obvious that in the t — J model the electron-phonon 
coupling is different from the HF one, since the former 
contains an additional term coming from the many body 
effects, which are not comprised by the HF (LDA-DFT) 
calculations. The first term in Eq. ([79|) describes the 
hopping of a 3d-ho\e into the O 2p-states and this term 
exists also in the LDA-DFT coupling constant - see Eq. 
(l78l) . However, the second term in Eq. ((79)) . which is due 
to many body effects, describes the hopping of a O 2p- 
hole into the (already) single occupied Cu 3d-state and it 
does not exist in the LDA-DFT approach. Since the cor- 
responding dimensionless coupling constant A^f, is pro- 
portional to \ghb\ one obtains that the bare t-J coupling 
constant is almost three times larger than the LDA-DFT 
one A^^"' ~ 3A^{,^. This example demonstrates clearly 
that LDA-DFT method is unreliable for calculating the 
EPI coupling constant in HTSC cuprates. 

Note, that there is also a covalent contribution to EPI 
which comes from the change of the effective hopping (i) 
in the t—J model Eg. (175)) and the exchange energy (J) 
in the presence of atomic displacements 
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Here, we shall not go into details but only stress that 
since |e^| 3> \tij \ then the covalent term in the effective t- 
J model is m uch smaller than the ionic term - see more in 
[3 , Si , (133| | and References therein, and in the following 
only the term H'fJ" will be considered [|, [HI, [9^. 



B. Controllable X-method for the quasi-particle 
dynamics 

The minimal model Hamiltonian for strongly corre- 
lated holes with EPI (discussed above) is expressed via 
the Hubbard operators which obey " ugly" non-canonical 
commutation relations. The latter property is rather 
unpleasant for making a controllable theory in terms 
of Feynmann diagrams (for these " ugly" operators) and 
some other approaches are required. There is one very 
efficient and mathematically controllable approach for 
treating the problem with Hubbard operators without 
using slave boson (or fermion) techniques. The method 
(we call it the X-method) is based on the general Baym- 
Kadanoff technique where the expansion for the 

Green's functions in te rms of Hubbard operators was 
first introduced in [l64j and refined in [5l[, [hI]. In the 
paramagnetic and homogeneous state (with finite dop- 
ing) the Green's function Ga^a^{l — 2) is diagonal, i.e. 
C,,(l - 2) = 5a,aM^ - 2) where 



G(l 



(rl"-(l)l-"(2)) =g(l-2)g, (81) 
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with the Hubbard spectral weight Q = (X"") + {X° 

The function ^(1 — 2) plays the role of the quasi-particle 
Green's function - see more in [5l|, [1^, [97.] . It turns 
out, that in order to have a controllable theory (1/iV 
expansion) one way is to increase the number of spin 
components from two to N by changing the constraint 
Eg. ([71)1 into the new one 
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(82) 



In order to reach convergence of physical quantities in the 
limit iV — )• oo the hopping and exchange energy are also 
re-scaled, i.e. Uj = to^ij/N and Jij = Jo,ij/N. In order 
to eliminate a possible misunderstanding we stress one 
important fact, that in the case N > 2 the constraint in 
Ea.(|82|) spoils some projection properties of the Hubbard 
operators. Fortunately, these (lost) projection properties 
are not used at all in the refined theory. As the result one 
obtains the functional integral equation for G{1,2), thus 
allowing unambiguous mathematical and physical treat- 
ment of the problem. In , foB] , [23] it is developed a 
systematic expansion for the quasi-particle Green's 
function 5(1 - 2)(= go + gi/N + ...), Q(== Nqo + qi + ...) 
(also for G(l — 2)) and the self-energy. For large A^(— > oo) 
the leading term is Go(l - 2) = 0{N) = go{l - 2)Qo 

with 50 = 0(1) and Qa = = N6/2. Here, S 

is the concentration of the oxygen holes (i.e. of the 
Zhang-Rice singlets) which is related to the chemical 
potential by the equation 1 — 6^2 np (p) with 

np{p) — (^e'^«(^)~f^ + l) ^. The quasi-particle Green's 
function go(k, i^) and quasi-particle spectrum eo(k) in the 
leading order are given by 



go(k,w) = 



uj - [eo(k) - fj] ' 



eo(k) = ec - (5 • t(k} - ^ Jo(k + p)nF(p). 



(83) 



(84) 



The level shift is Ec = e*^ + 2 J2p i{p)^F{p) and t{p) is the 
Fourier transform of the hopping integral tij - see more 
in i3j. 

Let us summarize the main results of the X-method 
in leading 0(l)-order for the quasi-particle properties in 
the t-J model (1, [HI, [ol], i9I]: (i) The Green's function 
(/o(k, w) describes the coherent motion of quasi-particles 
whose contribution to the total spectral weight of the 
Green's function Go(k, w) is Qq = N5/2. The coherent 
motion of quasi-particles is described in leading order by 
Go(k, w) = (3o5o(k, w) and the quasi-particle residuum 
Qo disappears in the undoped Mott insulating state {S = 
0). This result is physically plausible since in the Mott 
insulating state the coherent motion of quasi-particles, 
which is responsible for finite conductivity, vanishes; (ii) 
The quasi-particle spectrum eo(k) plays the same role as 



the eigenvalues of the ideal band structure Hamiltonian 
/io(x, y) (it contains the excitation potential Vf_Bs(x, y) 
which includes high-energy processes due to the Coulomb 
interaction). So, if we would consider etf,(k) = ^^(k) 
as the tight-binding parametrization of the LDA-DFT 
band-structure spectrum which takes int account only 
weak correlations (with the local potential Vxc{x)S{x — 
y)) then one can define a non-local excitation potential 

^iBsi.^^ y) = ^iBsi^^ y) + K;c(x)(5(x - y) which mimics 
strong correlations in the t-J model 
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(x,y)«K,<5(x-y) + (l-(5)t(x-y) 



J(x-y). 

(85) 

Here, Vq = 2 ^p t(p)ni?(p) and i(x — y) is the Fourier 

transform of t(k) while J(x — y) is the Fourier trans- 
form of the third term in Eq. ((84|) . The relative exci- 
tation potential Vj^g (x, y) is due to strong correlations 
(suppression of doubly occupancy on each lattice site) 
and as we shall see below it is responsible for the screen- 
ing of EPI in such a way that the forward scattering 
peak appears in the effective EPI interaction - see dis- 
cussion below, (iii) For very low doping eo(k) is dom- 
inated by the exchange parameter if Jq > StQ. How- 
ever, in the case when Jq <^ Sto there is a band nar- 
rowing by lowering the hole doping 5, where the band 
width is proportional to the hole concentration S, i.e. 
Wb = z-5-to; (iv) The 0(l)-order quasi-particle Green's 
function 5o(k, w) and the quasi-particle spectrum eo(k) 
in the X-method have similar form as the spinon Green's 
function go,/(k, w) — — (T/cr/J)k,w and the spinon en- 
ergy 65 (k) in the 5*5- method. However, in the SB 
method there is a broken gauge symmetry in metallic 
state which is characterized by (6^) ^ 0. This broken 
local gauge symmetry in the slave-boson method in 0(1) 
order, which is due to the local decoupling of spinon and 
holon is in fact forbidden by to the Elitzur's theorem. On 
the other side the local gauge invariance is not broken in 
the X-method where the Green's function Go(k, w) de- 
scribes motion of the composite object, i.e. simultaneous 
creation of the hole and annihilation of the spin at a given 
lattice site, while in the SB theory there is a spin-charge 
separation because of the broken symmetry {{bi) 7^ 0). 
The assumption of the broken symmetry (bi) 7^ gives 
qualitative satisfactory results for quasi-particle energy 
for the case = 00 in D > 2 dimensions. However, 
the analysis of response functions of the system and of 
higher order 1/A^ corrections to the self-energy in the 
SB theory is very delicate and special techniques must 
be implemented in order to restore the gauge invariance 
of the theory. On the other side the X-method is in- 
trinsically gauge invariant and free of spurious effects in 
all orders of the 1/A^ expansion. Therefore, one expects 
that these two methods may deliver different results in 
0(1) order in response function. This difference is al- 
ready manifested in the calculation of EPI where charge 
vertices are peaked at different wave vectors q, i.e. at 
q in the X-method and g 7^ in the SB method - 
see Subsection E; (v) In the important paper |l65| it is 
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shown that in the superconducting state the anomalous 
self-energies (which are of 0(l/A^)-order in the ex- 
pansion) of the X- and 5-B-methods diffe r substantially. 
As a consequence, the S'B-method [166] predicts false 
superconductivity in the t — J model (for J = 0) with 
large Tc (due to the kincmatical interacti on), w hile the 
X-method gives extremely small Tc{^ 0) '165*1. So, al- 
though the two approaches yield some similar results in 
leading 0(l)-order they are different at least in the next 
leading 0(l/A^)-order. 



C. EPI effective potential in the t-J model 

The theory of EPI in the minimal t — J model based 
on the X-method predicts that leading term in the EPI 
self-energy Egp is given by the expression [5l|, [11] 



Sep(l, 2) = -Kp(l - 2)7,(1, 3; l)5o(3 - 4)7,(4, 2; 2), 

(86) 

where the screened (by the dielectric constant) EPI po- 
tential Vep{l - 2) = £ji(l - I)Kp (1 - 2)£r^(2 - 2) and 
^^0(1 - 2) = -(r4(l)d(2)) is the "phonon" propaga- 
tor which may also describe an anharmonic EPI. It is 
obvious that Ea. (|5B|) is equivalent to Eg. ([5(11) in spite 
the fact that the theory is formulated in terms of the 
Hubbard operators. The charge vertex 7,(1, 2; 3) — 
— SgQ^ (1,2)/ Sueff (3) corresponds to the the renormal- 
ized vertex Tc^r m Eg. dSH)) and it describes the "screen- 
ing" by strong correlations. It depends on the relative 
excitation potential Vjgg{x,y). The electronic dielec- 
tric function £,(1 — 2) describes screening of EPI by the 
long-range part of the Coulomb interaction. Note, that 
in the harmonic approximation $(1) contains the bare 
EPI coupling constant and lattice displacement u, 

i.e. $ ~ ffep'"' " see more in (Note, that in the above 
equations summation and integration over bar indices are 
understood.) The self-energy Eep(k, w) due to EPI reads 
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withi?(w,fl) = -2Tri{nB{n) + l/2) + il;{l/2 + i)-il;{l/2- 
i{Vl+LL:) /2ttT) where nsi^) is the Bose distribution func- 
tion and '0 is di-gamma function. The Eliashberg spectral 
function is given by 



a2^^(k,k',c.)=iV(0)^ 



5i,(k,k- k') 



X (5(L^-w,(k-k'))7c(k,k-k') 



where Qi, (k, p) is the EPI coupling constant for the v-t\ve 
mode, where the renormalization by long-range Coulomb 
interaction is included, i.e. (/i/(k, p) = ^,(k, p)/ee(p)■ 
(...)k denotes Fermi-surface average with respect to the 
momentum k and A^(0) is the density of states renormal- 
ized by strong correlations. The effect of strong correla- 
tions (in the adiabatic limit) is stipulated in the charge 



vertex function 7,(k, k — k') which, as we shall see be- 
low, changes the properties of Vep(q, ^2) drastically com- 
pared to weakly correlated systems. In fact the charge 
vertex depends on frequency but in the adiabatic limit 
[ujph ^ W) and for qvp > i^ph it is practically frequency 
independent, i.e. 7i'"''' (k, q, cj) « 7,(k, q, w ~ 0) where 
the latter is real quantity. For J = the 1/N expan- 
sion gives A^(0) = A^o(0)/qo where qq = S/2 in the t — t' 
model. For J 7^ the density of states N{0) does not di- 
verge for J ^- where 7V(0)(- 1/Jo) > A^o(O). The bare 
density of states Nq(0) is calculated in absence of strong 
correlations, for instance by the LDA -~ DFT scheme. 

Depending on the symmetry of the superconducting or- 
der parameter A(k, u) (s— , d—wave pairing) various pro- 
jected averages (over the Fermi surface) of (k, k',a;) 
enter the Eliashberg equations. Assuming that the su- 
perconducting order parameter transforms according to 
the representation F^ of the point group Civ of the square 
lattice (in the Cu02 planes) the appropriate symmetry- 
projected spectral function is given by 
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k and k' are momenta on the Fermi line in the irre- 
ducible Brillouin zone (1/8 of the total Brillouin zone). 
Tj , = 1, ..8, denotes the eight point-group transforma- 
tions forming the symmetry group of the square lattice. 
This group has five irreducible representations which we 
distinguish by the label i = 1,2,. ..5. In the following 
we discuss the representations i — 1 and « = 3, which 
correspond to the s— and d — wave symmetry of the full 
rotation group, respectively. Di{j) is the representation 
matrix of the j — th transformation for the representation 
i. Assuming that the superconducting order parameter 
A(k, w) does not vary much in the irreducible Brillouin 
zone one can average over k and k' in the Brillouin zone. 
For each symmetry one obtains the corresponding pairing 
spectral function o^Fiiio) 



a'F,iuj)^{{a'F,ik,k,uj))r)r, 



(90) 



which governs the transition temperature for the order 
parameter with the symmetry F^. For instance a^F^^uj) 
is the pairing spectral function in the d-channel and it 
gives the coupling for d—wave superconductivity (the ir- 
reducible representation F3 - sometimes labelled as Big). 
Performing similar calculations for the phonon-limited 
resistivity one finds that the resistivity is related to the 
transport spectral function a'^Ftr{tu) 

a^.F(.)=<<-^^(^''-:'-)[-Wr^^'^^^'^^"'- (91) 



2((v2(k))) 
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The effect of strong correlations on EPI was discussed 
in [96| within the model where (/i/(k, p) and the phonon 
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frequencies (k — k ) are weakly momentum dependent. 
In order to elucidate the main effect of strong correlations 
on EPI and a^Fi(w) we consider the latter functions for 
a simple model with Einstein phonon, where these func- 
tions are proportional to the (so called) relative coupling 
constant A,- 
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^((|7c(k,k-r,k')P»kk'A(j) (92) 



Similarly, the resistivity p(T){'^ Xtr ^ ^tr) is renormal- 
ized by the correlation effects where the transport cou- 
pling constant Atr is given by 
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(93) 

As we see, all projected spectral functions afF{uj) de- 
pend on the charge vertex function 7^ (k, q) which de- 
scribes the screening (renormalization) of EPI due to 
strong correlations (suppression of doubly occupancy) 
[Slj . |96| . This important ingredient (which respects also 
Ward identities) is decisive step beyond the MFA renor- 
malization of EPI in strongly correlated systems which 
was previous ly st udied in connection with heavy fermions 
- see review [l67| . 



D. Charge vertex and the EPI coupling 

The charge vertex function 7^(k, q) (in the adiabatic 
approximation) has been calculated in |5l|, [9^, in 
the framework of the 1/N expansion in the X-method - 
see also and here we discuss only the main results. 
Note, that 7^(k, q) renormalizes all charge fluctuation 
processes, such as the EPI interaction, the long range 
Coulomb interaction, nonmagnetic impurity scattering 
etc. In fact 7^(k, q) describes specific screening due to 
the vanishing of doubly occupancy in strongly correlated 
systems. Note, the latter constraint is at present impossi- 
ble to incorporate into the LDA — DFT band structure 
calculations, thus making the latter method unreliable 
in highly correlated systems. In [5l|, [o^] 7c (k, q, w) 
was calculated as a function of the model parameters 
t, t/, 5, J in leading 0(1) order of the t — J model 



7,(k,g) = l-^^^^„(k)[i 

a=l /3 = 1 



x{Q)VpX^2{'l), (94) 



where Xapil) = Ep ^'^(p, g)^>(p), F„(k) = 
[t(k), 1, 2 Jo cos kx, 2 Jo sin k^, 2 Jq cos fcj,, 2 Jo sin ky], and 
Ga{p,q) = [l,t{p + q),coapa;,smpa:,cospy,sinpy]Il{p,q). 
Here, Il{k,q) = -g{k)g{k + q) and q = (q,ig„), 
qn = 2imT, p = {p,ipm), Pm = 7rr(2m + 1). The 
physical meaning of the vertex function 7^(k, g) is 
following - in the presence of an external or internal 
charge perturbation there is a screening due to the 



J=0.2, t'=0. 



5=0.50 




FIG. 36. Adiabatic (w = 0) vertex function 7(kF,q) of the 
t-J model as a function of the momentum aq with q = (q, q) 
for three different doping levels 5. From [9^. 



change of the excitation potential Vfgg{x,y), i.e. of 
the change of the band width, as well as of the local 
chemical potential. The central result is that for 
momenta k laying at (and near) the Fermi surface the 
vertex function 7^(k, q, w = 0) has very pronounced 
forward scattering peak (at q = 0) at very low doping 
concentration S{<^ 1), while the backward scattering 
is substantially suppressed, as it is seen in Fig. 1361 
where 7^(ki?,q, w = 0) is shown. The peak at q = is 
very narrow at very small doping since its width qc is 
proportional to doping 5, i.e. qc ^ 5{ti / a) where a is the 
lattice constant. It is interesting that 7c(k, g), as well as 
the dynamics of charge fluctuations, depend only weakly 
on the exchange energy J and are mainly dominated by 
the constraint of having no doubly occupancy of sites, 
as it is shown in 
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The existence of the forward scattering peak in 7c(k, q) 
at g = is confirmed by numerical calculations in the 
Hub bard model, which is very pronounced at large U 
[l68j |. This is important result since it means that 
expansion in the X-method is reliable method in study- 
ing charge fluctuation processes in strongly correlated 
systems. The strong suppression of 7^(k, g) at large 
q(~ kp) means that the charge fluctuations are strongly 
suppressed (correlated) at small distances. Such a be- 
havior of the vertex function means that a quasi-particle 
moving in the strongly correlated medium digs up a gi- 
ant correlation hole with the radius ^cft,(~ '"'/'Zc) ~ a/ 5, 
where a is the lattice constant. As a consequence of this 
effect the renormalized EPI becomes long ranged which 
is contrary to the weakly correlated systems where it is 
short ranged. 

By knowing 7c(k, q) one can calculate the relative cou- 
pling constants Ai = A^, A3 = A^, Ktr etc. In the ab- 
sence of correlations and for an isotropic band one has 
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FIG. 37. Normalized s-wave As, d-wave Ad and transport 
Air ■ S coupling constants as a function of doping S for = 
and J = 0. From 111. 



Ai = Atr = 1, Ai = for i > 1. The av erag es in A^, 
Ad and Atr were performed numerically in [96| by using 
the realistic anisotropic band dispersion in the t—tf—J 
model and the results are shown in Fig. [37l For con- 
venience, the three curves are multiplied with a com- 
mon factor so that As approaches 1 in the empty-band 
limit (5 — >■ 1, when strong correlations are absent. Note, 
that the superconducting critical temperature Tc in the 
weak coupling limit and in the i — th channel scales like 
~ exp(-l/(AoA, - where Aq is some effective 
coupling which depends on microscopic details and is 
the effective residual Coulomb repulsion in the i-th chan- 
nel. We stress here several interesting results which come 
out from the theory and can be also extracted with the 
help of Figs. 

(1) In principle the bare EPI coupling constant 
(k, q) depends on the quasi-particle momentum k and 
the transfer momentum q. In the t — J model the EPI 
coupling is dominated by the ionic coupling Hl°" (see 
Eg. ([77)1 ) and corresponding EPI depends only on the mo- 
mentum transfer q, i.e. g^(k, q) « 5°(q) while for the 
much smaller covalent coupling if™" depends on both 
k and q 0, 0- However, the EPI coupling for most 
phonon modes are renormalized by the charge vertex and 
since the latter is peaked at small momentum transfer 
q = |k — k'l one expects that the maxima of the corre- 
sponding effective potentials are pushed toward smaller 
values of q. This vertex renormalization is important 
result since if in the absence of strong correlations the 
bare EPI coupling |(7'^(k, q)| for some phonon modes 
(which enters in the effective t — J model) is detrimen- 
tal for d-wave pairing it can be less detrimental or even 
supports it in the presence of strong correlations. To 
illustrate this let us consider the in-plane oxygen breath- 



ing mode with the frequency ujbr which is supposed to 
be important in HTSC cuprates. The bare coupling 
constant (squared) for this mode is approximately given 
by |gO^(k,q)|^ = \g°^\^ [sm"^ (q^a / 2) + sm'^{qy a/ 2)] which 
reaches maximum for large q = (7r/a,7r/a). By extract- 
ing the component in the d-channel one has 



with 
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stant 



9br 



[1 - (l/4)^,(k)V',(k') ...] (95) 



■!/'c;(k) = COS k^a 



cos kya. 



(96) 



gives 
in 



the 
the 



repulsive 
d-channel. 



coupling con- 



i.e. 



(2M,)^Vd(k)|ff,°,(k-k')| ^rf(k')) < 0. However, in 

the presence of strong correlations one expects that the 
effective coupling constant is given approximately by 

g,t^(k,k-k')| « |5t(k-k')|S2(k^,k- k') which 
is at small doping 5 suppressed substantially at large 
q since 7^ falls off drastically dX q ^ qc ^ 5{n/a). The 
latter property makes the effective coupling constant 
(in the d-channel) A^j-'^-'^for these modes less negative or 
even positive (depending on the ratio ^^/i/" ^ i-^- 
one has X'^J^ > X^. We stress again that this analysis is 
only qualitative (and semi-quantitative) since it is based 
on the t — J model while the better quantitative results 
are expected in the strongly correlated {Ud ^ t, Aj^) 
three-band Emery model - see Appendix D in [Sj]. 
Unfortunately, these calculations are not done until now. 

(2) In weakly correlated systems (or for instance in the 
empty-band limit (5 — > 1) the relative d — wave coupling 
constant A^ is much smaller than the s — wave coupling 
constant A^, i.e. A^ ^ Ag as it is seen in Fig. |371 
Furthermore, Ag decreases with decreasing doping. 

(3) It is indicative that independently on the value 
of i' ^ or = the coupling constant and A^ 
meet each other (note A^ > A^ for all 5) at some small 
doping 5 ~ 0.1 — 0.2 where A^ « A^. We would like 
to stress that such an unique situation (with A^ sa A^) 
was practically never realized in low temperature and 
weakly correlated superconductors and in that respect 
the strong momentum dependent EPI in HTSC cuprates 
is an exclusive phenomenon. 

(4) By taking into account residual Coulomb repul- 
sion of quasi-particles then the s — wave superconduc- 
tivity (which is governed by A^) is suppressed, while the 
d — wave superconductivity (which is governed by A^) 
stays almost unaffected, since /.t* ^ fj,'^. In that case the 
d — wave superconductivity which is mainly governed by 
EPI becomes more stable than the s — wave one at suf- 
ficiently low doping d. This transition between s- and 
d-wave superconductivity is triggered by electronic cor- 
relations because in the model calculations 51 [, [1^ the 
bare EPI coupling is assumed to be momentum inde- 
pendent, i.e. the bare coupling constant contains the 
s — wave symmetry only. 
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(5) The calculations of the charge vertex 7^ are per- 
formed in the adiabatic limit, i.e. for w < q • (q) 
the frequency w in 7^ can be neglected. In the non- 
adiabatic regime, i.e. for cj > q • (q) , the function 
7^(ki?, q,iuj) may be substantially larger compared to the 
adiabatic case because 7^(k_F,q,w) tends to the bare 
value 1 for q = 0. This means that EPI for different 
phonons (with different energies lo) is differently affected 
by strong correlations. For a given u! the EPI coupling 
to those phonons with momenta q < q^^ — cu/vp will 
be (relatively) enhanced since j^{'kp,c[,uj) «1, while the 
coupling to those with q > = uj/vp will be substan- 
tially reduced due to the suppression of the backward 
scattering by strong correlations [l^. These results are 
a consequence of the Ward identities and generally hold 
in the Landau- Fermi liquid theory ^169i] . 

(6) The transport EPI coupling constant Atr is sig- 
nificantly reduced in the presence of strong correlations 
especially for low doping where A^^ < A/3. This mathe- 
matical result is physically plausible since the resistivity 
is dominated by the backward scattering processes (large 
q ^ kp) which are suppressed by strong correlations - the 
suppression of 7^(ki;', q,u;) at large q. 

The theory based on the forward scattering peak in 
EPI is a good candidate to explain the linear tempera- 
ture behavior of the resistivity down to very small tem- 
peratures T(~ 8d/30) « 10 K in some cuprates with 
low Tc(« 10 ii') [3, [Hi. A physically rather plausi- 
ble model (based on the forward scattering peak in EPI) 
is elaborated in It takes into account: (i) the quasi- 
particle scattering on acoustic (a) optic (0) phonons; (ii) 
the extended van Hove singularity in the quasi-particle 
density of states iV(^) which in some cuprates is very near 
the Fermi surface; (iii) umklapp and "undulation" (due 
to the flat regions at the Fermi surface) processes with 
Vk' = — Vk - this condition can partly increase the EPI 
coupling. The transport Eliashberg function a^^F{Ld) 
is calculated similary to Ea. (|9ip by using the definition 
of Q;^i^(k, k', a;) in Ea. (j89p with the renormalized cou- 
pling constant g^y^ (k — k') = (7^(k — k')7j,(k — k') of the 
V — a,o mode, respectively. The phenomenological form 
for the forward scattering peak in 7^ (k — k ) with the 
cut-off qc kp (and which mimics the exact results 
from [m, di], [93) can be found in Q. Since the scat- 
tering of the quasi-particles on phonons (with the sound 
velocity Vs) is limited to small-q transfer processes (with 
q < qc) then the maximal energy of the acoustic branch is 
not the Debye energy Qoi^ kpVg) but the effective De- 
bye energy 6a(~ qcVs) <^ Qd- In the case of Bi — 2201 
in 9^ it is taken (from the numerical results in [5l| , [9^ , 
[93) that qc « kp /lO which gives 6^ ~ (30 - 50) K. As 
the result the calculated a1^F[(jj) gives that Pab{T) ^ T 
down to very low T(- 0.2 6^) « 10 K. The slope 
idp^fy/dT) is governed by the effective EPI coupling con- 
stant for acoustic phonons. In systems with the extended 
van-Hove singularity (in Ni^)) near the Fermi surface, 
what is the case in Bi-2201, the effective coupling con- 
stant for acoustic phonons can be sufficiently large to give 



experimental values for the slope (dp^f^/dT) ^ (0.5 — 1) 
pQcm/K - for details see [9^. This physical picture is 
applicable also to cuprates near and at the optimal dop- 
ing but since in these systems Tc is large the linearity of 
Pab(T) down to very low T is "screened" by the appear- 
ance of superconductivity. 

(7) The forward scattering peak in the charge vertex 
7g(ki?,q) with the width qc ^ Sijr/a) is very narrow in 
underdoped cuprates which may have further interesting 
consequences. For instance, HTSC cuprates are charac- 
terized not only by strong correlations but also by rel- 
atively small Fermi energy Ep, which is in underdoped 
systems not much larger than the characteristic (maxi- 
mal) phonon frequency oj'^^^, i.e. Ep ~ 0.1 — 0.3 eV, 
i^max ^ gg ^ ^j^g appreciable magnitude of 
ijJd/Ep in HTSC oxides it is necessary to correct the 
Migdal-Eliashberg theory by the non-Migdal vertex cor- 
rections due to the EPI. It is well-known that these ver- 
tex corrections lower Tc in systems with isotropic EPI. 
However, the non-Migdal vertex corrections in systems 
with the forward scattering peak in the EPI coupling 
with the cut-off qc « kp may increase Tc which can 
be appreciable. The corresponding calculations |l70l | 
give two interesting results: (i) there is a drastic in- 
crease of Tc by lowering Qc — qd'^.kp, for instance 
Tc(Qc = 0.1) « 4rc(Qc = 1); (m) Even small values of 
Aep < 1 can give large Tc. The latter results open a new 
possibility in reaching high Tc in systems with apprecia- 
ble ratio ud/Ep and with the forward scattering peak. 
The difference between the Migdal-Eliashberg and the 
non-Migdal theory can be explained qualitatively in the 
framework of an approximative McMillan formula for Tc 
(for not too large A) which reads Tc ~ (u;)e~^/[^~^ 1 . The 
Migdal-Eliashberg theory predicts A w A/(l + A) while the 
non-Migdal theory [170] gives A « A(l-I-A). For instance 
Tc ~ 100 K in HTSC oxides can be explained by the 
Migdal-Eliashberg theory for A ^ 2, while in the non- 
Migdal theory much smaller coupling constant is needed, 
i.e. A ~ 0.5. 

(8) The existence of the forward scattering peak in 
EPI can in a plausible way explain the ARPES puzzle 
that the anti-nodal kink is shifted by the maximal super- 
conducting gap Amax while the nodal kink is unshifted. 
The reason is (as explained in Subsection III.C), that the 
EPI spectral function a^F{\i, k', 9) « a^F{ip - ip' , Vl) is 
stron gly peaked due to strong correlations dX Lp — ip' = 

[m. 

(9) Finally, the scattering potential on non-magnetic 
impurities is renormalized by strong correlations giving 
also the for ward scattering peak in scattering potential 
(amplitude) [l20j |. The latter effect gives large d-wave 
channel in the renormalized impurity potential what is 
the reason that the d-wave pairing in HTSC cuprates is 
robust in th e pre sence of non-magnetic impurities (and 
defects) [3, [l20t . 
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E. EPI and strong correlations - other methods 

The calculations of the static (adiabatic) charge- vertex 
7^(k F, q ) in the X-method are done for the case U — oo 
[51| , [9^, [131 where it is found that it is peaked at g = 
- forward scattering peak (FSP). This result is confirmed 
by the numerical Monte-Carlo calculations for the finite- 
U Hubbard model I681] , where it is found that FSP exists 
for all U, but it is especially pronounced in the limit 
t ^ U. The se re sults are additionally confirmed in the 
calculations |l7l| within the four slave-boson method of 
Kotliar-Riickenstein where 7f,(kF,q) is again peaked at 
q = and the peak is also pronounced at t <C C/. 

There are several calculation s of thecharge v ertex 
in the one slave-boson method fl6^ . [l72| . (iTSt . (l73 | 
which is invented to study the limit U — >■ 00. It is inter- 
esting to compare the results for the charge vertex in the 
X-m etho d [5l| . [96j . [97| and in the one slave-boson the- 
ory jl72l | which are calculated in 0{{1/N)'^) order. For 
instance, for J = one has 



1 + 6(q) - a(q)t(k) 
"[l + 5(q)]2-a(qMq) 



(97) 



[l + 5(q)]2-a(q)c(q) 

The explicit expressions for the "bare" susceptibilities 
a(q)i ^(q) a-nd c(q) can be found in [5l|, [9^. It is obvious 
from Eqs.(|n31Ml) that 7^^'(k,q = 0) =7i^'^^(k,q = 0) 
but the calculations give that max{7c (k, q)} is for q — 
0, while max{7i'^^'(k,q)} is for g 7^ [1^. So, the 
SB vertex is peaked at finite q which is in contradiction 
with the numerical Monte Carlo results for the Hubbard 



model |168l | and with the four slave-bosons theory [l7l|. 
The reason for the discrepancy of the one slave-boson 
(SB) in studying EPI with the numerical results and the 
X-method is not clear and might be due to the symmetry 
breaking of the local gauge invariance in leading order of 
the SB theory. 



to the suppression of doubly occupancy at the Cu lattice 
sites in the Cu02 planes, which drastically weakens the 
screening effect in these systems. The pronounced non- 
locality and suppression of the screening is mathemati- 
cally expressed by the charge vertex function 7^(ki?, q, lo) 
which multiplies the bare coupling constants. This func- 
tion is peaked at 5 = and strongly suppressed at large 
q, especially for low (oxygen) hole doping 5 near the 
Mott-Hubbard transition. Such a structure of 7^ gives 
that the d-wave and s-wave coupling constant are of the 
same order of magnitude around and below some opti- 
mal doping 5 op « 0.1, i.e. « Ag. This is very peculiar 
situation never met before. Since the residual effective 
(low-energy) Coulomb interaction is much smaller in the 
d-channel than in the s-channel, i.e. /i* 3> /i^, then the 
critical temperature for c?-wave pairing is much larger 



than for the s-wave one, i.e. ri'*'' ^ Tc*^ Since all charge 
fluctuation processes are modified by strong correlations 
then the quasi-particle scattering on non-magnetic im- 
purities is drastically changed, that the pair-breaking ef- 
fect on c?-wave pairing is drastically reduced. This non- 
local effect, which is not discussed here - see more Q 
and References therein, is one of the main reasons for 
the robustness of d-wave pairing in HTSC oxides in the 
presence of non-magnetic impurities and numerous local 
defects. The development of the forward scattering peak 
in 7^(ki?,q) and suppression at large g(3> (?c = ^{t^Io)) 
gives rise to the suppression of the transport coupling 
constant \tr making it much smaller than the self-energy 
coupling constant A, i.e. one has \tr ~ A/3 near the opti- 
mal doping 5 — 0.1 — 0.2. Thus the behavior of the vertex 
function and the dominance of EPI in quasi-particle scat- 
tering resolves the experimental puzzle that the transport 
and and self-energy coupling constant take very differ- 



ent values, A 



tr^ep 



<^ Aep. This is not the case with the 



SFI mechanism which is dominant at large q thus giving 

Atr,s/ ~ As/. 

We stress that the strength of the EPI coupling con- 
stants Aep, Xep.d is at present impossible to calculate since 
it is difficult to incorporate strong correlations and nu- 
merous structural effects in a tractable microscopic the- 
ory. 



VIII. SUMMARY OF PART II 

The experimental results in HTSC cuprates which are 
exposed in Part I imply that EPI coupling constant is 
large and in order to be conform with d-wave pairing 
this interaction must be very nonlocal (long range), i.e. 
weakly screened and peaked at small transfer momenta. 
In absence of quantitative calculations in the framework 
of the ab initio microscopic many-body theory the effects 
of strong correlations on EPI are studied within the mini- 
mal t— J model where this pronounced non-locality is due 
to two main reasons: (1) strong electronic correlations 
and (2) mixed metallic-ionic layered structure in these 
materials. In the case (1) the pronounced non-locality of 
EPI, which is found in the t — J model system, is due 



IX. DISCUSSIONS AND CONCLUSIONS 

Numerous experimental results related to tunnelling, 
optics, ARPES, inelastic neutron and x-ray scattering 
measurements in HTSC cuprates at and near the optimal 
doping give evidence for strong electron-phonon interac- 
tion (EPI) with the coupling constant 1 < Aep < 3.5. The 
tunnelling measurements furnish evidence for strong EPI 
which give that the peaks in the bosonic spectral function 
a^Fluj) coincide well with the peaks in the phonon den- 
sity of states Fph{uj). The tunnelling spectra show that 
almost all phonons contribute to Tc and that no par- 
ticular phonon mode can be singled out in the spectral 
function Q!^_F(w) as being the only one which dominates 
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in pairing mechanism. In light of these results the small 
oxygen isotope effect in optimally doped systems can be 
partly due this effect, thus not disqualifying the impor- 
tant role of EPI in pairing mechanism. The compatibility 
of the strong EPI with d-wave pairing implies an impor- 
tant constraint on the EPI pairing potential - it must 
be nonlocal, i.e. peaked at small transfer momenta. The 
latter is due to: (a) strong electronic correlations and (b) 
the mixed metallic-ionic structure of these materials. If 
EPI is the main player in pairing in HTSC cuprates then 
this non-locality implies that at and below some optimal 
doping {Sop ^ 0.1) the magnitude of the EPI coupling 
constants in d-wave and s-wave channel must be of the 
same order, i.e. Xep,d ~ Xep,s- This result in conjunc- 
tion with the fact that the residual effective Coulomb 
coupling in d-wave channel is much smaller than in the 
s-wave one, i.e. fi* S> fi^ gives that the critical temper- 
ature for d-wave pairing is much larger than for s-wave 
pairing. 

The tunnelling, ARPES, optic and magnetic neutron 
scattering measurements give sufhcient evidence that the 
spin fluctuation interaction (SFI) plays a secondary role 
in pairing in HTSC cuprates. Especially important evi- 
dence for the smallness of SFI (in pairing) comes from the 
magnetic neutron scattering measurements which show 
that by varying doping slightly around the optimal one 
there is a huge reconstruction of the SFI spectral function 
/mx(Q,aj) (imaginary part of the spin susceptibility) 
while there is very small change in the critical tempera- 
ture Tc- These experimental results imply important con- 
straints on the pairing scenario for systems at and near 
optimal doping: (1) the strength of the d-wave pairing 
potential is provided by EPI (i.e. one has Xep,d ~ Aep,s) 
while the role of SFI, together with the residual Coulomb 
interaction, is to trigger d-wave pairing; (2) the Migdal- 
Eliashberg theory, but with the pronounced momentum 
dependent of EPI, is a rather good starting theory. 

The ab initio microscopic theory of pairing in HTSC 
cuprates fails at present to calculate Tc and to predict 
the magnitude of the d-wave order parameter. From that 
point of view it is hardly to expect a significant improve- 
ment of this situation at least in the near future. How- 
ever, the studies of some minimal (toy) models, such as 
the single band t— J model, allow us to understand part 
of the physics in cuprates on a qualitative and in some 
cases even on a semi-quantitative level. In that respect 
the encouraging results come from the theoretical studies 
of EPI in the t — J model by using controllable mathe- 
matical methods. This theory predicts dressing of quasi- 
particles by strong correlations which dig up a large scale 
correlation hole of the size ^^/i ^ a/S ior S 1. These 
quasi-particles respond to lattice vibrations in such a way 
to produce an effective electron interaction (due to EPI) 
which is long ranged, i.e. the effective pairing potential 
Ve//(q, w) is peaked at small transfer momenta q - for- 
ward scattering peak. This theory (of the toy model) 
is conform with the experimental scenario by predicting 
the following results: (i) the EPI coupling constants in 



d-wave and s-wave channels are of the same order, i.e. 
Xep,d ~ Aep,s, at some optimal doping 6 op ^ 0.1; («i) 
the transport coupling is much smaller than the pairing 
one, i.e. Atr ~ A/3; (Hi) due to strong correlations there 
is forward scattering peak in the potential for scattering 
on non-magnetic impurities, thus making d-wave pairing 
robust in materials with a lot of defects and impurities. 
Applied to HTSC superconductors at and near the op- 
timal doping this theory is a realization of the Migdal- 
Eliashberg theory but with strongly momentum depen- 
dent EPI coupling, which is conform with the proposed 
experimental pairing scenario. This scenario which is also 
realized in the t — J toy model may be useful in making a 
(phenomenological) theory of pairing in cuprates. How- 
ever, all present theories are confronted with the unsolved 
and challenging task - the calculation of Tc- From that 
point of view we do not have at present a proper micro- 
scopic theory of pairing in HTSC cuprates. 
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X. APPENDIX - SPECTRAL FUNCTIONS 

A. Spectral functions o^F(k, k', oj) and a'^F{uj) 

The quasi-particle bosonic (Eliashberg) spectral 
function a^i^(k, k', w) and its Fermi surface average 
a^F(cj) = ((a^F(k, k', cj)))j^ j^, enter the quasi-particle 
self-energy E(k, oj), while the transport spectral func- 
tion a^Ftriui) enters the transport self-energy I](r(k, w) 
and dynamical conductivity cr[uj). Since the Migdal- 
Eliashberg theory for EPI is well defined we define the 
spectral functions for this case and the generalization to 
other electron-boson interaction is straightforward. In 
the superconducting state the Matsubara Green's func- 
tion G(k, a;„) and E(k, a;„) are 2x2 matrices with the 
diagonal elements Gu = G'(k, aj„). En = E(k, a;„) and 
off-diagonal elements G12 = i^(k, aj„), E12 = $(k, aj„) 
which describe superconducting pairing. By defin- 
ing iujn [1 - Z{k,ujn)] = [i;(k,w„) - i;(k, -w„)] /2 and 
X(k,w„) = [I](k,w„) S(k, -w„)] /2, the Eliashberg 
functions for EPI in the presence of the Coulomb inter- 
action (in the singlet pairing channel) read [i^, [175| . 
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Nifi) 



D{p, LOr, 



p,m 



(100) 

^(P,^m) 

iV(M) -"Pj i?(p,c.™)' 

(101) 

where -/V(/i) is the density of states at the Fermi sur- 
face, ujn = 7rr(2n + l), $(k,a;„) = Z(k, a;„)A(k, a;„) and 
D = ujI^Z^ + {e- H + xf + (For studies of optical 
properties - see below, it is useful to introduce the renor- 
malized frequency j£l;„(iw„)(= iuJnZ{LUn)) — LOn ~ S(a;„) 
(or its analytical continuation uj{uj) — Z{u})uj — uj — 
5](w)). These equations are supplemented with the elec- 
tron number equation n(/i) is the chemical potential) 



^ _ 2T ^ e(p) - ^1 + x(p, w„) 



L»(p,a;„) 



(102) 



//. Usually in the case of low-temperature superconduc- 
tors (LTS) with s-wave pairing the anisotropy is rather 
small (or in the presence of impurities it is averaged out) 
whic h all o ws an averaging of the Eliashberg equations 

S, [iTi, [iTi 



Z{uJn) 



E 



(106) 



^(CJ„)A(CJ„) = 7rT^[A(w„ - LOm) 



H{uJc)OioJc - km I)] 



A(w„0 



A(w„ - UJra) = 



dv 



+ [Un - Wm)^ 



(107) 



(108) 



Here a'Fiuj) = ((^^^(k, k', l^)))|^_|^, where ((...))k,k' is 
the average over the Fermi surface. The above equations 
can be written on the real axis by the analytical contin- 
uation iLOm — > oj -|- 1(5 where the gap function is complex 
i.e. A(w) = ^r{uj) +iAi{u}). The solution for A(cj) al- 
lows the calculation of the current- voltage characteristic 
I(V) and tunnelling conductance Gisrs(y) — dl^s/dV in 
the superconducting state of the NIS tunnelling junction 
where lNs{y) is given by 



Note that in the case of EPI one has A^p(i^„) — 

^kpi'^n){= Akp(j^«)) (with = TrTn) where Akp(j^,i) 
is defined by 



r vat 
Akp(i^„) - 2 / 
Jo 



2 



F(i.) = iV(M)^|.g— J'i3.(k-p,j.) 



(103) 



(104) 



where Bi^^k. — p;j^) is the phonon spectral function of the 
K-th phonon mode related to the phonon propagator 



-B^{(l,v)dv. 



(105) 



However, in systems with is measured the general- 
ized phonon density of states GPDS{l^){= G{uj)) 
(see Part I Subsection D) defined by G{uj) = 

where ai and A'U are the 
cross section and the mass of the i-th nucleus and Fi {w = 



(i/A^)E« 



is the amplitude-weighted den- 



sity of states. The renormalized coupling constant 



?^t,kp 



?K kp^e ^7) comprises the screening effect due 



to long-range Coulomb interaction (~ - the inverse 
electronic dielectric function) and short-range strong cor- 
relations (~ 7 - the vertex function) - see more in Part 



iV) 



k,p 



7^k,p P 



duj 
2^ 



AN{k,u:)As{p,u: + ey)[/(^) - f{u; + eV)]. (109) 

Here, Ajv,s(k, w) = —2ImGN^s(^,^) are the spectral 
functions of the normal metal and superconductor, re- 
spectively and f{uj) is the Fermi distribution function. 
Since the angular and energy dependence of the tun- 
nelling matrix elements | Tk,p P is practically unim- 
portant for s — wave superconductors, then in that case 
the relative conductance aNs{V) = GNs{V)/Gj\rN{V) is 
proportional to the tunnelling density of states Nxi^^) = 
/As(k,a;)d^fc/(27r)3, i.e. crArs(w) ~ Nt{uj) where 



Nxiio) = Re 



(w-Hz7(cj))2-Z2(cj)A(w)2 



(110) 



Here, Z{uj) = Z{Lu)/ReZ{uj), j{uj) = -f{uj)/ ReZ{uj), 
Z{uj) = ReZ{oj) -|-J7(a;)/w and the quasi-particle scat- 
tering rate in the superconducting state ^g{uj,T) = 
—2Im'E{uj,T) is given by 



7^(w,T) = 27r J dua^F{v)Ns{v + uj){2nB{v) 
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+ nFiv + io) + nFiiy-oj)} + -f'"'P, (111) 

where Ns{u}) ~ Re{u}/{uj'^ — A^(w))^/^ is the quasi- 
particle density of states in the superconducting state, 
nB,F{i^) are Bose and Fermi distribution function respec- 
tively. Since the structure of phononic spectrum is con- 
tained in a'^F{uj), it is reflected on A(a;) for uj > Aq 
(the real gap obtained from Aq = ReA{uj = Aq)) 
which gives the structure in Gs{V) at y = Aq -f uipt- 
On the contrary one can extract the spectral function 
a^F{u) from Gns{V) by the invers ion procedure pro- 
posed by McMillan and Rowell [1], [Ttt} . It turns out 
that in low-temperature superconductors the peaks of 
—(Pl/dV^ at eVi = A + ujph,i correspond to the peak 
positions of a'^F{uj) and F{uj). However, we would like 
to point out that in HTSC cuprates the gap function is 
unconventional and very anisotropic, i.e. A(k, iw„) ^ 
cosfca;a — cos kya. Since in this case the extraction of 
a^i^(k, k', w) is difficult and at present rather unrealis- 
tic task, then an "average" a'^F{uj) is extracted from the 
experimental curve Gs{V). There is belief that it gives 
relevant information on the real spectral function such as 
the energy width of the bosonic spectrum (0 < < Wmax) 
and positions and distributions of peaks due to bosons. It 
turns out that even such an approximate procedure gives 
valuable information in HTSC cuprates - see discussion 
in Section III D. 

Note that in the case when both EPI and spin- 
fluctuation interaction (SFI) are present one should make 
difference between X^piivn) and A^p(iJ^n) defined by 



'^kp(«^n) = As/,kp(il'n) + Aep,kp(Wn), (112) 



Vkp = KpMpiiVn) - As/,kp(ii^n)- (113) 



u! > luq) and then aini{V) is calculated. In the second 
step the functional derivative 6(j{uj)/6a'^F{uj) {ui = eV) 
is found in the presence of a small change of a'^ Fini{uj) 
and then the iterated solution a^F(i)(cj) — a'^Fini{uj) + 
6a'^F{uj) is obtained, where the correction 6a'^F{uj) is 
given by 



Sa^F{u) = / diy[ 



^da^F{iy)^ ' 



Ml- (114) 



The procedure is iterated until a^i^(„)(cLi) and /i^^^^ con- 
verge to cP'FiLu) and /x*which reproduce the experimen- 
tally obtained conductance (J°^^{y). In such a way the 
obtained a^Fiuj) for Pb resembles the phonon density 
of states Fpb{uj), that is obtained from neutron scatter- 
ing measurements. Note that the method depends ex- 
plicitly on II* but on the contrary it requires only data 



on aNs{V) up to the voltage V-m 



ph 



Aq where 



ph 



is the maximum phonon energy (a^i^(w) — for 
uj > uj'^i^^) and Aq is the zero-temperature superconduct- 
ing gap. One pragmatical feature for the interpretation 
of tunnelling spectra (and for obtaining the spectral pair- 
ing function a^F(w)) in LTS and HTSC cuprates is that 
the negative peaks of cPl/dV^ are at the peak positions 
of a'^F{uj) and -^(0;). This feature will be discussed later 
on in relation with experimental situation in cup rates. 

The second method has been invented in 178l | and it 
is based on the combination of the Eliashberg equations 
and dispersion relations for the Greens functions - we call 
it GDS method. First, the tunnelling density of states 
is extracted from the tunnelling conductance in a more 
rigorous way |l79l | 



Nt{V) = 



1 



(7NS{V) 

<jnn{V) <j*(V) 



du 



In absence of EPI \^p{ivn) and A^p(ii'„) differ by sign 

i.e. A^p(«z/„) — ^X^p{ii'n) > since SFI is repulsive 
interaction in the singlet pairing-channel. 



1. Inversion of tunnelling data 

Phonon features in the conductance CNsiV) at eV = 
Ao + ujpfi makes the tunnelling spectroscopy a power- 
ful method in obtaining the Eliashberg spectral function 
a^F^uj). Two methods were used in the past for extract- 
ing a'^F{uj). 

The first method is based on solving the inverse prob- 
lem of the nonlinear Eliashberg equations. Namely, by 
measuring (7Ars(y), one obtains the tunnelling density of 
states Nt{uj){^ (^JLsi^)) and by the inversion procedure 
one gets a-^F{u}) |l77| . In reality the method is based on 
the iteration procedure - the McMillan-Rowell (MR) in- 
version, where in the first step an initial a'^Fini{uj), /i*„j 
and Aini{uj) are inserted into Eliashberg equations (for 
instance Ai„i(a;) = Aq for a; < wq and Ami(w) = for 



da*(u) .^^ 
X [Nt{V- 
du 



u) - Nt{V)] 



(115) 



where (t*{V) = exp{—f3V}(TNN{V) and the constant P 
are obtained from aNN{V) at large biases - see jlTSf . 
Nt{V) under the integral can be replaced by the BCS 
density of states. Since the second method is used in 
extracting a'^F{u!) in a number of LTSC as well as in 
HTSC cuprates - see below, we describe it briefly for 
the case of isotropic EPI at T=0 K. In that case the 
Efiashberg equations are [H], [l7i|, ^7^, [iTSj 



Z{uj)A{lu) 



duj'Re 



A(^') 



- A2(w')]^^^ 



X [K+{ijj, uj') — ijl*9{lUc — uj)] 



1 - Z{uj) = - 



dijj' Re 



Ao 



- A2(w') 



1/2 



(116) 

K_{uj,uj') 
(117) 
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where 



K±{uj,uj')^ dva^F{v){— — 

^ ' Ja„ ^^lj' + 1^ + ^ + 10+ 



± 



Lj' — uj + ly — iO+ 



(118) 



Here ^* is the Coulomb pseudo-potential, the cutoff ujc is 
approximately (5 — 10) uj'^^^, Ao — A(Ao) is the energy 
gap. Now by using the dispersion relati on fo r the matrix 
Greens functions G'(k, w„) one obtains (l78l | 



where S{uj) = w/ [w^ - A'^{uj)] ^'"^ . From Eqs. ()ll6llll7|) 
one obtains 



i?eA(cj) r-^"- , N lml\(u:) 



+ 



ImA{uj) 



du'Nr{^') / d. , 



(120) 

Based on Eqs. (|115m20|) one obtains the scheme for 
extracting a^F{ijj) 

aNsiy),aNN{V) NriV) (121) 



^ /m5(w) ^ A(cj) ^ a2F(cj). (122) 

The advantage in this method is that the explicit knowl- 
edge of ^* is not required ^178i] . However, the integral 
equation for a'^F{uj) is linear Fredholm equation of the 
first kind which is ill-defined - see the discussion in Sec- 
tion II.B.2. 




FIG. 38. (a) Model phonon density of states F{uj) with the 
peak at uq. (b) The real (solid) Ah and imaginary (dashed) 
part A/ of the gap A(a;). (c) The normalized tunnelling den- 
sity of states Nt{i^) /N{0) (solid ) co mpared with the BCS 
density of states (dashed). From [l80l ]. 



It follows that for uj < uo + Aq most phonons have 
higher energies than the energy w of electronic charge 
fluctuations and there is over-screening of this charge by 
ions giving rise to attraction. For w « wq + Aq charge 
fluctuations are in resonance with ion vibrations giving 
rise to the peak in Aniuj). For cjo + Ao < w the ions move 
out of phase with respect to charge fluctuations giving 
rise to repulsion and negative Afl(w). This is shown in 
Fig. [^Slb). The structure in A(cj) is reflected on Nt{lij) 
as shown in Fig. ISSf c) which can be reconstructed from 
the approximate formula for Nt{uj) expanded in powers 

of A/U! 



2. Phonon effects m Nriijj) 

We briefly discuss the physical origin for the phonon 
effects in Nxito) by considering a model with only one 
peak, at wq, in the phonon density of states F{uj) by as- 
suming for simplicity fj,* — and neglecting the weak 
structure in Nt{l^) at nua + Aq, which is du e to the 
nonlinear structure of the Eliashberg equations |l80j |. In 
Figl38]it is seen that the real gap Afl(w) reaches a max- 
imum at uq + Aq then decreases, becomes negative and 
zero, while A/(cj) is peaked slightly beyond ujq -j- Aq that 
is the consequence of the effective electron-electron inter- 
action via phonons. 



N{0) 2 



Afl(w) 



Aiiio) 



(123) 



As Afl(a;) increases 
iVsc5(cj), while for uj 



above Aq this gives Nxi^^j) > 
ujq + Aq the real value Aj^(cj) 



> 



decreases while A/(w) rises and Nt{u}) decreases giving 
rise for iVT(aj) < Nbcs{^)- 



B. Transport spectral function atrF{uj) 

The spectral function a'^^F{uj) enters the dynamical 
conductivity fjij{uj) = a, 6, c axis in HTS systems) 



56 



which generally speaking is a tensor quantity given by 
the following formula 



a.y(a;) = - — 



l'i.{q,k + q) 



X G{k + q)T,{q,k + q)G{q), 



(124) 



where q — (q, i^) and k ^ {h ~ O.lo) and the bare current 
vertex 'f^iq, fc + k = 0) is related to the Fermi velocity 
VF,i, i-e. Jiiq^k + g; k = 0) — VF,i- The vertex function 
Fj{q, k + q) takes into account the renormalization due to 
all sc attering processes responsible for finite conductivity 
[l8l| . In the following we study only the in-plane conduc- 
tivity at k = 0. The latter case is realized due to the long 
penetration depth in HTSC cuprates and the skin depth 
in the normal state are very large. In the EPI theory, 
^j{q, k + q) = icun, i'^n + i^m) is a solution of an 

approximative integral equation written in the symbolic 
form 83] Tj — Vj + Vef fGGTj . The effective potential 
K// (due to EPI) is given by V^fj = Ek I 5^" P ^k, 
where is the phonon Green's function. In such a case 
the Kubo theory predicts ct™*''"(ci;) {i = x, y, z) 



dv[th{- 



2T 



2T 



-thi — )]S"\u;,,.)}, 



(125) 



■ imp 



and 



I](V(w + - Str(i^) 

contribution. In the 
omit the tensor index ii in (7ii{uj). In the presence of 



where S{u!,h') 

imp 



is the impurity contribution. In the following we 



several bosonic scattering processes the transport self- 
energy Titr{(jj) = Re'Str{u}) + iImTitj.(ui) is given by 

/"OO 

Str(w) = - ^ / dval^iFi{v)[K^{u,v) + iK2{Lo,v)], 

(126) 

r ,1 Ld -\- U , ,1 LU — V , 

Xi(^, v) = i?e[vl'(- + z^-^) - vl,(- + (127) 



-)] (128) 



2T 

Here a1^iFi{v) is the transport spectral function which 
measures the strength of the ^-th (bosonic) scattering 
process and ^' is the di-gamma function. The index I 
enumerates EPI, charge and spin-fluctuation scattering 
processes. Like in the case of EPI, the transport bosonic 
spectral function a^^iFiyi) defined in Eq. (|9ip is given 
explicitly by 



iV2 



F(w). 



_F.k'^F,k 



(129) 



We stress that in the phenomenological SFI theory Q 



one assumes a 



kk 



,F{u;) 



7V(/i)g2^Imx(k 



p, ui), which, 



as we have repeated several times, can be justified only 
for small gsf, i.e. gsj <^ Wh (the band width). 

In case of weak coupling (A < 1), (7(w) can be written 
in the generalized (extended) Drude form as discussed in 
Section III.B. 
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